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Chapter I 
INTRODUCTION 


Switching-mode dc-dc regulators are coming into increasing 
use as power supplies because of the significant reduction 
of weight, size and increase in equipment efficiency that can 
be attained. The dissipative regulators used earlier had 
the advantage of design simplicity, but suffer from low ef- 
ficiency and higher weight due to the large transformer and 
filter needed to achieve the required voltage/current regu- 
lation. 

Since the late sixties, the rapid expansion of the com- 
puter and communication industries and increasing complexity 
and sophistication of various systems necessitated the de- 
velopment of higher performance switched mode power sup- 
plies. The incentive for performance improvement prompted 
the initial development of multiple loop control schemes, 
such as the standardized control module (SCM) for dc-dc con- 
verters [2,7,8] and the current-injected control scheme [9]. 
The ever continuing search for performance improvement forms 
the underlying theme of this dissertation. 

The work presented in this dissertation is mainly con- 
cerned with developing a control scheme to alleviate the 
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problem brought about due to the use of an input filter in 
switching regulators. The switching regulator input current 
has a substantial ripple component at the switching frequen- 
cy end this necessitates the use of an input filter to 
smooth the pulsating current drawn from the supply. Furth- 
ermore the input filter also serves to attenuate noise pre- 
sent in. the supply voltage. from being propagated through the 
regulator to the payload, downstream. The presence of the 
input filter, however, often results in various performance 
difficulties such as loop instability, degradation of tran- 
sient response, audiosusceptibility (closed loop input-to- 
output gain) and output impedance characteristics 
(3, 4, S, 6,101. These problems are caused mainly by the inter- 
action between the resonant peaking of the output impedance 
of the input filter amd the regulator control loop. Conven- 
tional single-stage and two-stage input filters can be de- 
signed such that the peaking effect is minimized, however 
such a design is often accompanied with a penalty of weight 
and loss in the input filter [3,4,5,6,10]. This disserta- 
tion presents a different approach via a feedforward control 
scheme to mitigate the undesirable interaction between the 
input filter and the regulator control loop. 

The concept of pole-zero cancellation is used, in this 
dissertation, to develop a novel feedforward control loop 
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that senses the input filter state variables and processes 
this information in a manner designed to csuicel the detri- 
mental effect of peaking of the output impedance of the in- 
put filter. The feedforward loop working in conjxinction 
with the feedback loops developed earlier [2,7,8] constitute 
a total state control scheme that eliminates the interaction 
between the input filter and the regulator control loop. Em- 
ploying the novel feedforward compensation scheme presented 
in this dissertation, a high performance converter together 
with an effective input filter design (minimum weight and 
loss) can be accomplished concurrently. A buck regulator 
employing a feedforward control loop working in conjionction 
with the feedback loops was used to obtain measurements that 
showed significant improvement in the following performance 
categories ; 

1. Loop stability (open loop gain and phase margins); 

2. Audio susceptibility; 

3. Output impedance; and 

4. Transient response. 

In this dissertation the problem caused by input filter 
interaction and conventional input filter design techniques 
are discussed in Chapter 2 followed in Chapter 3 by the con- 
cept of pole-zero cancellation developed earlier [2,7,8]. 
Chapters 4 and 5 present the modeling of the power stage 
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with input filter and the implementation of the feedforward 
for a buck regulator respectively. Measurements of open loop 
gain and phase that confirm the analytical prediction of 
performance improvement using the feedforward scheme devel- 
oped, are discussed in Chapter 6 along with other measure- 
ments of audiosusceptibility, output impedance and transient 
response. Chapter 7 presents experimental and analytical re- 
sults pertaining to transient response while Chapter 8 dis- 
cusses the use of the feedforward loop in stabilizing a re- 
gulator system made xinstable due to input filter 
interaction. Chapter 9 is concerned with extending the con- 
cept of feedforward compensation to other types of control 
and to other types of regulators i.e. the buck-boost regula- 
tors. Finally, Chapter 10 presents the conclusions and sug- 
gestions for future work. 



Chapter II 


INPUT FILTER RELATED PROBLEMS AND CONVENTIONAL 
DESIGN TECHNIQUES 

2.1 INPUT FILTER RELATED PROBLEMS 

An input filter is often required between a switching regu- 
lator and its power source. A buck type switching regulator 
with a single-stage input filter is shown in Figure 1 . The 
regulator input current has a substantial pulsating current 
component at the switching frequency as a result of the 
opening and closing of the switch and this component should 
be prevented from being reflected back into the source ; an 
input filter is required to provide high attenuation at 
switching frequency and thus smooth the current drawn from 
the source. The input filter also serves to isolate source 
voltage disturbances from being propagated to the switching 
regulator payload downstream. 

A presumably well-designed input filter, satisfying the 
above mentioned requirements, when used with a switching re- 
gulator can often cause significant performance degradations 
[3,4,5,6,10]. This is due primarily to the complex interac- 
tion between the switching regulator control loop, the input 
filter and the regulator output filter [3,4,5,6,10]. 
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DUTY CYCUE 



Figure 1: Buck Converter with an Input Filter 


V( =12v Vq=5v Rl= 0.86 ohm L=200 microH 

R^=20mohm C=1540microF Rc=7mohm Tp=50microsecs 
L^=77 microH R =39.6 mohm =412 microF 
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The interaction between the control Loop and the input 
filter is illustrated in Figure 2 . The switching regulator 
has been shown to have a nonlinear negative resistance, as 
illustrated in the figure, [3]. The input current ij. to the 
switching regulator is related nonlinearly to the input vol- 

<iir 1 

tage e- and the input resistance --- = . Under certain 
conditions . the input filter-switching regulator combination 
can become a negative resistance oscillator, producing large 
amplitude voltage excursions across capacitor C. When this 
happens serious degradation of regulator performance could 
occur, [3], including loss of stability. 

The, effect of the input filter is more clearly seen us- 
ing a small signal model. 

The averaging technique [1] is used to relate the low 
frequency modulation component of the source voltage and 
control signal to the corresponding frequency components of 
the converter output voltage. Using the continuous inductor 
current buck regulator with input filter of Figure 1 , as an 
example, a small signal model using the dual-input describ- 
ing function can be developed, as shown in Figure 3 , [6]. 
In this model the effect of the input filter is character- 
ized by the following two parameters: the forward transfer 
characteristic of the input filter H(s) and the output impe- 
dance of the input filter 2(s). 




Figure 2: Negative resistance oscillation. 



Figure 3 ; 


Small-signal mode 
describing functi 
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Fp(S) 



using the dual -input 
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In Figure 1 the output filter is made up of R^, L, and C, 
Rj^ is the load resistance, D is the steady state duty cycle 
ratio D = Tq^/T, Vj and Ij are the steady state regulator 
input voltage and current respectively, and the lower case 
letters with a caret cibove them denote modulation signals. 

The small signal model of Figure 3 is used to illus- 
trate briefly the complex interaction between the input fil- 
ter, output filter and the control loop and the problems 
caused by the interaction. For detailed auialysis please re- 
fer to [ 6 ] , 

2.1.1 Input Filter Interaction -- Loop Stability and 
Transient Response 

The stability of a switching regulator can be examined by 
the open loop gain G^(s): 

0 ^( 3 ) = Fc(3)Fp{s)Fj;{s)F5j(s) (2-1) 

where F(j(s)Fp(s) is the duty cycle- to-output describing 
fiinction v^/d , and Fg(s), F 5 j(s) are the transfer functions 
of the error processor and the pulse modulator respectively. 
The peaking of the output impedance of the input filter Z(s) 
has the following effects:' 

(1) The duty-cycle power stage gain F^,(s) includes the 

output impedance Z(s) — 


F^,(s) = Vj-Z(s)Ij or 


( 2 - 2 ) 
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F^(s) 



(2-3) 


h 


j is the negative input 


The first term in the brackets 
impedance of the regulator. At the input filter resonant 
frequency, Z(s) reaches a peak value and if this value is 
large enough the result could be a reduction in loop gain or 
even worse a negative duty cycle power stage gain F^j(s). Re- 
duction in loop gain could lead to loop instability, whereas 
a negative Fq(s) together with the negative feedback loop 
will result in a positive feedback unstable system. 

(2) The power stage transfer function Fp(s) includes the 
output impedance 2(s) — 


[R- + l/sC]//K. 

F (s) = — = where (2-4) 

P D'^Z(s) + Zj^(s) 


Zj^(s) = Rj^ + sL + [R^, + 1 /sC]//Rl (2-S) 

= input impedance of the regulator. 

Excessive Z(s) at the input filter resonant frequency can 
significantly reduce Fp(s), and thus the loop gain. 

Figures 4 [6], illustrate the effect of peaking of Z(s) 
on the duty cycle- to-output transfer function F(^(s)Fp(s) if 
an improperly designed input filter is employed. 
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<tb 




Figure 4: Duty cycle-to-output voltage characteristic with 

and without input filter (a) gain (b) phase, 
(parameter values are as shown in Fig.l) 
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At the input filter resonant frequency the peaking of the 
output impedance Z(s) causes a sharp change in the gain and 
phase of the duty cycle-to-output transfer ftinction. This 
could result in loop instability and degradation of tran- 
sient response from a presximably well damped system to an 
oscillatory one; control of the peaking effect of the output 
impedance Z(s) is necessary to avoid these problems. 

2.1.2 Input Filter Interaction — Audiosusceptibility and 
Output Impedance 

The audiosusceptibility refers to the switching regula- 
tor's ability to attenuate small signal sinusoidal distur- 
bances present at the input so as not to affect the regulat- 
ed output voltage. The audiosusceptibility performance is of 
considerable importance, as the regulator generally shares 
the input bus with other on-line equipment. The operation of 
this equipment generates noise voltages on the input line 
which must be attenuated by the closed- loop regulator so 
that operation of the various payloads at the regulator out- 
put will not be adversely affected. The audiosusceptibility 
is expressed in terras of the closed loop input- to-output 
transfer function G (s): 

G^(S) 


_ Vfi 


V^ (s) 


F_(s)F„(s) 

I ^ 


l+Fc(s)Fp(s)FE.(s)Fj^(s) 


(2-6) 

1 + G^(s) 
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where Fj(s) = DH(s) = input voltage gain of the power stage. 
G (s) and thus the audiosusceptibility are affected by the 
reson 2 Uit peaking of the output impedance Z(s) and of the 
forward transfer function of the input filter with the regu- 
lator disconnected H(s), because Fj(s) is a function of H(s) 
whereas Fp{s) and Fp(s) are fxinctions of Z(s). The reduc- 
tion of loop gain at the resonant frequency can thus severe- 
ly degrade the audiosusceptibility . Figure 5 , [6], illus- 

trates the audiosusceptibility of the buck regulator with 
euad without an input filter. 

The output impedance of the regulator should be small 
so that the regulator behaves like an ideal voltage source, 

however the output impedance is increased by the peaking of 

/ • 

the output impedance of the input filter. 


2o(s) 


Zp(s) 

1 + G^(s) 


(2-7) 


where Zp(s) is the output impedance of the power stage with 
the control loop open. 

At the resonant frequency the output impedance of the regu- 
lator Zq(s) is increased. This is a consequence of the loss' 
of loop gain G^(s) as a result of. peaking. 
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db 



Figure 5: Audiosusceptibility with and without input 

filter. 

(parameters values are as shown in Fig.l) 
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The peaking of H(s) and 2(s) of the input filter thus 
results in a reduction in the loop gain, this in turn af- 
fects stability, transient response, audiosusceptibility and 
the output impedance of the regulator - 

A buck type switching regulator with a two stage input 
filter is shown in Figure 6 . Figure 7 , [6], shows the mea- 
sured values of open loop gain and phase as a function of 
the frequency. ( The input filter damping resistance is 
not employed to purposely illustrate the effect of input 
filter interaction with the regulator control loop. ) Signi- 
ficant changes in the open loop gain and phase characteris- 
tics at the resonant frequencies of both the first stage and 
the second stage of the input filter are observed, [6]. Thus 
it is seen that the peaking of the output impedance at reso- 
nant frequency of the two stage input filter can also cause 
serious performance degradation. 

2.2 INPUT FILTER DESIGN CONSIDERATIONS 

The design of the input filter is made more complicated by 
the necessity of satisfying the following constraints, which 
result from the interaction between the input filter and the 
regulator control loop discussed in section 2.1 : 

1. The amount of regulator switching current reflected 
back into the source should be limited (conducted in- 
terference requirement) . 



input filter. 


r Vo=20v Rl=7 o] 
I C=455microF 
. mH Rj =0 . 05 mohm 

)2 ohm Cj =225 mien 
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Figure 7; Open loop gain and phase characteristics without 
the input filter damping resistance 

(parameter values used are as shown in Fig. 6) 


PHASE » .(f GJS) 
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2. The peaking of the output impedance of the input fil- 
ter Z(s) should be limited to a safe value to avoid 
significant loop gain reduction. 

3 . The peaking ' of the transfer function of the input 
filter H(s) should be limited to achieve a satisfac- 
tory rejection rate of audio signals propagating from 
input to output. 

4. The input filter weight and energy loss should be 
limited to low values. 

5. The Nyquist stability criterion has to be satisfied; 
thus the closed loop poles should be in the left half 
plane for stable operation - 

|1 + F^,(s)Fp(s)Fjj(s)Fj^(s)| > 0 

6. The closed-loop input- to-output transfer characteris- 
tic (audiosusceptibility) and transient response due 
to a sudden line/load chamge should not be degraded 
by a notice2j3le amount. 

An input filter design that satisfies one constraint 
may often result in violating some other constraint. For ex- 
ample, an input filter design that limits performance degra- 
dation (degradation of stability, transient response and au- 
diosusceptibility) often results in higher weight and 
increased losses in the input filter. A satisactory input 
filter design trades off one or more of the performance deg- 



20 


radations for size, weight and loss. A near-to- optimal de- 
sign thus requires many trial and error design attempts. 
Some of the conventional input filter design techniques are 
presented next. 

2.2.1 Conventional Incut Filter Design Techniques 
A single stage input filter as shown in Figure 8(a) can be 
designed to avoid performance degradation — but this would 
result in larger filter and thus resulting in weight 
and size increase. The filter is simple and commonly used 
but it cannot often satisfy the stringent requirement on au- 
dio susceptibility without size/weight penalty. Resonant 
peaking of the filter of Figure 8(b), [3], is lowered by ad- 
ding resistance R, but this lowers efficiency because the 
pulse current flowing through increases losses. Another 
design uses a resistance R in parallel across , [4], but 

this results in a large Cj . 

The optimal design of a single stage input filter thus 
is rather difficult without tradeoff between performance 
degr adat i ons ahdTthe weight and loss limitations. 

The degradation of the power stage transfer function 
Fp(s) due to peaking of Z(s) can be avoided if there is suf- 
ficient separation of the input filter resonance frequency 
^ and the output filter resonant frecuency 
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cOq = —r , [4,5,10]. Figure 9 , [4,5,10], shows three pos- 
fLC 

sible -combinations of (Jq and . ^p(s) is related to both 

Z(s) and the input impedance of the regulator Z£(3) thus 

[R^ + l/sC]//R- 

F (s) = — i i (2-8) 

P D Z(s) + 2^(s) 

where Zj^(s) = R^ sL + l/sC]//R^ 

— input impedance of the regulator. 

Z(s) and Z^(s) peak at the frequencies cjq snd respective- 
ly. If the two resonance frequencies are the same as in Fig- 
ure 9 then both Z(s) auid Z^(s) peak at the Seune frequency 
and thus at that frequency the transfer fvmction Fp(s) would 
be affected, i.e. reduced, to the maximum possible extent. 
Shifting -the two frequencies cJq and apart as shown in 
Figure 9 will result in reducing -the effect of peaking on 
Fp(s). Reducing 4J^ would result in increasing the size and 
weight of the input filter. A high value of is desirable 
from the point of view of weight and size reduction but this 
can re sult in __ severe perf Qrmance_degradati on_ -- _ f rom^ -Figure- 
4 it is clear that the gain of the duty cycle- to -output 
transfer function F^(s)Fp(s) decreases with increasing fre- 
quency and thus the effect of peaking of Z(s) on the gain of 
Fq(s)Fp(s) would be more pronounced if Z{s) peaks at a high- 

1 * 


er ui 
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The reduction of the loop gain at higher input filter 
resonant frequency cjj often results in poor audio suscepti- 
bility, oscillatory transient response or even an unstable 
system. The choice of thus involves a trade off between 
meeting performance specifications and size/weight. 

2.2.2 Optimal Conf i qur ati on 

A two-stage input filter configuration has been described, 
[3,6] and is shown in Figure 10 . The first stage consisting 
of Lj,Cj,R 3 and Rj controls the resonant peaking of the fil- 
ter. The second stage consisting of L^, supplies most of 
the pulse current required by the regulator. As shown in the 
literature, [6], the two-stage input filter is capable of 
reducing H(s) and 2(s) at resonant frequency without signi- 
ficantly increasing weight and loss, unlike the single-stage 
input filter. Computer optimization techniques have .been 
utilized to optimally design the two-stage filter[6j. It has 
been shown that the two- stage filter is much lighter than 
its single-stage counterpart iinder identical design const- 
“■ra’int^ AI^^Yt'haa been shown that for the same filter 
weight the single stage filter has a significantly higher 
peaking of H(s) and 2(s). Figure 11 shows the gain and phase 
of the duty cycle-to-output describing function of a power 
stage with a two stage input filter, [6j. 


Figure 10: 
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REGULATOR 


Two- stage input filter. 
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db 





a» 

cJ> 


Figure 11: Duty cycle- to-output describing function of 

power stage with two- stage input filter. 

(parameter values used are as shown in Fig.l with the fol 
lowing two-stage input filter parameters -- 
L,=232 microH R^=0;0276 ohm '^^,=11 microH 

=0.0119 ohm C =100 microF R2=l-73 ohm 

02^=30 microF ) 
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Figure 4 shows the gain and phase of the duty 
cycle- to -output describing function of a power stage with a 
single-stage input filter, and the two-stage filter of Fig- 
ure 11 was designed to have the same weight as the single- 
stage input filter of Figure 4 . Comparing the two figures 
the improvement in performance regarding the duty cycle-to- 
output transfer function is dramatic. 

It can therefore be concluded that the two-stage filter 
provides the best compromise among the conflicting require- 
ments of an input filter. 

2.2.3 Input Filter Compensation Via a Feedforward Loop 
Limiting interaction between the input filter and the regu- 
lator control loop is possible with the addition of a feed- 
forward control loop. At- this point it is important to em- 
phasize that such a scheme is designed to eliminate the 
effect of peaking of the output impedance of the input fil- 
ter Z(s), since the peaking of Z(s) interacts with the con- 
trol loop. The forward transfer function H(s) does not in- 
teract with the regulator control loop, as is evident from 
Figure 3 and therefore the peaking of H(s) cannot be cont- 
rolled in any way by adding a feedforward loop. The peaking 
of H(s) can only be controlled by proper filter design. In 
this dissertation a feedforward loop is implemented for a 
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switching buck regulator employing multiple control loops. 
The feedforward loop will be designed to cancel the detri- 
mental effect of input filter interaction and thus improve 
the regulator performance and stability. The use of the 
feedforward control loop thus removes some of the conflict- 
ing design constraints mentioned above and makes an optimal 
input filter design more easily attainable. 

2.2.4 Oblectives of Feedforward Loop Design 

The proposed feedforward loop will be designed such that : 

1. It will eliminate input filter interaction with the 
regulator loop. This will result in improvement in 
stability margins, audiosusceptibility, output impe- 
dance and transient response of the regulator. 

2. It will allow the input filter to be optimized. Some 
of the design constraints that make an optimal filter 
design difficult to attain are removed with the addi- 
tion of feedforward. 

3. It will eliminate equipment interaction. Figure 12 
shows a switching regulator and _ijt^ pr^egul^tor 
which may be a rectifier and a filter. The dynamic 
output impedance of the preregulator will interact 
with the switching regulator and may cause problems 
like loop instability, degradation of audiosuscepti- 
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bility, output impedance and transient response. The 
addition of a feedforward loop will eliminate such 
interaction, thus isolating the switching regulator 
from equipment upstream. 
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Preregulator 


Figure 12: 


Switching 

Regulator 


Interaction between the switching 
regulator and its preregulator. 




Chapter III 


CONCEPT OF POLE- ZERO CANCELLATION USED IN TOTAL 

STATE CONTROL 


3.1 CONCEPT OF POLE- ZERO CANCELLATION . 

The concept of pole-zero cancellation that was developed 
earlier [2,7,8] was implemented by feedback loops that sense 
the regulator output filter state variables and process this 
information to achieve better performance, and a control 
adaptive to filter parameter and load changes. 

Figure 13 [2], shows a two-loop controlled switching 
buck regulator. The dc loop senses the converter output vol- 
tage and compares ii with the reference voltage to generate 
a dc error signal for voltage regulation. The ac loop senses 
the ac voltage across the output filter inductor to generate 
an ac signal. Both ac and dc signals are processed through 
an operational amplifier summing junction to provide a total 
error signal at the output of the operational amplifier 
integrator. It is apparent that the error signal at the out- 
put of the integrator contains information regarding the 
output filter state variables -- the inductor current and 
the capacitor voltage. 
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Fi^re 13: Two loop controlled switching buck regxilator 
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It was shown, [2], that the feedback control loops when 
properly designed can provide complex zeros to cancel com- 
pletely the complex poles presented by the low-pass output 
filter of the power stage; It was also shown that the feed- 
back control loop has the ability to sense filter parameter 
cheuages and automatically provide pole-zero cancellation. To 
examine the adaptive nature of the control loops the open 
loop regulator transfer function G^(s) is used 


G^(s) 


KZ ( jdi) 
sP(]aj) 


(3-1) 


where K is a constant determined by the power stage and con- 
trol loop parameters and 

Z(jw) * 1 + “/“nl * “^'^“^nl (3-2) 

P(j«) = 1 + j2^2 (3-3) 

P(j<o) has complex poles corresponding to the output filter 
and Z(jcj) has complex zeros produced by the two loop feed- 
back control. 


w , = /— 

nl / LC 


0) 


n2 



(3-5) 


(t) 

r - T (3-6) 

'^l 2 ^z 






R,C) 


(3-7) 
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I± 

NR^ 


(3-8) 


T 


Z 


(Ry + ^5)02 + 



(3-9) 


Ry = (R3_/R2) + ^^ 3 ! 


(Rl + R2) 


(3-10) 


L, C, and R^^ form the output filter as in Figure 13 . The 
control parameters can be chosen such that 


“nl “ “n2 


(3-11) 

?! » ?2 


(3-12) 

thus resulting in 

P(ju) = z(jw) 

(3-13) 

and 

St(s) = f 

(3-14) 


The open loop transfer function is of first order and is 
completely independent of output filter parameters. The 
ad'^tiVe nature of the coni^ol loop is apparent from the 
fact that the complex zeros imitate the change in the com- 
plex poles due to component tolerance, aging or temperature 
variations, thus preserving the pole-zero cancellation. 
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3.2 TOTAL STATE CONTROL 

The concept of pole-zero cancellation of the output filter 
characteristics led. to the idea of using similar means to 
control the effect of peaking of the output impedance of the 
input filter. The objective of the work reported in this 
dissertation is thus to develop a feedforward loop that 
senses the input filter state variables and uses the infor- 
mation contained therein to eliminate the interaction bet- 
ween the input filter and the regulator control loop. Such a 
feedforward loop working in conjunction with existing feed- 
back loops forms a total state control scheme, which is il- 
lustrated in Figure 14 . The f eedf orward . loop senses the in- 
put filter state variables and feeds this information to the 
error processor. 

The other inputs to the error processor are the ac vol- 
tage across the output filter inductor and the output vol- 
tage - as shown in section 3-. 1 these contain information 
regarding the output filter state variables. 

The pulse modulator thus has as its input information re- 
garding the state variables of the output filter and also 
the input filter. The duty cycle signal d(t), which controls 
the switch in the power stage, is thus also affected by the 
input filter state variables. It is shown later in this dis- 
sertation, in Chapters 4 and 5, that the feedforward loop 
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Figure 14: Total State Control. 
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can be designed such that interaction between the input fil- 
ter and the regulator is eliminated. 



Chapter 17 

MODELING OF THE POWER STAGE WITH INPUT FILTER 


The first step in the design and analysis of the feed*° 
forward loop is to develop the small signal model of the 
power stage with input filter for the buck-boost, buck and 
boost type of switching regulators, using the averaging 
technique, [1]. The modeling is carried out in the continu- 
ous conduction operating mode, in which the inductor current 
is always nonzero. This mode is the prevalent operating mode 
for most dc-dc converters. The discontinuous conduction op- 
erating mode in which the inductor current is zero for some 
time during the cycle occurs at light loads and is seldom 
used as the intended design at full load. 

The modeling is carried out in the following steps — 

1. State space equation formulation during Tqpj and Tqj,j.. 

2. State space averaging and perturbation. 

3. Linearization and derivation of the small signal 
equation ^r^_-Uie sn^alJ,_signal__,equi valent circuit -and -- 
state space model. 
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4.1 BUCK-BOOST CONVERTER SMALL SIGNAL MODEL DERIVATION 
The buck-boost converter is shown in Figure 15 . In the 

buck-boost converter shown the input filter is composed of 
R^j, LI, Rjjj and Cl. The load is represented by Rj^ while Vj 
and Vq are the input and output voltages respectively. Dur- 
ing Tqj^, the switch S is on and the circuit as shown in Fig- 
ure 16 { a) . 

The equations describing the circuit are 




<p » flux in core 


= .R l£i 

dc LI ^ *X1 ^Cl^ LI Lp " LI 


:ii 

LI 


ii » _£i j. "Cl . _L 

dt Np Np ^1 ’ Lp 


'^Cl . ®C1 


'Cl 


+ Rp) 


(4-1) 


(4-2) 


(4-3) 


^/ii A.. 

dt Cl ■ ClLp ^ 


dt “ C(Rj^+R^) 

Vc • 


(4-4) 


(4-5) 


(4-6) 



Figure IS: Buck-boost converter power stage. 
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During Topp the switch S is off and the circuit is as shown 
in Figure 16 (b) . 

The equations describing the circuit are 


^s- 




(4-7) 


^1 f 

_ _ - "'^Cl 



dt 

“TT < 

■®L1 ■ ®Cl' ' U 

LI 

(4-8) 

ii , 


* w* 

«l"c 


dt 


«c * ®L> 

•>s«c * 


^^^Cl 

^1 




Cl 

Cl 



(4-10) 






dt 

^s^^c 



(4-11) 


Vc 

, Ws* 


(4-12) 

^0 ^ 





The following vectors are defined 


— 

^1 




U - [Vj] 

X » 




'^Cl . 

Z- [Vq] 


(4-13) 


resulting in the following state space equations 
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X = X + u 

Z - X 


where 


T 

OFT 

X » Aj X + u 
Z “ ^2 X 


(4-14) 


C 


*1 


^\l‘‘‘^Cl^ 

\l"p 

-1 

LI 

, LI 

LlLp 

^C1 

-(Rci +Rp) 

1 




1 

-Np 

A 

Cl 

ClLp 

U 

0 

0 

0 


0 

0 

0 

-1 


LI 

0 

1 

Cl 

0 


0 



0 


CLgCRc+V 


2l 

LI 

0 

0 

0 


0 

0 

-1 

C(Rc+V 



0 


B 


2 


B. 



(4-15) 


(4-16) 

(4-17) 

(4-18) 





(4-19) 
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Ws 






( 4 - 20 ) 


^ — 


( 4 - 21 ) 


The state space averaged model over the entire period T is 
X - [dA^ + d'A2] X + [dBj_ + d’B2l u 
X - [dCj^+d’C2] - 

where d = duty cycle ratio = Tojg/ ( T on'^Tq^j) 


( 4 - 22 ) 


d' s 1-d 

T » T-„+T-_ 

ON OFF 

The state space averaged model is perturbed thus - 
d = D + d 
d' = D' - d 
u = U + u 

y = Y + y 
X = X + X 


( 4 - 23 ) 


( 4 - 24 ) 


Assximing that the perturbation is small 
d X 

— << 1, — << 1 etc. leads to the following small 
D X 


signal linearized model 
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0 = [DA^ + D'A2] X + [DBj^ + D'Bj] 

Vq - [DC^ + D'C^] X 

« 

X » [DA^ + D'A2] X + [DB^ + D’B2] (4-25) 

+ [A^ - A 2 ) X + - B 2 Vj] d 

Vq = [Cj^ - C 2 ] X d + [DCj_ + D.'C2] X 

Defining 

A » DA^ + D'A2 

B - DBj^ + D'B2 (4-26) 

C » DCj_ + D'C2 

results in 

X » A X + B Vj 

+ [A^ - A 2 ) X + (Bj^ - B 2 ) Vj] d ' (4-27) 

Vq - [C^ - C 2 ] X d + C X 

Using Laplace transforms results in 

x(s) - [SI - A]"^ Vj(s) 

+ [SI - A]“^[(A^ - A 2 ) X + (B^ - B 2 ) Vj] d(s) 

Vq(s) = [Cj^ - C 2 ] X d(s) + C x(s) 


(4-28) 
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The state space model is derived from the two equations 
above, and is shown in Figure 17 . 

To derive the small signal equivalent circuit, it is 
first necessary to use the equation for Vq to substitute for 

A A 

v_ in terms of v_ , in the. small signal equations described 

w O 

above. Simplifying, the four equations are obtained 




''i ■ «01 - -lP ♦ * 


LI 


dt 


(4-29) 


^Cl^P - 

+ Vci - ^ d 

dv ^ DN . N 


Lp - I, 


(4-30) 




(4-31) 


Ici^ 


Ngd^ RciNgD 
dt 






(4-32) 


— R-N. DD iRpIL-N- ^ 


- d + — 7 d 




h> 


RcW^(P-D’) 



4' 



Figure 17; Buck-boost conve: 
state space model 




^Cl^l^S ^ ^Cl^S 


Figure 18: Buck-boost converter power stage small signal 

equivalent circuit 
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Using a fictitious current i described by 

Lgi = Ng^ (4-33) 

an equivalent circuit can be made up that is described by 
the four equations given above. This circuit will use the 
current i flowing through Lg and it is thus the small signal 
equivalent circuit for the buck-boost converter, as shown in 
Figure 18. 


4.2 BUCK CONVERTER SMALL SIGNAL MODEL DERIVATION 
The procedure used in deriving the small signal model for 
the buck converter is exactly similar to that used for the 
buck-boost converter. The buck converter is shown in Figure 
19 . 

During Tqj^^ the switch S is on and the circuit is as 
shown in Figure 20 (a) . 

The equations describing the circuit are 


*^^1 ^ -(^Ll ^Cl^ ^ ^ ^ 

dt “ LI T.1 LI h. “ LI LI 

^ h ■ +IC1 Vc 

dt * L T.1 " L h. L ” L(R^+R^) 


(4-34) 


(4-35) 


dt “ Cl " Cl 


(4-36) 


^^C ^C 

dt * C(R^ ~ 


(4-37) 
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(4-38) 


During Tqpj,, the switch S is off^ and the circuit is shown 
in Figure 20 (b) , The equations describing the circuit are 


dt " LI Hi “ LI LI 


(4-39) 


fx Vc 

dt " L ■ L(Rj^ + R^O 


(4-40) 


Hi 

dt “ Cl 


(4-41) 


Vl H 

dt “ C(R^ + 


HHH , Vc 
W 

^2 " R^ + 


(4-42) 


(4-43) 


(4-44) 


H “ Hi *'■ ^Z R(, + 


(4-45) 


The following vectors are defined 
,ii = [Vj-1 i£ 


Hi 

H 


(4-46) 


V = [Vq] 





S3 


resulting in the following state space aquation 


ON 


^ X + B^u 

2 ’ C- X 


^OFF 

X = A2 X + B2U 

2- ^2 X 


(4-47) 


where 




LI 

L 

1 

Cl 


31 

LI 

- 8 ., 


Cl 


C(R^ + R^) 


Ll 

1 

L 


0 


L(Rc+V 

0 


(4-48) 




Ll 


_1_ 

Cl 


:h 

L 


2I 

Ll 


0 

0 


C(Rj^ + R^,) 


0 


L(R^+R^) 

0 


-1 


C ( R ^+ Rc ) 


(4-49) 


0 


0 
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0 0 



(4-SO) 





0 


_A_1 


(4-51) 

(4-52) 

(4-53) 


The state space averaged model over the entire period T is 


i = [dAj + d'AgJs + (dBj + (4-54) 

^ = (dCj + d'Cals 

where d = duty cycle ratio = Tqjj/(Ton+Tqj.p) 


d' a 1 - d 


T 


= T. 


ON 


OFF 


(4-55) 


The state space averaged model is perturbed and linearized 


in exactly the same way as for the buck-boost converter. The 
resulting small signal linearized model is 


£ =■ ' A 2 lX+ [DB~+ ~D^2^I ^ 

X = [DA^ + D'A2]x + [DB^ + 

+ [(A^ - A2)X + - B2)Vj]d (4=56) 

Vq - [DCj^ + D'C2]X 

Vq = [C^ - C2]Xd + [DC^ + D'C2 ]x 
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C 


Defining 

A = DAj + D'Ag 
B = DBj + D'Ba 
C = DCj + D'Ca 
results in 

X » + Bv^ 

+ [(A^ - A^)X + (B^ - B2)Vj]d 

^0 " - C2]M + Cx 

Using Laplace transforms gives 

x(s) - [SI - A]"S Vj(s) 

+ [SI - A]"^[(A^ - A2>X + (Bj_ - B2)Vj]d(s) 

Vq(s) - [C^ - C2]X aCs) + C x(s) 


(4-57) 


(4-58) 


(4-59) 


The state space model is derived from the aUsove two equa- 
tions and is shown in Figure 21 . 

The procedure for deriving the small signal equivalent 
circuit is exactly similar to the one used for the buck- 
boost converter. The four equations that result are 

^i " ^\i ^ci^^i ■ °\i ^ 

^^1 - ^C1 V 

DVci = (DR^^ " °^C1^1 "'O 

dll 

^ “dF " 


(4-61) 
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dv" v" 

r» k - 0 - - 

‘^IT'^L-R^ (*-63) 

The small signal equivalent circuit is described by the four 
equations above, as in the buck-boost converter, and is 
shown in Figure 22 . 


4.3 BOOST CONVERTER SMALL SIGNAL MODEL DERIVATION 
The procedure used in deriving the small signal model for 
the boost converter is exactly similar to that used for the 
buck-boost converter. The boost converter is shown in Figure 
23 . 

During the switch S is on, and the resulting cir- 

cuit is shown in Figure 24 (a). 

The equations describing the circuit are 




'Cl 


^C1 . "'l 


dt 

fici 

_ dt__ 

^ R, 

dt 
dv„ 


LI 

ii ht 




fhl- 


(Rci + R,) 




dt C(R^ + R^) 

V 


(4-64) 

X4^5)_ 

(4-66) 

(4-67) 


(4-68) 




Figure 21; Buck converter power stage small signal state 
space model. 
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During Tqj,j, the switch S is off and the resulting circuit is 
shown in Figure 24 (b). 

The equations describing the circuit are 


dt ’ LI T.1 LI T. " L 1 LI 

^ ^ ^L 
dt “ Cl " Cl 

dt L T.1 L ^^Cl ^ ^ 

Vc , ^C1 
■ l(Rc + V ^ 

fic \h \ 

dt “ C(R^, + \) “ C(R^ + R^) 

» ■ ^■■— f +!i!=_ 

° \ ^ ‘‘c + "4. 

The following vectors are defined 


^1 

"^Cl 


U =» [Vj] 
1 ” 


(4-69) 

(4-70) 

(4-71) 


(4-72) 

(4-73) 


(4-74) 


resulting in the following state space equations. 


ON 


X » A^x + 


T 

OFF 


X = A^x + 


(4-75) 


1 = Cj_x 


Z = C 2 X 


Figure 23 : Boost converter power stage 
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where 






LI 

!ci 

L 

Cl 


Cl 

LI 


.(Rci+V 


Cl 


LI 

1 

L 

0 

0 


0 

-1 


C(Rc + R^) 


(4-76) 


A. 


-(Rj^l+Rci) R, 


LI 

!ci 

L 

X 

Cl 


Cl 

LI 


-R, 


Cl 




-1 

LI 


7T 0 


1 

L 




L(Rc+Rj^) 

0 

CCR^Tr^ 


(4-77) 




^1 “ ^C1 r^ +k^ 




nT 


h - h 


\ 




(4-78) 

(4-79) 

i 

(4-80) 
( 4-81) 
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( 4 - 82 ) 


The state space averaged model over the entire period T is 
X » [dA^+d'A2]x + [dB^+d’B2]u 

2 “ tdC^+d'C2]x 


( 4 - 83 ) 


where 


d 


d' 


T 


duty cycle ratio = Tqjj/(Tqjj+Tqfp) 

1 - d 

T + T 
ON ^OFF 


( 4 - 84 ) 


The state space averaged model is perturbed and linearized 
in exactly the same way as for the buck-boost converter. The 
resulting small signal linearized model is 


_0 - A X + B Vj 

I ’ ^ 

X ■ A 3c + B Vj 

+ C(A^-A 2 ) X+ (Bj^-B 2 )Vj]d 


^0 " 

Vq = -C 2 ] X d + C X 

where A = DA^ + D'A 2 

B - DB^ + D'B2 
C = DC^ + D’C2 

Using Laplace transforms gives 


( 4 - 85 ) 


( 4 - 86 ) 

( 4 - 87 ) 

( 4 - 88 ) 
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x(s) - [SI -A]“^ B Vj(s) 

+ [SI-A]“^[(A^-A2)X + (B^-B2)V^]d(s) (4-89) 

Vq(s) - d(s) + d^(s) 

The state space model is derived from the above two equa- 
tions and is shown in Figure 25 . 

The procedure for deriving the small signal equivalent 
circuit is exactly similar to the one used for the buck- 
boost converter. The four equations that result are — 


LI- 


d£ * " ^C1 


^ dt " \l^l" 


DD 


, ^ch 

R.,+1 




, RR-Ld(D-D') 
+ - D'v^ -(- ^ 


V„d 


Cl “O' + \ ’ '0 


dv 


Cl 


‘1 d£ 


^1 " ^ 


*^^c •* ""o 

= — ■ “■‘l - r;; - V 


(4-90) 

(4-91) 


(4-92) 

(4-93) 


The small signal equivalent circuit is described by the four 
equations, above as in the buck-boost converter, and is shown 
in Figure 26 . 
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The power stage small signal models developed include 
the input filter state variables, whereas earlier models 
[3, 4, 5, 6] had treated the input filter only in terms of its 
output impedance and transfer fxinction. The models developed 
in this chapter are used to analyze and design a feedforward 
loop that includes the input filter state variables. 



66 








Chapter V 


IMPLEMENTATION OF THE FEEDFORWARD LOOP FOR A 

BUCK REGULATOR 

This chapter first presents an analysis that leads to a de- 
sign of the feedforward loops for a buck regulator. A small 
signal model that includes a general form of the feedforward 
loops is developed first. Analysis of the small signal model 
leads to a design of the feedforward loop. Implementation of 
the design is next discussed and two feedforward circuits 
are presented. The buck regulator alone is treated in this 
chapter, however the analysis and design procedure would be 
similar for the boost and the buck-boost regulators. 

5.1 STATE SPACE MODEL OF BUCK REGULATOR 

The state space model of the buck regulator is shown in Fig- 
ure 27 . In the figure the feedforward loop has as its in- 
put the input filter state variables ij^j and v^^ (the input 
filter inductor current and capacitor voltage respectively) . 
These two inputs are multiplied by the transfer functions 
c^Cs) and 03 ( 3 ) whose properties are yet to be determined. 
The feedback loop has as its inputs the output voltage and 
the output filter inductor current. The feedback control in 
this case is the two loop control (the standardized control 
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module or SCM) developed earlier [2,7,8] and discussed in 
Chapter 3. The error processor in Figure 27 is thus com- 
posed of the blocks labelled c^, Cg and the feedback, and 
has as its input information regarding the output filter and 
input filter state variables. The pulse modulator is repre- 
sented by its transfer function Fj^ [2,8]. The rest of Figure 
27 is the state space model of the buck power stage devel- 
oped in Chapter 4. 

The feedforward and feedback signals are added and fed 
to the pulse modulator. In physical terms this means sens- 
ing the small signal variations in input filter inductor 
current and capacitor voltage, processing these variations 
(as represented by the blocks C 2 (s) and Cg(s) in Figure 27) 
and adding the processed variations to the feedback signal. 
The total state feedforward/feedback error signal is then 

used to modulate the duty cycle of the switch for loop gain 

/ 

correction. 

The trcuisfer function Vq{s)/v^(s), Figure 27, is used 
to design the feedforward because it expresses clearly what 
the feedforward does; also the resulting design is indepen- 
dent of the feedback loop parameters. The generalized small 
signal model for the buck regulator is developed next and 
used to write the transfer function v^(s)/v^(s). 



< > 
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Figure 27 : State space model of the buck regulator 
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5.2 GENERALIZED SMALL SIGNAL MODEL OF BUCK REGULATOR 
The generalized small signal model of the buck regulator is 
shown in Figure 28 . The regulator is modelled according to 
the three basic functional blocks: power stage, error pro- 
cessor and duty cycle pulse modulator. The power stage model 
consists of two inputs: disturbances from the line and 

A. 

the duty cycle control d, and four outputs: the output vol- 

A ^ 

tage Vq , the output filter inductor current ij^, the input 
filter capacitor voltage and the input filter inductor 

A 

current ij^^ . The error processor has as its input informa- 
tion regarding the output filter and the input filter state 
variables. The transfer functions Fg, F^q and Fjj^^ constitute 
the two loop standardized control module (SCM) developed 
earlier [2,7,8], whereas the feedforward loop gains Cg and 
C 3 are as yet unknown. The error processor processes infor- 
mation regarding the state variables of the input and output 
filters and feeds a total error signal to the pulse modula- 
tor, whose transfer function Fj^ was developed earlier 
[2,7,8] . 

The power stage transfer functions F^^ etc. are written 
using the following equations: 


■'ll''! 

+ 

Ti^d = 

V 

0 

’’ 21 ''! 

+ 

= 



+ 

’’32^ - 

^C1 

’’41''! 

+ 

T42d = 

^L1 


(5-1) 
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Generalized Small Signal Model 


Figure 28: 
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Thus it can be seen that 


and 



0 


(5-2) 



H F 


11 



(5-3) 


5.2.1 Development of Power Stage Transfer Functions 
As can be seen from equations (5-1) , Tg^ , Tgj and T 41 

A A 

can be evaluated with d = 0 and the other four. with Vj = 0. 
The starting point for the evaluation of the transfer func- 
tions is the small signal equivalent circuit model for the 
buck regulator power stage developed in Chapter 4, Figure 22 


5.2. 1.1 Evaluation of T^j , "^31 

A 

These transfer functions are evaluated with d = 0 in Figure 
22 . The resulting circuit is shown in Figure 29 . In Fig- 
ure 29 the input filter has been replaced by its forward 
transfer function H(s) and its output impedance Z(s). 

s^Ll Cl + sCKR^j^ + R^^) + 1 


H(s) 


( 5 - 4 ) 
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Z (s) 


s LICIR^^ + ^1 

S^LlCl + sCKRj^^ + R^^) + 1 


(5-5) 


From the equivalent circuit of Figure 29 the following can 
be derived; 


v^ DRl( 1 + scR^)H 


11 


ij^ D(1 + sCRj^) 


H 


■21 


31 


V, 


Vci a^H 


V, 


(5-6) 


- LI 


A - a^H 


■41 - (R + sLl)A 

Vj LI 


where 


ai - s LCRj^ + sCRj^(Rj^ ''' CR^^ ^ 

A = + D^Z (1 + sCRj^) 

V 

D = duty cycle = = supply voltage 


(5-7) 
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Using equations (5-3) and (5-6) the following are derived: 


Fii = DR^d + sCR^) 
F2JL = D(1 + sCRj^) 

^31 " ^1 

A/H - a^ 

^41 “ Rj^j^ + SLI 


(5-8) 


In the derivation of equations (5-6) the resistance R^^ has 
been assumed negligibly small. This is not an unrealistic 
assumption since the ESR of the input filter capacitor (Rqj) 
can be assumed negligibly small compared with the other re- 
sistances; also in the derivation of the following is 
used: 


v. 


- V 


Cl 


^1 


+ 


sLl 


= 1 


LI 


(5-9) 



0 
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5.2. 1.2 . Evaluation of T 32 and T 4 ^ • 

These transfer functions are evaluated with Vj = 0 in Figure 
22 . The resulting circuit is shown in Figure 30 . 

From the equivalent circuit of Figure 30 the following 
are derived: 






V 

0 _ 

Vq(Rj^ - D^Z) (1 + 

SCR^) 

d 

DA 


k ^ 

Vq(I^ - D^Z) (1 + 

sCRj^) 

/V 

d 



^C 1 

_ + '^<1 - 

• sCRj^) ] 

A 

d 

<3 

1 


^L 1 

. 2V^(ai + + 

sCRj^) ] 

A 

d 

(Rj^j^ + SLI) 



(5-10) 


Using equations (5-1) and (5-10) the following are der- 
ived: 


12 


22 


V (Rt - D Z) (1 + SCR^) 


Vo(R|^ - D Z) (1 + sCRj^) 

drI 


32 


-ZV [a, + R_ (1 + sCR, ) ] 

O X Ij Li 


(5-11) 




Figure 30 
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ZV^[ai + + SCR^)] 

^LDR^ 


where ^ and a^ are as defined in equation (5-7) and 



R 


Cl 


= 0 


(5-12) 


5.2.2 Feedback Transfer Functions 

Transfer functions F^ , ®a(j/ c:oi^stitute the feed- 

back. A two loop standardized control module (SCM) cont- 
rolled [2,7,8] buck regulator was used to obtain experimen- 
tal results that are discussed later in this dissertation. 
The buck regulator used is shown in Figure 31 , and with re- 
ference to that figure the following are defined [2,7,8] : 

= snL 


AC sC|R^ 

7 = _i_( 2 

DC sCi^R^4 + R^ 

\ = ^ii//^i2 ' ^ = 


(5-13) 


^C ^13 sC. 
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F 


M 


fVl 

nM 


(Pulse Modulator Transfer Function) 


M = Constant depending on the type of control used. 


5.3 DESIGN OF THE FEEDFORWARD LOOP . 

A A 

The transfer function v^/v^ is used to design the feedfor- 

A 

ward. With = 0 the following equations are derived from 
Figure 28 ; 


[V^ =2Vl ^ ^ 


(5-14) 


and 


dF 


12 


o 

A, 

dF 


V 


32 


Cl 


■LI 


A 

dF 


(5-15) 


42 


/^ A A 

Substituting for v„ , i and d from (5-15) in (5-14) re- 

Cl LI 

suits in 


V 


^12 
(^^) F 
^ A ^ M 


V 1 - c.,/ 32,F,. - c.,/42,F„ 


(5-16) 
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In the absence of feedforward i.e. with Cg = Cg = 0 it can 
be seen that: 


^ ^ Vq(Ri, - D^ZHl ^ 

V + D^Z(1 + sCR^) 

X 1 L 


(5-17) 


The effect of peaking of the output impedance of the input 
filter Z is to cause a reduction in the term (Rj^ -D Z) and 
also an increase in the denominator, thus resulting in a 
substantial loss of loop gain. 

With feedforward the detrimental effect of peaking of Z 
could be avoided by a proper choice of the feedforward loop 
gains Cg and Cg . Choosing 


c., = 


-D 


V 

o M 


= 0 


(5-18) 


leads to 


V 


V 


Vm^^L - D Z) (1 t SCR^) 
DA 


1 + 


D‘ 


fZV^Cai+R^d + SCR^)] 


o M '' 




(5-19) 
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which can be simplified to 

^ ^ ^ SCR^)R^A 

DAai(I^ - D^Z) 


(5-20) 


Cancellation of the two terms leads to 


V V F„(l + SCR )R, 
O O M c L 


V 


Da, 


(5-21) 


Thus a proper choice of the feedforward loop gains has re- 
sulted in the transfer function Vq/v^ being completely inde- 
pendent of the input filter output impedance 2. It is also 
noted from equations (5-11) and (5-16) that at frequencies 
other than the resonant frequencies at which Z peaks, the 
gain of Fgg is fairly small since Z would be small at those 
frequencies. Thus the addition of feedforward would not af- 
fect, in any noticable manner, the open loop gain and phase 
margin at any frequency other than those at which Z peaks. 

The following points regarding the feedforward loop de- 
sign can be made: 

1. It has been shown analytically that a proper choice 
of feedforward loop gains results in eliminating com- 
pletely the effect of peaking of Z on the loop gain. 
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2. The gain C 3 = 0, thus the inductor current informa- 
tion is not needed, only the input filter capacitor 
voltage information is used. 

3 . The feedforward loop gains are independent of the in- 
put filter parameter values, and are free of any fre- 
quency dependent term. 

4. The feedforward loop gains are independent of the 
type of feedback control used. The pulse modulator 
transfer function Fjyj is, however, an integral part of 
the design and thus the compensation depends on the 
type of duty cycle control used. The feedforward loop 
design process is independent of the particular type 
of control used and thus the same design can be used 
for other types of control, for example for single 
loop control, current injected control and others. 

5.4 IMPLEMENTATION OF FEEDFORWARD 

The buck regulator used to obtain experimental results that 
are discussed later is shown in Figure 31 . The feedforward 
circuit processes the small signal variation across the in- 
put filter capacitor and adds this processed information to 
the feedback signal. 

Two circuit implementations of the feedforward design 
were used in making measurements and they are discussed 


next . 
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Figure 31: Buck regulator with feedforward used to obtain 

exprimental results 
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5.4.1 Nonadaptive Feedforward Circuit 

It was shown earlier that the feedforward loop gain is 

CgCs) = -D^/VqFj^ , equation (5-18) 

The nonadaptive feedforward circuit was developed for the 
buck regulator of Figure 31 . The key parameters of the re- 
gulator are as follows - 

Input-Output Parameters 

Vj = 25-40 volts Vq = 20 volts Pq = 40 watts 
Power Stage Parameters 

L = 230 micro H C = 300 micro F = 0.2 ohm 

= 0,067 ohm (nominal) Rj^ = 20 ohm (load) 

Pulse Modulator Parameters 

M = VjTqj^= 0.88 * 10^ V-sec 

Control Circuit Parameters 


E„ = 6.7 volts 

.K 

Rii 

’= 33.3 Rohm 

Ri 2 = 16.7 

Rohm 

■ = 2 Rohm 

Ri4 

= 47 Rohm 

R4 = 40.7 

Rohm 

n = 0.65 

1 

= 5600 picoF 

Cg = 0.01 

microF 


The buck regulator was operated in a predetermined duty cy- 
cle control mode (constant Vj T^^^ control ) , [2,7,8]. Substi- 
tuting for F , [2,8] and for D leads to 

M 
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-V nM 

c-(s) = — (5-22) 

where M = Vj Tqjj is constant. 

For the nonadaptive design the input voltage was kept cons- 
tant at Vj = 30 volts. Substituting in equation (5-22) it 
is calculated that C 2 ^(s) = -0.03. The nonadaptive feedfor- 
ward circuit implementation is shown in Figure 32 . The in- 
put to the circuit is the input filter capacitor voltage and 
a series capacitor (27 microF) blocks out the dc component. 
•The input is then multiplied by the gain of 0.03 implemented 
by the 5.1 Kohm and 164 ohm resistances. The feedforward 
signal available at the potential devider network is then 
subtracted from the feedback signal available at the output 
of the integrator in the feedback loop. The result is then 
fed to the pulse modulator. The capacitor voltage fed into 
the operational amplifier subtracting circuit consists of 
two components - a small signal variation and a component 
corresponding to the switching frequency. The feedback sig- 
nal is also at the switching frequency, but the amplitude of 
the feedback signal is large compared to the switching fre- 
quency information in the capacitor voltage, and thus the 
second component has negligible effect. It is to be noted 
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that the circuit of Figure 32 constitutes the feedforward 
circuit and also the summing junction shown in Figure 31 . 

The nonadaptive circuit has the advantage of being ex- 
tremely simple and easy to implement. The gain of the poten- 
tial devider in the feedforward circuit is, however, a func- 
tion of supply voltage and thus the circuit of Figure 32 
cannot be used at any other value of supply voltage. Mea- 
surements made using this circuit are discussed in the next 
chapter. 

5.4.2 Adaptive Feedforward Circuit 

The adaptive feedforward circuit is shown in Figure 33 . 
From equation (5-22) it is clear that changes in supply vol- 
tage Vj will change the gain of the feedforward loop, the 
circuit of Figure 33 implements the feedforward of equation 
(5-22) and adjusts the gain automatically as Vj changes. 

The input voltage in Figure 33 is allowed to vary bet- 
ween 25v and 40v. It is fed to a voltage devider and then 
squared. The input filter ' capacitor voltage consists of a 
large dc component and this is blocked out by the 27 microF 
capacitor in series with the feedforward path. The small 
signal variation and the small magnitude component at 
switching frequency are then devided by the squared input 
voltage. A pair of resistances provides the final gain; the 
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Figure 32; Nonadaptive Feedforward Circuit 




Figure 33: Adaptive feedforward implementation for a 

variable voltage supply. 
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feedfoirward signal now available is then added to the feed- 
back signal at the output of the integrator. As for the no- 
nadaptive circuit the switching frequency component in the 
input filter capacitor voltage is small compared to the cor- 
responding component in the feedback signal and thus its ef- 
fect on the duty cycle implementation is negligible. 

The gain of the feedforward circuit is thus a function 
of input voltage Vj and is made adaptive to changes in Vj . 
The feedforward is thus capable of tracking any variations 
in supply voltage. 



Chapter VI 


ANALYTICAL AND EXPERIMENTAL VERIFICATION OF THE 
FEEDFORWARD DESIGN 


This chapter presents extensive measurement data that 
was made to verify experimentally the feedfoirward design 
outlined earlier. The feedforward control was designed for a 
buck regulator and the same regulator is used in obtaining 
measurements . 

Measurement data pertaining to the open loop gain and 
phase margin are presented first; data made using the adap- 
tive feedforward circuit designed earlier confirm the adap- 
tive nature of the circuit. The audiosusceptibility and out- 
put impedance measurements presented next show that the 
feedforward significantly improves performance in both cate- 
gories. Lastly, measurements of transient response are in- 
cluded that show that the feedforward improves the transient 
response. 
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6.1 MEASUREMENTS OF Tffi NQNADAPTIVE FEEDFORWARD CONTROL 
The buck regulator with feedforward used to obtain experi- 
mental results is shown in Figure 34 , and is the same as 
that presented earlier in Chapter V. The parameters of the 
regulator are the same as given earlier in Chapter V with 
the single-stage input filter parameters specified as: 

Rj^j^ = 0.2 ohm LI = 116 micro-H Cl = 20 micro-F 
The feedforward circuit used was the nonadaptive feedforward 
circuit for a fixed input voltage discussed in Chapter V, 
with Vj = 30 V. The small signal open loop transfer function 
of the multiloop controlled buck regulator of Figure 34 
without feedforward can be expressed as [2,7,8] 

In equation (6-1) Fj^ 2 , ^2,2 power stage transfer 

functions, Fjj^, Fg and F^^, are the feedback control loop 

transfer functions, as discussed in Chapter 5 (section 5.2). 
The peaking of the output impedance of the input filter, 
2(s), affects the transfer functions Fj 2 , snd F^ 2 , seen 

below: 

Vo(R- - D^Z) (1 + sCR(,) 

p = 

D[ D*’Z( 1 + sCR^) + aj ] 

Vo(Rl - D^2)( 1 + sCRj^) 

p _ 

DRj^[ D Z( 1 + sCRj^) + a^ ] 


(6-3) 




Figure 34: Buck regulator with feedforward used to obtain 

experimental results 
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The peaking of Z(s) reduces the gain of both and and 
thus reduces the open loop gain at the resonant frequency of 
the input filter. 

The addition of the feedforward loop modifies the open 
loop transfer function as shown below -- 


G't(s) = 




D a, 


D a. 


(6-4) 


The addition of feedforward modifies the transfer functions 
^12. ^22 it can be seen ftom equation (6-4) that 
the open loop gain with feedforward is not affected by the 
peaking of the input filter output impedance Z(s). G'>p(s) is 
now independent of Z(s) and is a function only of the feed- 
back loop parcuneters and the power stage parameters, unlike 
Gj(s) of equation (6-1) which is affected by the peaking of 
Z(s). Equation (6-4) is also the open loop gain of the buck 
regulator without input filter, as can be seen from equa- 
tions (6-1), (6-2) and (6-3) by setting Z(s) = 0, and it can 
thus be concluded that the addition of the feedforward loop 
should eliminate completely the peaking effect of the input 
filter output impedance and that the open loop gain with 
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feedforward should be identical with the open loop gain 
without the input filter [10,11]. 

SINGLE STAGE INPUT FILTER 

Measurements of the open loop gain and phase margin of 
the buck regulator of Figure 34 were made with and without 
feedforward, and are presented in Figure 35 (a) and (b). The 
input filter resonates at around 3 KHz and results in dis- 
turbances in the open loop gain and phase margin at that 
frequency. The feedforward eliminates all these undesirable 
disturbances as is evident in the figures, thus providing 
close agreement with theory. It can also be seen from Figure 
35 that the characteristics with feedforward are almost 
identical to the gain and phase margin plots of the buck re- 
gulator without input filter, thus providing close agreement 
with the analytical prediction made earlier. 

TWO- STAGE INPUT FILTER 

The experiment was further extended to the buck regula- 
tor with a two-stage input filter. The two-stage filter of 
Figure 36 was used with the following parameter values: 


l~ 

0.2 ohm 

^1 = 

325 micro 

H 

Cj^= 200 micro F 

z~ 

0 . 02 ohm 


116 micro 

H 

C 2 i= 20 micro F 

3“ 

0.075 ohm 

(ESR 

o 

O 





DEGREES 
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The feedforward circuit used was the same nonadaptive feed- 
forward circuit used to obtain the measurements of Figure 35 
with the feedforward input being the voltage at capacitor C^. 

The other parameters of the circuit are the same as used to 
obtain Figure 25 . Measurements of the open loop gain and 
phase margin were made with and without feedforward and the 
results are shown in Figure 37 (a) and (b). The open loop 
gain and phase are affected at the resonant frequencies of 
the two stages of the input filter because the output impe- 
dcuice of the input filter Z(s) peaks at both resonant fre- 
quencies. The use of the feedforward circuit eliminates the 
detrimental effects of the input filter, as is evident from 
Figure 37 . 

REMARKS 

The following points regarding the above measurements 
are noteworthy: 

(1) Measurements of the open loop gain and phase margin 
show that the input filter output impedance causes distur- 
bances in the gain and phase margin at the filter resonant 
frequencies and the addition of feedforward eliminates these 
disturbances, providing close agreement with theory. 

(2) The feedforward compensation circuit is independent 
of the input filter parameter values. 
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Figure 36: Two- stage input filter. 
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(3) The feedforward compensation scheme is independent of 
the input filter configuration. It was demonstrated above 
that the same feedforward compensation network is equally 
applicable to a single stage input filter and a two-stage 
input filter. These observations lead to a stronger conclu- 
sion that the feedforward can provide effective compensation 
,for an unknown source impedance. For example, a preregula- 
tor which often has an unknown, dynamic output impedance can 
interact with a DC-DC converter downstream and result in 
system instability. The feedforward compensation scheme out- 
lined can be used to isolate the switching converter from 
the source thus preventing interaction between the switching 
converter and equipment upstream. 


6.2 ADAPTIVE FEEDFORWARD MEASUREMENTS . 

The buck regulator with feedforward used to obtain experi- 
mental results is shown in Figure 34 and the parameters of 
the circuit are as specified in section 6.1. The adaptive 
feedforward circuit for variable input voltage presented in 
Chapter 5 was used in obtaining measurements .The two-stage 
input filter of Figure 36 was used with the same parameter 
values as in section 6.1 with the feedforward input being 
the voltage at capacitor C 2 • The value of R 2 was changed 
to- 

R^= 0.2 ohm (ESR of C 2 plus external damping resistance) 


Measurements were obtained at four values of supply voltage 
using the same adaptive feedforward circuit in all cases - 
this was done to confirm the adaptive nature of the feedfor- 
ward circuit. 

A computer prograim was written to calculate the gain 
and phase margin of the open loop transfer function with and 
without two-stage input filter, at various input voltages. 
Ec[uation (6-1) was used to calculate the gain and phase mar- 
gin; setting Z(s) =0 in the equation gives the gain and 
phase margin without input filter. The following expression 
for the two- stage input filter was used: 


Z(s) = 


Z| + Rg ^ ^^2 
1 s 0 « ( Zy R9) L®C 


(6-5) 




where 


Zj(s) = (Ri + sL^) // (R 3 + 1/sCi) 


( 6 - 6 ) 


Figures 38 - 41 show the computed values of open loop gain 
and phase margin with and without the two-stage input filter 
for input voltages Vj = 2Sv, 30v, 35v and 40v. It can be 
seen that the two- stage input filter resonates at two fre- 
quencies and at each frequency the output impedance Z(s) 
peaks, thus causing sharp fluctuations in the open loop gain 
and phase margin. Measurements of the open loop gain and 
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phase margin at each of the above values of supply voltage 
were obtained with and without feedforward and are also 
plotted on the figures. It can be seen clearly that the ad- 
dition of feedforward removes the sharp fluctuations in open 
loop gain and phase margin caused by the input filter, pro- 
viding close agreement with theoretical prediction made ear- 
lier. The analytical prediction that the open loop gain and 
phase margin with feedforward are identical to the charac- 
teristics without input filter is also confirmed, as exami- 
nation Figures 38 - 41 show. 

The two stage input filter was modified so that 
Rj = 0.075 ohm (ESR of C 2 } and measurements of the open loop 
gain and phase margin with and without feedforward were 
made. Figure 42 shows the calculated values of open loop 
gain and phase margin together with the measured values at 
Vj = 25v, With the external damping resistance set to zero 
the effect of the input filter is seen to be more pro- 
nounced. Measurements without feedforward shown plotted on 
the figure also show the pronounced effect of the input fil- 
ter. The. addition of feedforward effectively eliminates the 
sharp fluctuations in gain and phase margin caused by the 
undamped input filter. 
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calculated values without 
input filter 


calculated values with input 
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(A) measured values with input 
filter, without feedforward 
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Figure 38; 


(a) Open loop gain at Vj=25v: Calculated 
values and measured values (A) without 
feedforward 
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Figure 38: (b) Open loop phase margin at 

Vj=25v: Calculated values and measured values 
(A) without feedforward 
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(c) Open loop gain at Vj=25v; Calculated 
values and measured values (A) with feedforward 


Figure 38: 
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Figure 38; (d) Open loop phase margin at 

Vj=25v: Calculated values and measured values(A) 
with feedforward 
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Figure 39: (a) Open loop gain at Vj=30v: Calculated 

values and measured values (z\) without 
feedforward 
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(b) Open loop phase margin at 

Vj=30v; Calculated values and measured values 

(A) without feedforward 


Figure 39: 


\ calculated values withou 
aV input filter 


calculated values with input 
filter, without feedforward 


(A) measured values with input 
filter and with feedforward 
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Figure 39: (c) Open loop gain at Vj =30v: Calculated 

values and measured values (A) with feedforward 
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Figure 39: (d) Open loop phase margin at 

Vi=30v: Calculated values and measured values 
(a) with feedforward 


Ill 


8 



FREQ (HZ) 


(a) Open loop gain at Vj=35v; Calculated 
values and measured values (A) without 
feedforward 


Figure 40: 
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Figure 40; (b) Open loop phase margin at 

Vj=35v: Calculated values and measured values 
(A) without feedforward 
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Figure 40: 


(c) Open loop gain at Vj=35v: Calculated 
values and measured values (A) with feedforward 
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Figure 40: (d) Open loop phase margin at 

Vj=35v: Calculated values and measured values 
(A) with feedforward 
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Figure 41: (a) Open loop gain at Vi=40v: Calculated 

values and measured values {^) without 
feedforward 
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Figure 41: (b) Open loop phase margin at 

Vj=40v: Calculated values and measured values 
(a) without feedforward 


calculated values without 
■ A input filter 


calculated values with input 
filter, without feedforward 


(A) measured values with input 
filter and with feedforward 


I I I I 


. . .J 


t r f .f I 


FREQ (HZ) 


Figure 41: (c) Open loop gain at Vj=40v: Calculated 

values and measured values (a) with feedforward 
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Figure 41: (d) Open loop phase margin at 

Vj=40v: Calculated values and measured values 
(a) with feedforward 
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The following observations regarding the measurements 
are made: 

(1) The peaking of the output impedance of the first 
stage is more pronounced and has a greater effect on the re- 
gulator than the peaking of the second stage. • 

(2) Both analysis and measurement results indicate that 
the open loop gain is higher at lower values of duty cycle 
D. This is further manifested by examining eqation (6-4). It 
shows that a lower value of D results in a higher gain. 

(3) The effect of input filter peaking varies with the 
input voltage. This is explained by noting that both and 

equations (6-2) and (6-3) depend on the duty cycle D. 

The effect of peaking of 2(s) is to cause a reduction in the 
1 

term (Rj^ -D 2) and the amoxint of reduction would be greater 
if D is larger. Consequently it is expected that at small 
values of Vj , when D is larger, the effect of peaking would 
be more pronounced. This is confirmed by examination of Fig- 
ures 38-41 . 

(4) The addition of the ’ feedforward loop effectively eli- 
minates the perturbation in open loop gain and phase caused 
by the input filter. The feedforward works effectively at 
all four values of supply voltage Vj . Since the same feed- 
forward circuit was used in making all the measurements of 
Figures 38 - 42 the adaptive nature of the feedforward cir- 


cuit is confirmed. 
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Figure 42: (a) Open loop gain at Vj=25v: Calculated 

values and measured values (a) without 
feedforward 
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Figure 42: (b) Open loop phase margin at 

Vi=25v: Calculated values and measured values 
(a) without feedforward 
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Figure 42: (c) Open loop gain at Vj =25v: Calculated 

values and measured values (A) with feedforward 
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(d) Open loop phase margin at 

Vj=25v: Calculated values and measured values 

{^) with feedforward 


Figure 42: 
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(5) From the results presented in section 6.1 and from 
measurements presented in "tois section it is logical to con- 
clude that the adaptive feedforward circuit can provide ef- 
fective compensation for an arbitrary, unknown source impe- 
dance for a variable supply voltage. The adaptive 
feedforward - circuit has been shown to be able to track the 
supply voltage and adjust the gain in accordance with supply 
voltage changes. The adaptive feedforward circuit can effec- 
tively isolate the switching regulator from its source impe- 
dance, thus preventing any interaction between the switching 
converter and equipment upstream. 

6.3 MEASUREMENTS OF CLOSED LOOP INPUT -TO -OUTPUT TRANSFER 
FUNCTION ( AUDIOSUSCEPTIBILITY ) . 

The closed loop input- to -output transfer function (audiosus- 
ceptibility) of a switching regulator is an important char- 
acteristic. It refers to the regulator's ability in attenu- 
ating small signal sinusoidal disturbances propagating from 
the regulator input to its output. The gain of the closed 
loop input- to-output transfer function should be as small as 
possible; thus the regulator will effectively attenuate 
noise at the input so as not to affect operation of the re- 
gulator payloads. Unfortunately, as pointed out in Chapter 
2, the peaking of the output impedance of the input filter 
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and the peaking of the forward transfer function of the in- 
put filter increase the audio susceptibility. Measurements 
of audio susceptibility with and without feedforward were 
made using the two-stage input filter of Figure 36 with a 
damping resistor so that = 0.2 ohm. Measurements were 
made by injecting a small sinusoidal signal at the input to 
the converter and then using a HP network analyser to mea- 
sure the audiosusceptibility [8]. The feedforward circuit 
used was the adaptive feedforward circuit of Figure 32 and 
the buck regulator used is shown in Figure 34 , with its 

parameters as specified in section 6.1. 

Figures 43 - 46 show the measured values of audiosus- 
ceptibility with and without feedforward at four values of 
supply voltage Vj = 2Sv, 30v, 35v and 40v using the same 

feedforward circuit in all cases. The top trace in each fig- 
ure is the plot without feedforward and it can be seen 
clearly that the audiosusceptibility is degraded 
at the two resonant frequencies where the outp- 
ut impedance 2(s) and the transfer function H(s) of the in- 
put filter peak , with the first stage resonating around 
600 Hz and the second stage around 3 KHz. It is also evident 
from the figures that the audiosusceptibility is dependent 
on duty cycle D or the supply voltage. At higher values of 
supply voltage when the duty cycle is low the audio suscepti- 
bility is lower, specially at the lower frequencies. 
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Figure 43: Measurement of audiosusceptibility [ Vq ( s )/v^ ( s ) ] 

with and without feedforward at Vj =25v 
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Figure 45: Audiosusceptibility with and without feedforward 

at Vj=35v (measurements) 
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The addition of feedforward substantially improves the 
audiosusceptibility, specially at the lower frequencies, as 
is evident upon examination of Figures 43 - 46 . The peaking 
effect of audiosusceptibility with feedforward loop is how- 
ever, more pronounced at the two resonant frequencies of the 
two-stage input filter; this is explained by noting that in 
equation (2-6), Chapter 2, the audiosusceptibility is shown 
to be affected both by the peaking of Z(s) and by the peak- 
ing of the forward transfer function H(s) of the input fil- 
ter. The peaking of H(s) cannot be controlled in any fashion 
by the addition of a feedforward loop since the control loop 
is not affected by H(s); thus the feedforward loop is effec- 
tive in cancelling 2(s) while the peaking effect of H(s) is 
manifested in audiosusceptibility. Figures 47 and 48 show 
the transfer functions H(s) and Z(s) of the two-stage input 
filter used. The peaking of H(s) at the two resonant fre- 
quencies is clearly seen. 

The two-stage input filter was modified by removing the 
damping resistor so that = 0.075 ohm. Measurement data of 
the audiosusceptibility with and without feedforward are 
shown in Figures 49-52 . The top trace in all these fig- 
ures is the closed loop input- to- output transfer function 
without feedforward and it can be seen that the gain is 
higher than in the earlier case (with damping resistance) . 
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Figures 53 and 54 show the forward transfer function H(s) 
and the output impedance Z(s) of the input filter used. Com- 
paring Figures 47 and 48 with Figures 53 and 54 it can be 
seen clearly that both Z(s) and H(s) peak at significantly 
higher values when the external damping resistance is re- 
moved . 

The addition of feedforward sxibstantially improves the 
audiosusceptibility as is evident from Figures 49-52 . The 
audiosusceptibility with feedforward peaks at the two reso- 
nant frequencies of the two- stage input filter as before, 
but in this case the peaks are higher. This is explained by 
noting that H(s) of the two-stage input filter without damp- 
ing resistance peaks at a significantly higher value, as 
shown in Figures 53 and 54 , than that with a damping resis- 
tance as shown in Figures 47 and 48 ; thus the effect of 
H(s) on the closed loop gain would be expected to be greater 
in the former case. 

It can therefore be concluded that the addition of 
feedforward significantly improves the audiosusceptibility, 
specially at the lower frequencies. The audiosusceptibility 
with feedforward is affected by the peaking of H(s) since 
the effect of peaking cannot be eliminated via any control 
means. The same adaptive feedforward circuit was used in 
making all the closed loop gain measurements mentioned 
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Figure 47; H(s) of two stage input filter, with R3=0.2 ohm 
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Figure 48: Z of two stage input filter, with R 3 -O .2 ohm 
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Figure 49; 


Audiosusceptibility with and without feedforward 
at Vj=25v (measurements) 
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Figure 51: Audiosusceptibility with and without feedforward 

at Vj=35v (measurements) 
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Figure 53 : 


H(s) of two stage input filter, 83=0.075 ohm 
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above, this confirms anew the adaptive nature of the feed- 
forward circuit. 

6.4 MEASUREMENTS OF OUTPUT IMPEDANCE . 

The closed-loop output impedance of a switching regulator 
should be as low as possible in order that the regulator be- 
have as much as an ideal voltage source as possible. Howev- 
er, as pointed out in Chapter 2, the peaking of the output 
impedance of the input filter increases the closed- loop out- 
put impedance of the regulator. 

Measurements of the regulator output impedance with and 
without feedforward were made, using the two-stage input 
filter of Figure 36 with = 0.075 ohm (ESR of C^) . Mea- 
surements were made by injecting a small signal sinusoidal 
disturbcuice at the output in parallel with the load of the 
regulator, and then using a HP network analyser to measure 
the corresponding voltage and current [8,9]. The feedforward 
circuit used was the adaptive feedforward circuit of Figure 
33 ,. Figures 55 - 58 show the measured values of output 

impedance with and without feedforward at four values of 
supply voltage Vj = 25v, 30v, 35v and 40v using the same 

feedforward circuit in all cases. 

It can be seen clearly from the figures that the output 
impedance is increased at the two resonant frequencies of 
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Figure 56: Output impedance with and without feedforward 

(X) at Vj=30v 
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Figure 58: Output impedance with and without feedforward 

(X) at Vj =40v 



the two- stage input filter. This is a consequence of the 
disturbances in the loop gain produced at those frequencies 
by the peaking of the input filter output impedance Z(s). As 
was seen earlier in sections 6.1 and 6.2 the effect of Z(s) 
on the open loop gain depends on the duty cycle D, at higher 
values of D i.e. lower supply voltages, the effect of Z(s) 
on the open loop gain is higher. Thus the effect of Z(s) on 
the output impedance would be greater at lower supply vol- 
tages, and this is experimentally confirmed as can be seen 
from Figures 55 - 58 . 

The addition of feedforward almost totally eliminates 
the undesirable perturbations in the output impedance char- 
acteristic at all supply voltages. 

6.5 MEASUREMENT OF TRANSIENT RESPONSE AND STARTING OF THE 
REGULATOR . 

In this section measurements of output voltage and other 
parameters for a step change in input voltage or load are 
presented with and without feedforward. 

6.5.1 Small Amplitude Transient Response Measurements . 
Photographs of the output voltage ripple and other parame- 
ters are presented with and without feedforward for two cas- 
es: 

(1) Step change in supply voltage. 


(2) Step change in load. 

The step changes mentioned aibove are small enough so that 
the output voltage remains in regulation throughout the 
transient response period. 

6.5. 1.1 Step Change in Input Voltage. 

The buck regulator used is shown in Figure 34 with the par- 
ameters as specified in section 6.1. The input filter used 
was a single-stage input filter with the following parame- 
ters: 

^Ll ~ -- 325 micro H Cl = 220 micro F 

The adaptive feedforward circuit of Figure 33 was used in 
making the measurements. 

The input voltage was abruptly switched from Vj = 30v 
to 40v euid photographs of output voltage ac ripple, input 
filter capacitor voltage, output filter inductor current and 
control voltage without using a feedforward loop were taken 
as shown in Figures 59 (a) and (b) , 60 (a) and (b), respec- 
tively. The control voltage is the input to the pulse modu- 
lator; without feedforward it is the output at the integra- 
tor in the feedback loop while with feedforward it is the 
sxim of the above signal and the feedforward signal. Figures 
61 (a) and (b) , 62 (a) and (b) show the photographs of the 
same variables with feedforward control. Comparing, for ex- 
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ample, the photographs of the output voltage ripple with and 
without feedforward. Figures 59 (a) and 61 (a), it can be 
clearly seen that the transient response is improved with 
the addition of feedforward. The amount of overshoot is 
less with feedforward. The magnitude of the oscillations in 
the output voltage caused by the interaction between the in- 
put filter and the regulator control are also lessened with 
the addition of feedforward. 

Comparison of the photographs of input filter capacitor 
voltage, output filter inductor current and control voltage 
do not reveal much difference between the two cases (with 
and without feedforward). This may be explained by noting 
that the gain of the feedforward loop is fairly small ( 0.03 
for 7j * 30v ) euid thus the feedforward signal would be 
fairly small in amplitude. The addition of such a small am- 
plitude signal to the fairly large amplitude waveforms re- 
corded on the photographs will not show very clearly. The 
output voltage ripple, however, is small in magnitude and 
the effect of adding feedforward shows clearly. A computer 
program was written to simulate the step change in voltage; 
results from the program are presented in the next chapter 
and show close agreement with the measurements. 

Thus it is concluded that feedforward improves tran- 
sient response for the case where the supply voltage is sub- 
jected to a step change. 
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CbJ 

Figure 59: Output voltage ripple (a) and input filter 

capacitor voltage (b) without feedforward 

X — 1 msec/di V 


( a ) Y — 0 . Iv/di V 

(b) Y — 5v/div 


X — 1 msec/div 
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^b) 

Figure 60: Output filter inductor current (a) and control 

voltage (b) without feedforward 

(a) Y— O.SA/div — 1 msec/di v 

(b) Y — 0.5v/div 


X — 1 msec/div 



ISO 



(b) 

Figure 61: Output voltage ripple (a) and input filter 

capacitor voltage (b) with feedforward 

(a) Y — O.lv/div X — 1 msec/div 

(b) Y — 5v/div 


X — 1 msec/div 



Figure 62: Output filter inductor current (a) and control 

voltage (b) with feedforward 

(a) y — 0.5 A/div ’ X— 1 msec/div 

(b) Y — 0.5v/div X — 1 msec/div 



6.5. 1.2 Step Change in Load. 

The buck regulator used is shown in Figure 34 , with the 
parameters as specified in section 6.1 with the following 
changes: 

Vj = 25v Rj 4 =3.69 Kilo ohm 

A single stage input filter was used with the following par- 
ameters: 

=0.2 ohm Ll = 325 micro H Cl = 100 micro F 
The adaptive feedforward circuit of Figure 33 was used in 
making the measurements. 

The load was switched repetitively between Rj^ =10 ohms 
and Rj^ = 20 ohms using a transistor switch. Figures 63 (a) 

and (b) show the photographs of the output voltage ripple 
without and with feedforward, respectively, as the load is 
switched. The output voltage ripple without feedforward. 
Figure 63 , shows distinct oscillations caused by the inter- 
action between the input filter and the regulator control 
loop. The oscillation frequency coincides with the input 
filter resonant frequency. These oscillations are eliminated 
with the addition of feedforward, as is evident from Figure 



(a) 



Figure 63: Output voltage without feedforward (a) and with 

feedforward (b) 
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6.5.2 Large Amplitude Transient Response Measurements . 
Photographs of the output voltage ripple and other parame- 
ters are presented with and without feedforward for two cas- 
es: 

(1) Large step change in supply voltage. 

(2) Starting of the converter. 

The cases mentioned above are large signal changes so that 
the voltage regulation is momentarily lost for part of the 
transient response period- 

6.5.2.1 Large Step Change in Supply Voltage. 

The buck regulator used is shown in Figure 34 with the par- 
ameters as specified in section 6.1. The input filter used 
was the single-stage input filter of section 6.5. 1.1. The 
adaptive feedforward circuit of Figure 33 was used in making 
the measurements with feedforward. 

The input voltage was abruptly switched from Vj = 40v 
to 25v and photographs of the output voltage were made with 
and without feedforward; Figures 64 (a) and (b) show the 
photographs . 

The step change in supply voltage .is large enough to 
cause the regulator to lose regulation -- the output voltage 
drops by about 1.5v before the regulator recovers. This may 
be explained by noting that for such a large change in sup- 
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ply voltage the input filter capacitor voltage drops down 
close to 20v, and since this value is lower than the design 
range of 25v - 40v for the regulator supply, the regulator 
loses regulation. Further details are given in the next 
chapter which presents results from a computer program writ- 
ten to simulate this large step change. The results show 
close agreement with the measurements presented here. 

Figures 64 (a) and (b) show that the behavior of the 
regulator is similar with and without feedforward for such a 
large step change in supply voltage. This is expected since 
the regulator control loop momentarily loses its control 
fiinction during this transient period. Since the feedforward 
is designed to compensate for the effects of input filter 
interaction via a duty cycle modulation scheme it is expect- 
ed that the feedforward does not contribute anything under 
these conditions. Computer based simulation results present- 
ed in the next chapter confirm the measurement results. It 
is also to be noted that the feedforward does not have any 
detrimental effect on the transient response. 

6. 5.2.2 Measurements of the Start Up Behavior of the 
Regulator 

This section investigates the start-up behavior of the 
switching regulator. The regulator used is shown in Figure 



Cb) 

Figure 64: Output voltage without feedforward (a) and with 

feedfoirward (b) 


(a) 

Y — O.Sv/div 

X — 0.5 msec/di V 

(b) 

Y— 0.5v/div 

X — 0.5 msec/div 



34 with the parameters as specified in section 6.1. The sin- 
gle stage input filter of section 6. 5. 1.1 was used with the 
supply voltage set at 30v. The adaptive feedforward circuit 
of Figure 33 was used in making the measurements with feed- 
forward. 

Prior to starting, the output voltage and the output 
filter inductor current were both zero. Photographs of the 
output voltage and output filter inductor current were made 
with and without feedforward suid are shown in Figures 65 (a) 
euad (b) 66 (a) and (b) . The output voltage without feedfor- 
ward builds up from zero to 20v (regulated output) in about 
3 msec. The inductor current rises sharply at starting but 
is limited to about 6.0A by the peak current protection cir- 
cuit built in with the regulator. It settles down to its 
steady state value in about 4 msec. 

The photographs with feedforward show similar behavior 
- thus the feedforward does not contribute significantly to 
this transient period nor does it present any detrimental 


effects. 
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(a) 



C^) 

Figure 65: Output voltage (a) and output filter inductor 

current (b) without feedforward 

(a) — 5v/div X — 0.5 rasec/div 

(b) Y — 1 A/div X — 0.5 msec/di v 
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(a) 



(b) 

Figure 66: Output voltage (a) and output filter inductor 

current (b) with feedforward 

(a) Y—5v/div 

(b) Y— 1 A/div 


X — 0.5 msec/di V 
X — 0.5 msec/div 
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6.6 CONCLUSIONS 

Extensive measurements made to verify experimentally the 
feedforward design are presented in this chapter. Measure- 
ments of the open loop gain and phase margin, closed loop 
gain, output impedeuice and transient response confirm the 
effectiveness of the feedforward in improving performance, 
as predicted in the analysis. The following points are. note- 
worthy ; 

1. The feedforward eliminates the detrimental effect on 
open loop gain and phase margin of the output impe- 
dance of the input filter. 

2. The feedforward circuit is shown to be independent of 
input filter parameters and also independent of input 
filter configuration. 

3. The feedforward effectively eliminates the interac- 
tion between an unknown dynamic spurce impedance and 
the regulator control loop. 

4. The closed loop input- to-output transfer function 
(audiosusceptibility) and output impedance are both 
improved significantly by the addition of feedfor- 
ward. 

5. The addition of feedforward improves transient res- 
ponse in those cases where the interaction between 
the input filter and the control loop degrades the 



response. Examples of the above are small step chang- 
es in supply voltage amd load, and results for these 
cases are included. For large signal transient behav- 
ior the feedforward does not in any way degrade the 
performance — examples of these are a large step 
change in supply voltage and the start up behavior of 
the regulator. 

No detrimental effects have been observed due to the 
use of the proposed feedforward compensation scheme 
through the course of study and through extensive ex- 
periments when the system was subjected to different 

I 

forms of small and large signal disturbances. 



Chapter VII 


DIGITAL SIMULATION OF BUCK REGULATOR FOR 
TRANSIENT ANALYSIS 


This chapter deals with a program developed to simulate 
the real time behaviour of the regulator with and without 
feedforward control under both steady state and transient 
conditions. Results for a step change in supply voltage are 
included that are in close agreement with the measurements 
presented in the last chapter. 


7.1 DESCRIPTION OF ^ METHOD 

The starting point is the definition of the buck regu- 
lator in terms of state equations t 12, 13, 14]. The switching 
regulator used was the one used in Chapter 5 and 6 and is 
shown in Figure 34 The buck regulator without feedforward is 
described by a set of state equations. The six state varia- 
bles used are: 

x(l) = ij^j (current ininput filter inductor) 

x(2) * (input filter capacitor voltage) 

x(3) = Bq (control voltage, input to pulse modulator) 

x(4) = ij^ (output filter inductor current) 

x(5) * V (output filter capacitor voltage) 
c 

x(6) - e (compensation loop capacitor voltage) 

XV 
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The control voltage is the output of the integrator in the 
feedback control loop if no feedforward is used and it is 
the stun of the integrator output and the feedforward signal 
if feedforward is used. The state variable- e is the voltage 
at the capacitor in the compensation loop. The pulse mo- 
dulator has as its input the control voltage e^ and it con- 
verts the voltage to a duty cycle signal d(t). The transfer 
fimction of the pulse modulator is a constant if the supply 
voltage is a constant, [8], for the constant volt-sec. 
(VjTqj^) control mode that was used. Thus it can be seen that 
the six state variables defined above completely describe 
the buck regulator system of Figure 34 . 

The state equations for the complete system without 
feedforward are developed [12,13,14] for the two time per- 
iods Tqj, and Topp. 

During Tq^j the transistor switch SI is on while diode 
S2 is off and the system is characterized by the following 
state equations: 

X = FI X + G1 u (7-1) 


where 
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(7-6) 


In equation (7-6) is the supply voltage, Ejj is the refer- 
ence voltage, Eq is the saturation voltage drop of the power 
transistor emd E^ is the conduction voltage drop across the 
diode . 

During Tqpp the transistor SI is off while diode S2 is 
conducting auid the equations are: 

X = E2 X + G2 u (7-7) 

where 

F2 =s Fi except for 

F2 (2,4) = 0 

F2 (3,2) = 0 (7-8) 

F2 (4,2) - 0 


and 


G2 = 0 except for 
G2(l,l) = l/Ll 


G2(3,2) 




G2(3,4) 


n 




(7-9) 
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G2(4,4) = - i 

Since u is not a function of time the solution to the state 


equations (7-1) 

and 

(7-8) can be written as: 


During 

x(t 

+ T) = * x(t) + Dp * u 

(7-10) 


X (t 

+ = <fN * 2 i(t) + * u 

(7-11) 


where are the state transition matrices defined for 

a small fixed step size T [12] and 

(7-12) 


Equations (7-11) and (7-12) are used to simulate the behav- 
ior of the regulator without feedforward. Starting at the 
beginning of the ON time period equation (7-12) is used to 
propagate the state. The step size T is defined as half the 
ON and OFF times and this determies the matrices 
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and Dp. Substitution into equation (7-11) along with 
the initial values at the beginning of the ON period propa- 
gates the state through time T. Similar substitution propa- 
gates the state through the ON period [12,13,14]. For the 
type of duty cycle contol used, constant VjTq^j control, the 
ON time is fixed whenever Vj is fixed. 

At the end of the ON period, equation (7-10) is used to 

propagate the state during the OFF period, with the initial 

condition being the state at the end of the ON interval. The 

end of the OFF interval is determined as the point when the 

/ 

control voltage e^ equals the comparator yoltage which 
was 7.0v in this case. Newton's iteration is performed to 
find the exact point in time when e^ equals [12,13,14]. 

The equations describing the system with feedforward 
are very similar to those without feedforward. The nonadap- 
tive feedforward circuit of Chapter 5 was used in the simu- 
lation program and it is shown in Figure 67 . The output of 
the capacitor is included as a state variable. The state 
variables are thus defined as: 

x(l) = i^i 
x(2) = Vj,j 

x(3) = (output of ) 
x(4)= Sj, 
x(5)= 


(7-13) 
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x(6)= 

h(7)- 

The equations during Tqj^ and TQpj,are the same as equations 
(7-1) euid (7-8) with the following definitions: 




i 
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K3 


'f2 


Rfl + 


"f2 


(7-16) 


In equation (7-15) kl, k2 and u are as defined earlier. 

F2= FI except for 
F2(2,5)= 0 
F2(3,S)= 0 
F2(4,2)= 0 

F2(4,S)= -Rj,*kl/ Cj (7-17) 

F2(5,2)= 0 

and 

G2= 0 except for 
G2(l,l)= 1/Ll 

G2(4,2)= l/C|.Rj^ (7-18) 

G2(4,4)= n/c;.R^ 

G2(5,4)= -1/L 

The solutions of the state equations and the procedure 
for propagation of the state from cycle to ccle is identical 
to the procedure for the case without feedforward. 

The discrete simulation technique outlined [12,13,14] 
lends itself easily to simulation of transient response for 
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Figure 67: Nonadaptive feedforward circuit 
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a step change in supply voltage. The change in supply vol- 
tage is made to occurr at the start of a new ON period; thus 
all that is necessary to do is to change the value of Tqj^ 
before the start of a new cycle. 

7.2 DESCRIPTION OF PROGRAM 

A flowchart of the program [12,13,14] is shown in Figure 68 
while a listing of the programs is included in the Appendix. 
The flowchart is for the simulation without feedforward, the 
simulation with feedforward proceeds in an identical fash- 
ion. 

A minimum off time duty cycle control is implemented in 
the program. The regulator also has peak current protection 
and this capability is also programmed. During the ON period 
the current in the switch may rise above the set value of 
the peak current. Newton's iteration [12,13,14] is performed 
to find the exact point in time where this occurr s and the 
program will automatically terminate the ON time calculation 
when the current exceeds the set value . 
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Figiire € 8 : Flow chart for discrete 

simulation 
























Vj = VISWIT (NSWIT) 

'"on “ ^iW^I 

u(l) = Vj 

Calculate <^p»Dp 
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^OFF 


switching instant 
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State vector 
propagation 
during ON time. 



x(t+T) = 

+ * a(t) 

t = t-KP 


Figure 68: continued 
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7.3 SIMULATION RESULTS 

The simulation program was used to simulate the behavior of 
the system with and without feedforward for a small step 
change in supply voltage from 30v to 40v and also for a 
large step change from 40v to 25v. Simulation results pre- 
sented are in close agreement with the experimental results 
presented in Chapter 6. 

7.3.1 Simulation of steady state operation 
The computer program discussed Ceui be used to simulate the 
steady state operation of the regulator system [12,13,14]. 
An accurate picture of the behavior of the system is ob- 
tained since the simulation method used is an exact method. 
The buck regulator of Figure 34 with feedforward was 

I 

simulated with the following single-stage input filter par- 
ameters : 

Rlj^= 0.2 ohm LI = 325 microH Cl =220 microF 
The supply voltage Vj was set equal to 30v and the other 
parameters were as presented in section 6.1. Figures 69 (a) 
- (g) show the plots of the output voltage, input filter in- 
ductor current, input filter capacitor voltage, output fil- 
ter inductor current, the feedforward voltage Vj , the con- 
trol voltage and the voltage er. It is noticed that the 
switching period is around 50 microsecs and that the output 
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voltage ripple is arovind 60 milivolts. The input filter in- 
ductor current has a small magnitude ripple component at the 
switching frequency, as does the input filter capacitor vol- 
tage. The voltage fed forward v^ consists of the switching 
frequency component, however its magnitude is around 50 mi- 
livolts. The switching frequency component of the control 
voltage is seen to be around 800 milivolts and thus the ef- 
fect of adding the switching frequency component from the 
input filter capacitor voltage to the control voltage would 
be negligible as discussed in Chapter 5 ; since the magni- 
tude added would be around 1.5 milivolts (v^ multiplied by 
the feedforward gain) . 

7,3.2 Simulation for a small step change in supply voltage 
The same buck regulator was used with its parameters as spe- 
cified in section 7.3.1. A value of 0.6 ohm for was 

used. The supply voltage Vj was abruptly switched from Vj = 
30V to 40V, Experimental results for this step change are 
given in Chapter 6 where a value of 0.2 ohm was specified 
for Rjjj- A value of 0.6 ohm was used in the simulation as 
this represents the combined winding resistance and source 
impedance. The values of the other parameters used are: 

Toppi( minimum) = 5 microseconds 

I-(maximum) - 6 amps, (set value of peak current) 

Eq = 0.2 volt Eq = 0.7 volt (7-19) 








VFCV 
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Figures 70 (a) - (d) are plots of the computed values 
of output voltage, input filter capacitor voltage, output 
filter inductor current and the control voltage, all without 
feedforward. 

The nonadaptive feedforward design of Figure 67 was 
used in the simulation of transient response with feedfor- 
ward. The feedforward circuit parameters used at Vj = 30V 
were; 

Cj = 27 microF =5.1 Kohm = 164 ohm 

when the supply voltage is switched to 40V the value of R^^ 
in the program is changed to 90.97 ohm at the switching ins- 
tant thus changing the feedforward circuit gain in accor- 
dance with the feedforward design. 

Figures 71 (a) -(d) are plots of the computed values of 
the same varieJsles with feedforward. Experimental measure- 
ments for the same step change in supply voltage were pre- 
sented in Chapter 6 and are repeated here for convenience; 
Figures 72 (a)-{d) are the measurement waveforms without 

feedforward loop while Figures 73 (a) -(d) are for the case 

with feedforward loop. 

Comparing the plots of the computed values with and 
without feedforward it is clearly seen that the analytical 
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Figure 70: (b) Input filter capacitor voltage simulation 

without feedforward 





inductor current simulation 
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Figure 71 


(c) Output filter inductor current simulation 
with feedforward 
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Figure 71: 


(d) Control voltage simulation with 
feedforward 
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(*) 



(b) 

Figure 72: Measurement without feedforward: output voltage 

(a) and input filter capacitor voltage (b) 

(a) Y— O.lv/div . X — 1 msec/div 


(b) Y — 5v/div 


X — 1 msec/div 




(d) 

Figure 72: Measureraen't without feedforward: output filter 

inductor current (c) and control voltage (d) 

(a) Y — 0.5 A/div X — 1 lasec/div 

(b) Y — 0,5v/div X — 1 msec/div 
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(») 



(b) 

Figure 73: Measurement with feedforward: output voltage (a) 

and input filter capacitor voltage (b) 

(a) Y — 0-lv/div X — 1 msec/div 

(b) Y — 5v/div X — 1 msec/div 
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results agree fairly closely with the experimental data. The 
analytical results show that the feedforward is effective in 
reducing the output voltage oscillations — the amplitude of 
the oscillation is reduced (Figures 70 (a) and 71 (a)); this 
is confirmed by the experimental results (Figures 72 (a) and 
73 (a)). The plots of the input filter capacitor voltage, 
output filter inductor current and the control voltage are 
also in excellent agreement with the experimental data. Also 
in excellent agreement with the measured data is the obser- 
vation that the plots of filter capacitor voltage, inductor 
current and control voltage do not show any noticable dif- 
ference between the two cases (with and without feedfor- 
ward) . An explanation for this observation was provided in 
Chapter 6 . 

It can thus be concluded that the addition of feedfor- 
ward improves the transient response for the step change in 
Vj , as demonstrated both from the analytical result as well 
as experimental data. 

7.3.3 Simulation for a Large Step Change in Supply Voltage 
The regulator used was the one used in Section 7.3.1 with 
the same parameters. The single stage input filter parame- 
ters used were: 

= 1 ohm LI = 325 microH Cl = 220 microF 
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The supply voltage was abruptly switched from Vj = 40V to 
25V. Experimental results for this step change are given in 
Chapter 6 where a value of 0.2 ohm was specified for A 
value of 1 ohm was used in the simulation as this represents 
the combined winding resistance and source impedance. The 
other parameters are as specified in Section 7.3.1. 

Figures 74 (a-d) are plots of the computed values of 
output voltage, input filter capacitor voltage, output fil- 
ter inductor current and the control voltage, all without 
feedforward. 

The nonadaptive feedforward design of Figure 67 was 
used also in the simulation of transient response. When the 
voltage is switched from 40V to 25V the value of Rj 2 is 
changed from 90.97 ohm to 230.27 ohm in order to change the 
gain of the feedforward circuit in accordance with the feed- 
forward design. Figures 75 (a-d) are plots of the computed 
values of the same variables _with feedforward. Experimental 
measurements for the same step change in supply voltage were 
presented in Chapter 6 and are repeated here for conveni- 
ence; Figure 76 shows the measured output voltage ripple 
without and with feedforward.. 

Comparing the plots of the output voltage ripple with 
the measured data it is clearly seen that the analytical re- 
sults agree closely with the experimental data. It can be 
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Figure 74: 


(c) Output filt« 
without feedforwai 
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Figure 74: (d) Control voltage simulation without 

feedforward 
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Figure 75: (a) Output voltage simulation with feedforward 
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Figure 75: (c) Output filter inductor current simulation 

with feedforward 
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Figure 75: (d) Control voltage simulation with 

feedforward 






(b) 


Figure 76: Measured output voltage ripple without (a) and 

with (b) feedforward 

(a) Y — 0.5v/div 


(b) Y— 0.5v/div 


X — 0.5 msec/div 
X — 0.5 msec/div 
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seen that the regulator loses regulation momentarily and the 
output voltage drops down to around 18.3V before regulation 
is regained. It should be noted when comparing the analyti-: 
cal plots without feedforward (Figures 74) with those with 
feedforward (Figures 75), the transient waveforms are almost 
identical. This is expected, since during this large signal 
transient the regulator has momentarily lost its voltage re- 
gulation. Therefore the ability of feedforward compensation 
via duty cycle modulation is also lost momentarily. The re- 
gulator is operating in an open loop fashion (under minimum 
off time control). This can be verified by observing the 
voltage waveform input to the regulator. 

The input filter capacitor voltage in Figures 73 (b) 
and 75 (b), which is also the input voltage to the regulator 
goes down close to 20V in both cases. For such a low input 
voltage to the regulator the duty cycle would have to be 
close to unity to maintain regulation, which is only possi- 
ble if the off time is decreased sharply. However a minimum 
Tqfp time of 5 microseconds is implemented for various mag- 
netic reset purposes and thus the regulator loses regulation 
till the input filter capacitor voltage rises. 

The analytical plots of the control voltage without and 
with feedforward both show the large disturbance in control 
voltage as a result of the loss of regulation. The plot with 


feedforward does show a somewhat smaller peak value of the 
overshoot than the one without feedforward but this differ- 
ence is about 9%. Negligible difference is noticed between 
the analytical plots of the inductor current without and 
with feedforward. The large drop in current as a result of 
loss of regulation is noticed in both plots, as is the over- 
shoot when regulation is regained. 

The analytical results thus strongly support the exper- 
imental data and show that the regulator momentarily loses 
its control function- It is also confirmed through simula- 
tion and testing that the feedforward does not have any de- 
trimental effect under large transient conditions.. 



Chapter VIII 


DISCRETE TIME DOMAIN STABILITY ANALYSIS OF BUCK 

REGULATOR 

An important concern is the interaction of an input filter 
with a switching regulator which often results in system 
instability. As described in detail in Chapter 2, the input 
filter interaction with the regulator control loop causes a 
reduction in the loop gain which may result in system inst- 
ability. When a switching regulator is acquired as a 'black 
box' for use in a system and an input filter is put in ex- 
ternally, the improper choice of filter parameters can cause 
the system to be unstable [4,5]. Work presented earlier 
[4,5] showed that the input filter design should be incorpo- 
rated in the regulator design itself in order to avoid loop 
instability and other input filter related problems. 

In this Chapter the stability of a buck regulator is 
examined by varying the input filter parameter values. Sys- 
tem instability is predicted analytically, and backed up 
with experimental data. It is then shown that the addition 
of the feedforward control enables one to stabilize the sys- 
tem. This is verified both analytically and experimentally. 
One can conclude that the addition of feedforward , in fact, 
isolates the regulator from the input filter; any input fil- 
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ter can be used with the buck regulator without causing loop 
instability. 

8.1 ANALYTICAL INVESTIGATION OF STABILITY 

An exact time domain analysis is performed to calculate the 
eigenvalues of the system with a single stage input filter; 
it was used to study how the eigenvalues move with a change 
in the input filter parameter (inductance). To correlate the 
'exact' discrete time domain analysis results with the 'ap- 
proximate' continuous time average analysis presented in 
Chapters 2 through 5, a computer program which calculates 
the closed loop poles of the continuous system was also 
written. The closed loop poles as calculated using the pro- 
gram track closely the eigenvalues derived from the exact 
time domain analysis. System instability occurs for some va- 
lues of filter inductance but the addition of feedforward 
eliminates the problem. 

8.1.1 Calculation of Eigenvalues . 

The method used in calculating the eigenvalues has been de- 
scribed earlier [15,16,17]. The state equations of the sys- 
tem are written for the buck regulator of Figure 34 . The 
state variables describing the system without feedforward 

x(l) = ij^j (input filter inductor current) 


are : 
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x(2) = v_ (input filter capacitor voltage) 
x(3) = ij^ (output filter inductor current) 
x(4) = Vq (output filter capacitor voltage) 
x(5) = 
x(6) = e^ 

As discussed in Chapter 7 is the voltage at the capacitor 
C *2 in the feedback loop while e^, is the control voltage (in- 
put to pulse modulator) . Except for the order in which they 
appear, these are the same state variables as used in Chap- 
ter 7, where it is noted that the system is completely de- 
scribed by the above six state variables. 

The state equations without feedforward are written as 
in Chapter 7: 

During (SI on S2 off) 

( 8 - 1 ) 

X = Fl*x + Gl»u 
During T__ (SI off S2 on) 

( 8 - 2 ) 

X = F2*x + G2*u 

The matrices FI and F2 are identical to the ones in Chapter 
7 except for the following reordering; the third rows of 
equations (7-2) and (7-8) are moved to the sixth place with 
the other rows being moved up one place to get FI, F2 of 
equations (8-1) and (8-2), while the third columns of equa- 
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tions (7-2) and (7=8) are moved into sixth place with the 
other columns being moved up one place in a manner similar 
to the rows. Matrices G1 and G2 of equations (8-1) and (8-2) 
are identical to the ones in Chapter 7 except that the third 
row of Chapter 7 is moved into sixth place with the other 
rows being moved up one place. The vector u of equations 
(8-1) and (8-2) is defined exactly as in Chapter 7. 

Figure 77 serves to establish notation regarding time 
instance t^ etc . and the solutions to the state equations 
are written as in Chapter 7: 

During 

x(t^ + T) = <|ip • x(tj^) + Dp « u 

During ,3.3, 

+ T) . . x(t^ i 

* • 2 

The matrices , Dp euid Dj^ are defined, as is time step 

T, exactly as in Chapter 7. 

From (8-3) the following is derived [15,16,17] : 

• ♦p • 2(tK> + • Dp • H 

* • 2 

For the constant volt- sec. (VjTon) duty-cycle control mode 
used in the stability investigation, the Tqjj is fixed if Vj 
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is fixed euid Tq^^, is determined by the point at which e^ in- 
tersects the threshold voltage : 

6 4 

E» = I <l»p(6,I) • x(I) + I D (6,J) • u(J) (8-5) 

^ 1=1 * J=1 ^ 

K 

^OFF thus a function of (8"4) is a nonlinear 

equation [15,16,17]. Equations (8-1) - (8-5) thus represent 
exactly the nonlinear buck regulator system without feedfor- 
ward . 

For the purpose of stability analysis the steady state 
or equilibrium state is necessary [15,16] and this is der- 
ived as 

2£(tK+l) = 2£(tK) = X* 

( 8 - 6 ) 

X* = ' X* + V°F 

The nonlinear system is linearized for a small perturbation, 
around the steady state operating point [15,16,17]: 


5x(tj^) = x(tj^) “ X* (8-7) 

and 

= K • 57 [*p • i"=K> * °P • h]} sxfV (3-8) 

X* 
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Equation (8-8) is derived assuming that the small perturba- 
tion in the state does not change the value of [15,16]. 

The expression in the curly brackets of equation (8-8) is 
us^d to investigate stability: 

• 5x(tj^) , (8-9) 

Equation (8-9) thus describes a linear system which will be 
stable if and only if all the eigenvalues of are abso- 
lutely less than unity — 

U^(,(,)( <1 i = 1,2,3,... (8“10) 

or in other words if all the eigenvalues X. lie inside the 
unit circle [15,16]. 

The buck regulator system with feedforward is described 
by a set of seven state variables as in Chapter 7. The nona- 
daptive feedforward circuit of Figure 67 was used in the 
analysis and the state variables then are: 

x(l) = ij^i 
x(2) = v^.^ 

x(3) = v^ (voltage at capacitor ) 

x(4) = (8-11) 

x(5)=v^ 

x(6) = e^ 
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x(7) = 

The state equations during and are written as: 

During Tqjj 

i = FI. 2 + Gl.ii (8-12) 

During Tqpp 

X = F2.X + G2.U (8-13) 

Matrices FI, F2, G1 and G2 of equations (8-12) and (8-13) 
are similar to the ones in Chapter 7 and are defined exactly 
as the ones in equations (8-1) and (8-2). The solution of 
the state equations, the definition of the equilibrium state 
and the matrix tjf are identical to the case without feedfor- 
ward. 

A computer program was written to calculate the matrix 
and its eigenvalues both without and with feedforward. It 
can be easily seen that by changing the single stage input 
filter parauneters a set of eigenvalues can be calculated 
from which inferences regarding system stability can be 
drawn. 

8.1.2 Description of the Program Used in Calculating 
Eigenvalues 

A flow chart for the program is given in Figure 78 and a 
computer listing of both programs is included in the appen- 


dix. 


There are essentially three steps involved in the 
calculation of eigenvalues. 

First the approximate steady state values are computed 
[15,16,17], using an approximate value of the off time TQ^p; 

"^OFF ■ *^0N ‘ (8-14) 


Equation (8-6) is used to calculate the approximate steady 
state values x* with ISML routine LEQT2F being used to solve 
the system of linear equations. It is to be noted that equa- 
tion (8-6) is not valid for the last state variable e since 

w 

both ^jj(6,6) and ^^ 6 , 6 ) equal unity, and equation (8-5) is 
used to calculate e^. 

c 

The aporoxiraate steady state calculated is used to find 
the exact sceady state values [15,16]. The best way of det- 
ermining the exact steady state is to determine the exact 
off time by iterative linearization ( Newton's method) on 
the cycle- to-cycle matching condition for e^, [15,16,17]. The 
vector ZT is defined : 


K 


^ = x(tK + Tjpp) 


ZT = 4'p • x(tj^) + Dp • u 


(8-15) 


with ZT(6) = 

The cycle-to-cycle matching condition is expressed as: 
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SMAT = x(6) 


5 

I <fr„(6,I) • 2T(I) + E_ 

1=1 1 

4 

+ I D„(6,J) • u(J) 

J=1 


(8-16) 


The value of SMAT is calculated using the available value of 

Topp; if SMAT is smaller than a certain error tolerance then 

the value of T___ used is accurate. If SMAT exceeds the to- 
OPP 

lerance then Newton' s method is applied to calculate a bet- 
ter estimate of TQpp. The steps involved in the calculation 
are: 


1. Calculate 2T and SMAT. 

2. If SMAT is not less than specified tolerance, then 
perturb T^pp by the amount DT^pp. 

3 . Calculate values of state variables corresponding to 
the perturbed value of off time, using equations 
(8-5) and (8-6). 

4. Calculate SMATN. 

5. Calculate a better estimate of using Newton's 

iteration [15,16,17]. 

6. Calculate a better estimate of the state variables. 

7 . Calculate SMAT and go to step 2 . 

The last step is to calculate the matix ‘\j; and its ei- 
genvalues using the exact steady state values computed 




Figure 78; Flowchart for stability analysis 
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above; Figure 79 is a flowchart for the calculation. The 
following are defined: 


FUN l (X* + Dx(I) ) 
FUN 2 (X* - Dx(I) ) 


• S + Op 

[♦p • 5 + Dp 



+ Dx(I) 
- DX(I) 


(8-17) 

(8-18) 


TEMP 2 = 


( FUN l - FUN 2) 
2 • Dx(l) 


( FUN l - FUN 2~ 

2 • Dx(6) 6X6 


(8-19) 


PSI = <|)jj « TEMP 2 


( 8 - 20 ) 


TEMP2 is thus the matrix of partial differentials used in 
equation (8-8). 

The matrix TEMP2 is evaluated one column at a time. For 
the I'th column a perturbation DX(I) is made in X(I). Two 
new sets of steady state values XA and XB are defined, as is 
shown in the flowchart. For the new set of values XA it is 
necessary to calculate cui exact value of off time TFFCl in 
order to obtain ^ , Dp. The error function ZETA is defined: 


6 4 

ZETA = I <()„(6,I) • XA(I) I D_(6>J) u(J) - E, 

1=1 ^ J=1 ^ 


( 8 - 21 ) 



1102 


1003 


Figure 79: Flowchart for calculating matrix 
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Figure 79: continued 
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Figure 79: continued 
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Figure 79: continued 
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The value of ZETA is calculated using the available value of 
TFFCl, if it is within the specified tolerance then the va- 
lue of TFFCl is the value needed. If ZETA is not within to- 
lerance then Newton's iteration [15,16] is used to obtain a 
better estimate of TFFCl, as is evident from the flowchart. 
The values of ^ , Dj, are obtained and then the vector FUNl 
is evaluated. 

A similar procedure is used to obtain the exact off 
time TFFC2 for the other set of values After calculation 
of FUN2 an estimate for the I'th column of TEMP 2 is ob- 
tained. 

In the calculation of the partial derivatives it was 
observed that a better estimate of the derivatives can be 
obtained using the Richardson's extrapolation method [18]. 
The derivative f' (a) is sometimes approximated as: 


A(h) 


_ f (a + h) - f (a - h) 
2h 


( 8 - 22 ) 


where h is the step size. 

A better estimate may be obtained by evaluating A(h) and 
A(rh): 


f (a + rh) - f (a - rh) 
. 2rh 


A(rh) 


(8-23) 
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where r is usually 0.5 [18]. 

Combining A(h) and A(rh) thus: 

B(h) = Mrh) - r^A(h) (8-24) 

1 - 

gives a better estimate, B(h), of the derivative f'(a) [18]. 
Thus by combining the value of A(h) obtained from equation 
(8-22) with two values for the step size h, it is possible 
to get a more accurate estimate of the derivative-. 

In the program four values of step size were used. The 
perturbation DX(I) is assigned four values and for each va- 
lue the vectors FUNl and FUN2 are calculated. This results 
in four estimates for the derivatives in the I'th column of 
TEMP2 . Combining these four estimates as in the program 
[18] gives an accurate estimate of the I'th column. 

The whole procedure is then repeated for the next co- 
lumn, thus an accurate set of values for the derivatives is 
obtained. Multiplying TEMP 2 by results in the matrix 
[15,16,17], Finally the eigenvalues X of the matrix are 
obtained using the IMSL routine EIGRF. 

The eigenvalues with feedforward are calculated in an 
identical fashion, the only real difference being the change 
in the order of the system from six state variables to sev- 


en. 
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8.1.3 Calculation of the Closed Loop Poles of the System 
The closed loop poles of the system of Figure 34 can be cal- 

A A 

culated. The expression for the closed loop gain can 

be easily written using Figure 33 for the case without feed- 
forward: 




(8-25) 


with Fjj , etc. as defined in Chapter 5. 

Equation (8-25) can also be used to calculate the closed 
loop poles of the system without input filter; all that is 
necessary to do is to set H = 1 and 2=0, as discussed ear- 


lier. 

Expressions for A etc. are svibstituted into equa- 
tion (8-28) for the case without input filter. The resulting 
expression in the denominator is then a polynomial of the 
fourth degree : 


Pj^(s) 


= S (K.^KgLC) + 

S^(K.7KgKg + K.7C^LC + K^K^CR^) + 
S^(K^C^Kg + K.yKg + KgK2CR^ + 

+ KgnL C Rj^) + 

S(K. 7 C[ + KgK 2 + KgKj^C R^ + nL Kg) 

+ KgK^ 


(8-26) 
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The constants etc. are defined as functions of the 
regulator parameters like inductance L etc. The appendix 
contains a listing of the program used and the constants are 
defined there. An IMSL routine ZPOLR is used to calculate 
the roots of Pj which are the closed loop poles of the sys- 
tem without input filter. 

For the case with input filter the same equation, 
(8-25), is used and sxibstitution leads to a sixth degree po- 
lynomial in the denominator: 


P2<s) - S^K^Kjig) + 

S“ (K^C^Y, + K^KgY^g + KjKgY^ + K^nLYj 
+ K^KjYg) 4- 

K^nLYg + K5K4Y.7) + 

(KtCI^U ^ ='7«6’'l2 ’'3 ''y’^3 

+ KgnLY.^ + K5K4Y4) + 

S(C^K.^Y^2 + ^3^2^4 


(8-27) 


+ K3K,Y^ 
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The constants Yj etc. used in equation (8-27) are again de- 
fined in the program. The same routine ZPOLR is used to ob- 
tain the six closed loop poles. 

■'8.1.4 Location of Eigenvalues and Poles as a Function of 
Input Filter Parameters . 

The computer programs written were used to locate closed 
loop poles auid eigenvalues as a function of input filter in- 
ductance. The buck regulator of Figure 34 was used and the 
parameters of the regulator were as defined in section 6.1 
except for the following changes; 

Vj = 25 V Rj 3 = 200 Kohms 

C 2 = 100 pf Rl “ 10 ohm 

The values of Rj^j and C 2 were changed so that the compensa- 
tion loop is largely ineffective, (7,8). This was done in 
order to facilitate making the measurements discussed in 
section 8.2. It was found difficult to achieve instability 
in the high loop gain system used in making the measure- 
ments. Making the compensation loop largely ineffective re- 
duced the loop gain so that with a large input filter the 
system was made unstable. 

A single stage input filter with the following parame- 
ters was used: 

Rj^^ = 0.2 ohm LI variable Cl = 220 microF 


V. 
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First the eigenvalues and closed loop poles without input 
filter and feedforward were calculated using the eigenvalue 
and the closed loop programs. Figures 80 and 81 are the 
plots of the eigenvalues and closed loop poles respectively. 
It is noticed that there are four poles and four eigenvalues 
since without input filter the system is of the fourth ord- 
er. 


The eigenvalues and poles are: 
Zj = 0.844 * 10"^ 

= 0.12 

Zj =0.98 + jO.103 
Z^ = 0.98 - j0.103 


s^ = -55,000 + jO 
s I = -50, 000 + jO 
S 3 = -346 + J2240 
S 4 . = -346 - j2240 


The eigenvalues and the closed loop poles are related 


[ 15 ] : 

sT 

z = e**^ (8-28) 

s = a + jo3 

Thus 

|z| = e'^'^ (8-29) 

^z = o)T 


Knowing the location of the closed loop poles it is thus 
possible to locate the eigenvalues if a value of T is known. 
In this case T is the switching period and its value was 
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foiind to be 41.85 microsec.; thus equations (8-28) and 
(8-29) could be used to crosscheck the results from the ei- 
genvalue program once the closed loop poles are known. Equa- 
tion (8-29) was used to calculate the eigenvalues and the 
calculated values were found to be very close to the values 
from the eigenvalue program. The two poles on the real axis 
correspond to the two eigenvalues close to the origin while 
the two complex poles result in the two complex eigenvalues. 

Next the single stage input filter was put in and the 
eigenvalues and closed loop poles calculated using the res- 
pective programs. Figures 82 and 83 are the plots of the ei- 
genvalues and the closed loop poles respectively. The values 
of the filter inductance LI used were -- 

50, 150, 325, 450, 650, 800, 1000, 1425 and 1800 microH 
Table 1 lists the eigenvalues and closed loop poles, for 
different values of LI. Starting at 50 microH it is seen 
that the complex eigenvalues move to the right, closer to 
the unit circle as inductance is increased. The two eigen- 
values on the real axis and Z 2 , , which correspond to 
and z 2 ^ on Figure 80 , do not move with change in inductance. 
The complex eigenvalues with a real part of 0.98 z^ and z«, 
which correspond to z^ and z ^ on Figure 80 , are pushed out 
of the unit circle for a value of LI = 800 microH. For va- 
lues above 800 microH the eigenvalues do not move much; the 
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Figure 80: Eigenvalues without input filter 
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\ 


Figure 81 : Closed loop poles without input filter 
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real part stays at around 1.003 for the ones outside the 
linit circle and around 0.97 for the ones inside. This may 
be used to explain the bunching up of the eigenvalues for 
large values of inductance. Thus the eigenvalue plot of Fig- 
ure 82 predicts that the system will be unstable for values 
of filter inductance above 800 microH. 

The plot of the closed loop poles. Figure 83 , shows 
similar behavior. The two poles on the real axis s and s 2 , 
are not changed as inductance increases, however two complex 
poles S 5 and s^ move into the right half plane, indicating 
instability. The poles move across the origin at 800 microH 
which checks well with the eigenvalue plot. Equation (8-29) 
was used to calculate the eigenvalues from the closed loop 
poles and the calculated eigenvalues agreed closely with the 
results from the eigenvalue program for all the values of LI 
that were used, thus providing a crosscheck. 

The other pair of complex poles s^ and in Figure 83 ini- 
tially move closer to the right half plane, then move away 
from it for values of L 1 greater than 450 microH. This be- 
havior should be reflected in the eigenvalue plot for eigen- 
values and Z|., however it is noticed than in equation 
(8-29) the value of T is so small that for a change in the 
real part from -300 to -500 the change in the eigenvalues 
will be about 1 %, which is too small a change to show up in 
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Eigenvalues 


“0.003 -i” jO 
0.1233 + jO 
0.971 ± 
j 0.132 
0.9897 i 
j 0,1026 


“0.003 + jO 


Closed Loop Poles 


-55,045 -t- jO 
-50,005 + jO 
-323 > 
j 3165 
-214 * 
j 2237 


-55,045 + jO 


0.1233 + jO 

-50,005 + 

0.97 ± 

-419 

j 0.1162 

j 2684 

0.9975 ± 

-49.8 ± 

j 0.0954 

j 2189 

-0.003 + jO 

-55,045 + 

0.1232 + jO 

-50,005 + 

0.97 i 

-483 ± 

j 0.111 

j 2546 

0.9996 ± 

43.6 ^ 

j 0.089 

j 2069 
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TABLE 1 
continued 



Eigenvalues 

Closed Loop Poles 

LI = lOOOyu-H 

-0.003 + jO 
0.1233 + jO 
0.971 ± 
j 0.1078 
1.0002 ± 
j 0.082 

-55,046 + jO 
-50,005 + jO 
-513 ± 
j 2455 
98.3 ± 
j 1909 

LI = 1425^ H 

-0.003 + jO 
0.1233 + jO 
0.972 ± 
j 0.1052 
1.003 + 
j 0.071 

-55,046 + jO 
-50,005 + jO 
-517 + 
j 2371 
133 ± 
j 1650 


-0.003 + jO 

-55,046 + jO 


0.1233 + jO 

-50,005 + jO 

LI = 1800 yUH 

0.972 ± 

-511 ± 


j 0.1055 

j 2339 


1.004 ± 

141 ± 


j 0.063 

j 1487 
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Figure 82; Eigenvalues with input filter, without 
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Figure 83: Closed loop poles with input filter, without 

feedforward 
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the eigenvalue plot. The bunching up of the closed loop 
poles for higher values of LI is also noticed; this provides 
another crosscheck with the eigenvalue plot. 

Thus it is seen that the eigenvalues and the closed 
loop poles match both qualitatively, as examination of Fig- 
ures 82 and 83 show, and quantitatively through the use of 
equation (8-29) to crosscheck t^e eigenvalues. It is also 
noticed that system instability is predicted for values of 
LI above 800 microH. 

Lastly the eigenvalues were calculated for the case 
with feedforward using the program. The nonadaptive feedfor- 
ward circuit of Figure 67 was used with the following param- 
eters: 

Cj, = 27 microF Rjj = 5100 ohms 
variable from 210 ohm to 300 ohm. 

Figure 84 shows the plot of the eigenvalues with feedforward 
for the following values of with the value of LI kept 
fixed at 1425 microH; 

R =210, 220, 230, 235, 240, 260, 270 and 300 ohms 

X ^ 

The value of Rj^^ was changed to see the effect of feedfor- 
ward on the eigenvalues. Table 2 lists the eigenvalues for 
different values of Rfji • It is noticed from Figure 84 that 
there are seven eigenvalues and that the three on the real 
axis , Z 2 and Z 3 do not change as R ^2 changes in value. 



As is increased one pair of complex eigenvalues and 
Zg move outside the unit circle while the other pair z® and 
z^ move inside. For the value of LI used the closed loop 
poles and eigenvalues without feedforward clearly predict 
instability, whereas with feedforward it is seen that the 
system will be stable for some values of Rj 2 . • regu~ 
lator under consideration the value of is calculated to 
be 220 ohm from the feedforward design as obtained in Chap- 
ter 5. A value less than the design value causes the eigen- 
values Zg and z^ to be outside the unit circle while it is 
noticed that a value higher than the design value also caus- 
es instability in the system. For the design value of 220 
ohms all the eigenvalues are inside the unit circle thus 
clearly indicating that the addition of feedforward stabi- 
lizes a system that was \instable due to the interaction bet- 
ween the input filter and the regulator control loop. 

As R^^ is increased further the real parts' of the com- 
plex eigenvalues do not change much but eigenvalues z^ and 
Zg do move outside the unit circle. For R^^ between 220 and 
260 ohm all the eigenvalues are inside the unit circle, with 
the optimum value being around 230 ohm. 

In order to see clearly the effect of changing the 
feedforward loop gain the open loop gain and phase margin 
were calculated for various values of R^^ • The parameters of 
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Figure 84: Eigenvalues with feedforward 
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the buck converter used in the calculation were the same as 
used in the eigenvalue calculation, with LI = 1425 microH. 
Figure 85 (a) - (f) shows the plots of the open loop gain 
and phase margin for three values of Rj 2 • Th® solid lines 
on the figures are the plots without feedforward and with 
feedforward ~ ohms), whereas the triangles repre- 
sent values with either lower than or higher than the 

design value of 220 ohm. Exact compensation of the input 
filter interaction is achieved with ohm, with any 

value other than 220 ohm the compensation is not complete. 
Decreasing the feedforward loop gain (210 ohm) results in a 
loss of open loop gain at the filter resonant frequency. In- 
creasing the feedforward loop gain results in a loss of 
phase margin at the resonamt frequency. For a value of 240 
ohm the open loop gain obtained is close to the ‘ gain with 
220 ohm but the phase margin is decreased. Increasing the 
feedforward gain further (300 ohm) makes the situation worse 
as Figure 85 (f) shows. The open loop gain in Figure 85 (e) 
is decreased to around 8 db at a frequency of about 300 Hz 
while at that frequency the phase margin is almost zero, in- 
dicating that the system is close to instability. This is 
confirmed on examination of Figure 84 which shows that the 
system is unstable for R^^ = 300 ohm. It should be noted 
that the open loop gain with an input filter and I^ 2~220 ohm 
is identical to that without input filter. 
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It is thus concluded that exact compensation is ob- 
tained at the design value of the feedforward loop gain. In- 
creasing or decreasing the feedforward gain from the design 
value results in incomplete compensation. There may be no 
compensation achieved at ail if the feedforward gain is far 
different from the design value. 

It is noticed in Figure 84 as in the other eigenvalue 
plot that the eigenvalues do not seem to move very much, for 
example in Figure 84 the real part of the pair and z j 
changes from about 0.9885 to 0.9975 while from Figure 82 it 
changes from about 0.98 to 1.004. It is noted however that 
from Figure 83 the above change in real part corresponds to 
a change in th<= real part of the poles from about -323 to 
+140. This may be explained by noting as before that T in 
equation (8-29) is fairly small so this change in real part 
will correspond to a very small change in the real and imag- 
inary parts of the eigenvalues. 

From the results presented in this section it is seen 
that the eigenvalue program is a powerful tool for the ana- 
lysis of stability. . It is also seen that for the regulator 
used, a high value of input filter inductance will cause 
instability, but the addition of feedforward restores sta- 
bility. Experimental data that strongly support this analyt- 
ical prediction are presented next. 


FREQ (HZ) 


h3 


Figure 85; (a) Open loop gain: calculated values without 

feedforward, with =220 ohm and Rfz =210 ohm 




FREQ ( HZ ) 


(b) Open loop phase margin: calculated values 
without feedforward, with Rj2=220 ohm and 
Rj 2,=210 ohm 


Figure 85 : 


FREQ (HZ) 


Figure 85: (c) Open loop gain: calculated values without 

feedforward, with Rfo=220 ohm and R.«=240 ohm 
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FREQ (HZ) 


Figure 85: (d) Open loop phase margin: calculated values 

without feedforward, with Rf.»=220 ohm and 
Rj^=240 ohm 




3 


• calculated values without 
=220 ohm and =300 ohm 

f 2f 



O i 

I .E+2 


1 .E+3 


I .E+4 


FREQ (HZ) 


Figure 85: (f) Open loop phase margin: calculated values 

without feedforward, with R,,=220 ohm and 
R- =300 ohm 

i2 


Eigenvalues 
0-00116 + jO 
0.1282 + jO 
= 210 ohm 0.9998 + jO 

0.9974 i j 0.0738 
0.9883 ± j 0.101 
0.00017 + jO 
0.1282 + jO 

= 220 ohm 0.9998 + jO 

0.9963 A j 0.0742 
0.9907 ± j 0.1006 
0.00122 + jO 
0.1282 + jO 

Rj^ = 230 ohm 0.9998 + jO 

0.9955 ± j 0.0743 
0.9922 ± j 0.10121 


/ 
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= 235 ohm 


Rj^ = 240 ohm 


R,^ = 260 ohm 

I 2. 


TABLE 2 
continued 

Eigenvalues 
=0.00025 + jO 
0.1282 + jO 
0.9998 + jO 
0.9952 ± j 0.0742 
0.9931 + j 0.1006 
-0.00011 + jO 
0.1282 + jO 
0.9998 + jO 
0.9948 + j 0.0743 
0.9936 + j 0.10049 
0.00039 + jO 
0.1283 + jO 
0.9998 + jO 
0.9946 + j 0.0744 
0.9941 + j 0.1011 
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TABLE 2 
continued 



Eigenvalues 

Rf2, ~ 

0-00023 + jO 
0.1283 + jO 
0.9998 + jO 
0.9933 + j 0.0742 
0.9964 + j 0.1008 

= 300 ohm 

0.00029 + jO 
0.1285 + jO 
0.9998 + jO 
0.98952 + j 0.0724 
1.0054 + j 0.1036 
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8.2 EXPERIMENTAL VERIFICATION OF STABILITY ANALYSIS . 

In this section experimental data pertaining to the stabili- 
ty of the system are presented. The buck regulator of Figure 
34 was used with the parameters as defined in section 8.1.4 
with the exception that the input filter inductance LI had a 
constant value of 1.425 mH. The adaptive feedforward circuit 
of Chapter 5 was used, however it is noted that if the input 
voltage is kept constant, as it was in the following mea- 
surements, the nonadaptive and the adaptive feedforward cir- 
cuits provide the Seune gain. 

The regulator was set up with the HP network analyser 
to measure open loop gain exactly like in Chapter 6 [7,8]. 
The open loop gain and phase margin were measured without 
input filter, with and without feedforward. Figure 86 shows 
the measured values of the open loop gain and phase margin. 
It can be seen that without the input filter the gain dips 
down around 320 Hz — this may be explained by noting that 
the values of R^^ and C^ are such that the compensation loop 
[7,8] is largely ineffective. Thus the open loop gain and 
phase margin correspond to a value of «< = 0.355 a^.. shown in 
the literature [7,8]. 

The. addition of the input filter causes the gain to 
fall below 0 db and when the gain is close to 0 db the phase 
margin is only about 10 degrees at a frequency of about 250 




Figure 86: 


(a) Open loop gain: measured values without 
input filter, with (•) and without . feedforward 
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Fic&4)oeMcr (Hz) > 


Figure 86: (b) Open loop phase margin: measured values 

without input filter, with («) and without 
feedforward 
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Hz. The system is thus marginally stable and a small signal 
250 Hz sinusoidal disturbance is injected into the loop to 
induce unstable operation and it was observed that the out- 
put voltage has oscillations. Figure 87 shows the output 
voltage ripple -- it is clearly seen that the oscillations 
are around 1.2 v in magnitude. It is concluded that the sys- 
tem is clearly unstable. Analysis presented earlier showed 
that the syatem would be unstable for LI = 1425 microH while 
experimental results show a marginally stable system, thus 
lending support to the analysis. 

Feedforward was added to the circuit and measurements 
of the open loop gain and phase margin were made. It is seen 
clearly from Figure 85 that the addition of feedforward re- 
moves the instability — this is conf ironed by recording the 
output voltage ripple. Figure 88 shows the output voltage 
ripple and it can clearly be seen that the oscillations that 
were present earlier without feedforward have been eliminat- 
ed, making the system stable. The 250 Hz oscillation in Fig- 
ure 88 is the injected disturbance. The analytical predic- 
tion that feedforward eliminates instability caused by input 
filter interaction is thus strongly supported by experimen- 


tal data. 
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(b) 


Figure 87: Output voltage ripple without feedforward: 

Y:0.2v/div X:5 msec/div. (b) Y:Q.2v/div 
msec/div. 


(a) 

X:10 
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8.3 CONCLUSIONS 

It is seen that the eigenvalue analysis program is a power- 
ful tool for the auialysis of stability. Using a buck regula- 
tor system a stability analysis was performed for various 
values of input filter inductance and it was observed that 
system instability was predicted for large values of induc- 
tance. A valuable crosscheck was provided by a program that 
calculated the closed loop poles and it was observed that 
the closed loop poles and the eigenvalues track very closely 
for all the values of inductance used, including those va- 
lues for which system instability is predicted. 

The addition of feedforward removes the instability at 
large values of inductance as the eigenvalue analysis pre- 
sented shows. It is noticed that complete compensation of 
the input filter interaction is obtained for the design va- 
lue of feedforward gain; any other value of feedforward gain 
results in incomplete compensation. It was demonstrated both 
from the eigenvalue analysis and the open loop gain and 
phase margin plots that a too large or too small value of 
feedforward gain would not be able to remove the instability 
caused by the input filter interaction. 

Experimental measurements were made to confirm the ana- 
lytical results presented. It was noticed that without feed- 
forward the output voltage has large magnitude oscillations 


265 


indicating system instability- These oscillations occur at a 
value of input filter inductance for which the eigenvalue 
analysis also predicts instability thus tying in the analy- 
sis with experimental measurements. The addition of feedfor- 
ward (with the gain set to its design value) completely re- 
moves these oscillations as measurement data show. Since the 
eigenvalue analysis also predicted stability with the addi- 
tion of feedforward, it is thus concluded that the addition 
of feedforward does result in removing the system instabili- 
ty caused by input filter interaction with the regulator 
control loop. 


Chapter IX 


EXTENSIONS OF THE CONCEPT OF INPUT FILTER 
COMPENSATION TO OTHER CONVERTER TYPES 

The concept of feedforward loop compensation to mitigate 
converter — input filter interaction is extended to several 
other cases, three of which are discussed in this chapter. 

9.1 EXTENSION TO OTHER TYPES OF REGULATOR CONTROL 

The concept of feedforward compensation for the input filter 
has been developed with reference to a buck regulator. Ex- 
tensive experimental data, together with analysis, confirm 
the validity of the concept. This was. done exclusively on a 
buck converter with multiple- loop control and constant 
volt-sec duty cycle modulation. Use of other types of con- 
trol are discussed in this section. 

9.1.1 Feedforward compensation for multiple- loop control 
using other types of duty- cycle control 

The feedforward design was presented in Chapter 5 where it 

was noted that the gain of the feedforward loop depends on 

the transfer function of the pulse modulator ,Fj^ . The 

transfer function Fj^ is dependent on the type of duty-cycle 

control used [2,8]. In this dissertation the constant volt- 

sec. type of duty-cycle control was used. However, it is 
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easy to see that for other types of duty-cycle control, for 
example for constant frequency control or constant off time 
control, the feedforward design procedure outlined in Chap- 
ter 5 is still valid. The feedforward design process pre- 
sented earlier is general since it treated the pulse modula- 
tor in terms of its transfer gain only. Thus other types of 
duty-cycle control can be easily incorporated in the feed- 
forward compensation after the transfer gain of the duty-cy- 
cle control is identified. 

As an example for multi-loop constant frequency control 
the following pulse modulator gain is identified [8]: 


where 



nM 


M = V_(1-2D)T„ 
I P 


(9-1) 


In equation (9-1) Tp is the constant switching period. Sub- 
stitution into the feedforward loop gain expression from 
Chapter 5 results in; 

^ -P^nM 

For a constant supply voltage D and M would be constants and 
thus the feedforward gain would be a constant. An adaptive 
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feedforward circuit that tracks the supply voltage can be 
realized using equation (9-2). 

For multi-loop constant off time control is identi- 
fied as [8] : 



nM 


V T 
I OFF 


(9-3) 


where TQpj,is the constant off time. Sustitution into the 
feedforward loop gain expression results in: 


"D^nV T 

V^2R.C,' 
O 4 1 


(9-4) 


For a constant supply voltage D would be a constant and thus 
the gain would be a constant. An adaptive feedforward cir- 
cuit that tracks the supply voltage variations can be real- 
ized from equation (9-4). 


9.1.2 Feedforward compensation for single- loop control 
A multi-loop feedback controlled buck regulator was used in 
the feedforward design discussed in Chapter 5. However, it 
is noted that in the design procedure the feedback loop 
transfer functions played no part in the feedforward design. 
This is evident since the feedforward gain is independent of 
the feedback loop transfer functions. 
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For single loop feedback control the regulator can be 
represented by a generalized block diagram as shovm in Fig- 
ure 89 . The feedback loop has as its input the output vol- 
tage ripple and it can be represented as Fg, as shown in 
Figure 89. Expressions for Fg and % are given in the litera- 
ture [8], and it is noted that the transfer function es- 
sentially processes the small signal variation in the output 
voltage and feeds it to the pulse modulator exactly like ^DC 
of Figure 28 . A feedforward loop that will sense the input 
filter state variedDles, process this information and feed it 
to the pulse modulator can thus be designed in a manner ex- 
actly as presented in Chapter 5. The transfer function F^^ 
would be a constant [ 6 ] and would be a part of the feedfor- 
ward gain, exactly as it is in the design presented. 

As an ex 2 unple the following is identified for single 
loop constant frequency control [6]: 


'■m = 


(9-5) 


where A^, Tp are constants. 

The error processor consists of an amplifier and a lead-lag 
compensation network and its transfer function % is availa- 
ble in the literature [6]. 
sign of feedforward loop is relatively simpler as substitu- 


Since F„ is a constant the de- 
M 
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tion into the feedforward gain expression of Chapter 5 
shows : 


-V A 

C = Q °p 

"=2 „ 2 


(9-6) 


For constant supply voltage the gain is a constant. An 
adaptive feedforward circuit that tracks the supply voltage 
variations could be developed exactly as shown in Chapter 5. 

It is thus seen that an extension of the feedforward 
concept to single loop control is easily affected. 


9.2 FEEDFORWARD COMPENSATION FOR THE BUCK-BOOST AND BOOST 
REGULATORS 

Extension of the feedforward corrept to the two other types 
of regulators — boost and buck-boost can be obtained follow- 
ing the outline presented in Chapter 5. A feedforward design 
for the buck-boost regulator is presented. 

The state space model of the buck-boost regulator is 
shown in Figure 90. In the figure the feedforward loop has as 
its input the input filter state variables (inductor current 
and capacitor voltage). These two are rilultiplied by the 
transfer functions c^ and c^ whose properties have yet to be 
determined. This state space model is based on the model 
for the buck-boost power stage developed in Chapter 4. The 


feedback loops form the standardized control module (SCM) 
developed earlier [2,7,8] and discussed in Chapter 5. The 
pulse modulator has as its transfer function F_. [ 8 ] and has 
as its input information regarding the state verifies of 
both the input filter and the output filter. The total 
state feedforward/feedback error signal is thus used to mo- 
dulate the duty cycle of the switch for loop gain correc- 
tion. 

The generalized small signal model for the buck-boost 
regulator is developed next and used to write the trainsfer 

A A 

function v /v^ . 

0 * 

9-2. 1 Generalized small signal model for the buck-boost 
regulator 

The generalized small signal model is shown in Figure 91. The 
regulator is modeled according to the three basic functional 
blocks : power stage, error processor and duty cycle pulse 
modulator, as in Chapter 5. It is noticed that in this case 
the flux is a state variable, because it is continuous un- 
like the inductor current. The transfer functions £"3 / and 
Fpc play no part in the feedforward design process, as ear- 
lier in Chapter 5. 

The power stage transfer functions F^^ etc. are written 
as in Chapter 5, using the following equations: 
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■^31^*1 ^32^ = ^C1 

■*’ "^42^ “ ^L1 


Thus it is seen that 


and 


T s — ^ 
^11 ; 


d = 0 


HF 


11 


■11 A 
P 


12 


12 

A 


(9-8) 


(9-9) 


The transfer ftinctions T^,T^,Tj 2 and T^j are developed 
with d=0 and the others with Vj * 0. The starting point for 
the derivation is again the small signal equivalent circuit 
model for the buck-boost regulator developed in Chapter 4. 


Evaluation of transfer functions with d = 0. 

A 

The equivalent circuit model with d = 0 is shown in Figure 
92 In Figure 92 the input filter has been replaced by its 
forward transfer function H(s) and its output impedance 
Z(s), the expressions for H and Z being the same as in Chap- 
ter 5. From the equivalent circuit of Figure 92 the follow- 
ing are derived, assuming that = 0 : 



— ^ 
H»Lg-Di(l + SCR^) 

n"a 

s 


(9-10) 



275 



ERROR PROCESSOR 


Figure 91: Generalized small signal model 
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■41 


H»a, 


(^1 


■+^iLTM 


where 

*1 ■ •" °'"c ■" 

+ RgDD'R^ + D'^R^ 

A = a^ + + sCRj^) 

DN 


D, = 


N. 


, where D is the duty cycle. 


(9-11) 


»L * \ “ A 

Using the above equations the following can be derived: 


^11 - Di D'R^d + sCR^) 


21 


31 


LgDid -f- sCR^) 


N 


= a. 


(9-12) 


A/H - a^ 

^41 “ (Rj^^ + sLl) 

The ESR of the capacitor R^^j has been neglected in the deri- 
vation for the same reasons as outlined in Chapter 5, while 
equation (5-9) is used in the derivation of ^41- 
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Evalua'tion of the transfer functions with =0 . 

The small signal equivalent circuit with Vj =0 is shown 
in Figure 93 The following are derived from the equivalent 
circuit: 


V^(l + SCR^) 


■12 


22 


■32 


■42 




Z + - DsLj 


DD' A 


.2, 


V^(l + SCR^) (-D'Dj^ Z + D Rj. ” 


o 


li 


NgDD* R^A 

Vq( 1 + SCR^) Dj^^Z (DR + DsLg - D ' ^R^) +V^a^ (DR - D^^Z) 


DD^D'^Rj^A 


= “^32/ ^\l ^ 


(9-13) 


where A and a^ are as defined earlier and 




I = 


V, = 


V 

V 

o 

D 


(9-14) 


Using the above equations the following are derived: 


Vod + SCR^) 


f 2 2 

' D^ Z + D' R, - DsL 


12 


' \ S 


i 


DD' 




22 


V^(l + sCRj^) (-D'Dj_^Z + D'^Rj^ - DsLg + DV^a^j 

2 


N„DD 
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p _ sCRj^)Dj^^Z(DR + DsLg - D'^R^^) + V^aj^(DR - D^^Z) 

dd-2djEj_ 

F^2 “ ■*■ 

(9-15) 

where R = R_-*-DD* R^ (9-16) 

s c 

The feedback transfer fimctions and consti- 

tute the two loop SCM feedback control scheme reported ear- 
lier (2,8). The transfer functions are defined in the liter- 
ature [2,7,8], and it is noted that they do not play any part 
in the design process, as in Chapter 5. The pulse modulator 
transfer function Fj^ depends on the type of duty cycle con- 
trol used [7,8] and it will be seen that the feedforward 
loop gain depends on E^, as earlier in Chapter 5 . 

9.2.2 Design of the Feedforward Loop 

With Vj = 0 the following equations are derived from 
Figure 91 ; 

and after substitution 




32 


42 


(9-18) 
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In the absence of feedforward i.e. with 02=0 = 0 ^ it can 
be seen that 

= V^g sCR^){-Di^Z 4- D'\ - DsLs|f„ 

V DD^"S 

X 

The effect of peaking of Z is to cause a reduction in the 
gain of the transfer function and thus also in the open loop 
gain. It is also noticed that equation (9-18) shows the po- 
sitive zero term in the transfer function. 

With the addition of feedforward the effect of peaking 
of 2 could be avoided by a proper choice of feedforward loop 
gains C 2 and C 3 . Two approximations are made to facilitate 
the design: 

1. The effect of the positive zero is assumed to be neg- 
ligible at the frequency at which Z peaks. Thus it is 
assumed that: 

- DsLg = (9-20) 

This is not an unrealistic assumption since Z peaks 
at a much lower frequency. 

2. It is also assiimed that 

2 2 
DR - Z = Z 


( 9 - 21 ) 
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at the frequency at which Z peaks. Again this is not 
an unrealistic assximption given the low value of R. 
Using the above assximptions, the following choice is made: 


03 

<=2 


= 0 



Vm 


(9-22) 


Svibstitution into equation (9-18) leads to: 




DD' A 


DD'^Rj^D^A 


n 2 2 ) 

®i°°i‘ ”°i 

d 


f (9-23) 


whi :h leads to 


lo , * 3CS^)D-y„ 

r. ~ Da, 


(9-24) 


It is seen from equation (.9-24) that the transfer function 
is now independent of the peaking effect of Z. It is also 
noted that equation (9-24) is valid only at low frequencies 
where the two assumptions made earlier hold. The feedforward 
loop gains are very similar to the gains developed for the 
buck regulator, and are identical if Ng = Np. It is also 
noted from equations (9-15) and (9-18) that the gain of 
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would be fairly small at frequencies other than the resonant 
frequencies at which 2 peaks. Thus the addition of feedfor- 
ward would not affect, in any noticable manner, the open 
loop gain and phase margin at any frequency other than those 
at which 2 peaks. Thus the feedforward compensation is ef- 
fective for most input filter designs because the filter re- 
sonant frequencies are relatively low in comparison to the 
positive zero. 

9.2.3 Analytical and Experimental Verification 
The buck-boost regulator shown in Figure 94 was used to ob- 
tain verification of the feedforward design. A multi-loop 
standardized control module (SCM) type feedback control was 
used [2,7,8]. The parameters of the converter are as fol- 
lows: 

Power Stage Parameters : 

Vj = 20 V Vq = 28 V D = 0.5833 Pq = 28 watts 
L 3 = 220 microH C = 300 microF 

= 0.05 ohm Rl = 28 ohm (load) 

Ng = Np = 33 Rg = 0.087 ohm 

Pulse Modulator Parameters : 

-3 

M = VjTqi^ 0.5*10 v-sec. 
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Control Circuit Parameters: 


7v 

Rg = 47.5 Kohm 
n = 0.667 


= 43-2 Kohm 
R^ = 40 Kohm 
Cj = 5600 pf 


= 15 Kohm 
Rj =1.1 Kohm 
c'^ = 32000 pf 


Input Filter Parameters 


0.2 ohm 

Lj = 325 

microH 

Cl 

= 200 microF 

0.02 ohm 

L 2 = 116 

microH 

C2 

=20 microF 


R 3 = 0.075 ohm ( ESR of Cj^ ) 


The regulator was operated in a predetermined duty cycle 
control mode (constant Vj Tqj^ control) [2,7,8]. 

The open loop gain can be written from Figure 91 as 


G^(s) 



12 

(_££) + F F ( 
' A ' 3-^AC' 






(9-25) 


It is noted that in deriving the above equation c^ is set 
equal to zero since that is the design value. Equation 
(9-25) was used to plot the open loop gain and phase margin 
without input filter, with and without feedforward. The 
transfer functions F^^ , F^ , F^^ and Fj^ constitute the feed- 
back and are defined in the literature [2,7,8], while 
F and F were defined earlier in this chapter. Equation 


31 


(9-25) is also valid for the case without input filter, all 
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that is necessary to do is to set Z = 0, while for the case 
without feedforward c^ is set equal to zero. Figure 95 shows 
the calculated values of open loop gain and phase margin. It 
is important to note at this point that equation (9-25) is 

an exact equation and that the assumptions made earlier in 

/ 

the feedforward design process are not made in deriving this 
equation. 

It is seen from Figure 95 (a) and (b) that the two 

stage input filter causes disturbances in the open loop gain 
euid phase at the two resonant frequencies, but the addition 
of feedforward removes these disturbances. The feedforward 
loop is seen to be effective at both the resonant frequen- 
cies though at the higher frequency there is a drop in phase 
margin. The two assumptions made in the feedforward design 
process have not been used in calculating the open loop gain 
with feedforward, and thus it is seen that Figure 95 indi- 
cates that the feedforward is effective in eliminating input 
filter interaction with the regulator control loop. 

The same regulator circuit was used to obtain experi- 
mental results using a single stage input filter: 

R_. =0.2 ohm L =116 microH Cl = 20 microF 

u 1 

The regulator was set up to measure the open loop gain and 
phase margin [8]. The nonadaptive feedforward circuit of 
Figure 67 was used in making the measurements. The feedfor- 


FREQ (HZ) 


Figure 95: (a) Open loop gain: Calculated values without 

input filter, with (A) and without feedforward 




O I 
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1 .E+4 
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Figure 95: (b) Open loop phase margin: Calculated values 

without input filter, and (A) and without 
feedforward 
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ward gain was calculated from equation (9-22) as -0.009 and 
this led to the following choice of feedforward circuit par- 
ameters : 

C- = 27 microF R., =5.6 Kohm R., = 50 ohm 
r f 1 I Z 

The input of the feedforward circuit is the input filter ca- 
pacitor voltage, and the output is subtracted from the out- 
put at the integrater in the feedback loop, as was done ear- 
lier. 

Figure 96 shows some preliminary measurements obtained. 
It is noticed that the input filter causes a disturbance in 
open loop gain and phase margin at the resonant frequency 
but the feedforward compensates for the input filter inter- 
action to some extent. The approximation of equation (9-20) 
is valid at low frequencies, thus if the input filter reso- 
nates at high frequencies the feedforward as designed may 
not be effective in completely compensating for input filter 
interaction. Measurements of Figure 96 were obtained using 
an input filter with a somewhat high resonant frequency 
which can be used to explain the fact that the feedforward 
is effective to a somewhat lesser degree in conpensating for 
the input filter interaction. Due to lack of time extensive 
experimental verification of the feedforward design was not 
possible; more measurement data is necessary to check the 
effect of feedforward on open loop gain and phase, and other 
performance specifications. 
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Figure 96: (a) Open loop gain: Measured values without 

input filter, with (•) and without feedforward 
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Figure 96: (b) Open loop phase margin; Measured values 

without input filter, with ( •} and without 
feedforward 
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The feedforward design process for the boost regulator 
would be identical to the one presented above; similar as- 
sumptions may be necessary to accomodate the positive zero 
in the boost regulator transfer function. 

9.3 REMARKS 

The following observations are made: 

1. The feedforward design process presented in Chapter 5 
is general; it can be used for any type of duty cycle 
control and for single loop or multi -loop feedback 
control schemes. 

2. The problem of input filter interaction with the re- 
gulator control loop is present in the buck-boost re- 
gulator also and the concept of f-sedforward compensa- 
tion can be extended to design a feedforward loop 
that eliminates the interaction. 



Chapter X 

CONCLUSIONS AND FUTURE WORK 
10.1 CONCLUSIONS 

A novel scheme for input filter compensation is presented in 
this dissertation- The input filter, while a necessary com- 
ponent of a switching regulator, interacts with the regulator 
control loop and this interaction can cause serious perfor- 
mance degradation such as loop instability, degradation of 
transient response, audiosusceptibility and the output impe- 
dance characteristics (3,4,5,6,10]. These problems are 
caused mainly by the peaking of the output impedance of the 
input filter and its interaction with the control loop. 

Conventional input filter design techniques for single 
stage and two stage input filters [3,5,6] minimize the peak- 
ing effect but this often results in a penalty of weight or 
loss increase in the input filter. A novel feedforward com- 
pensation scheme is developed for buck regulators that eli- 
minates the undesirable interaction between the input filter 
and the regulator control loop. The compensation scheme is 
implemented by a feedforward loop that senses the input fil- 
ter state variables. Extensive analytical and experimental 
data confirm that the regulator is made immune to the peak- 
ing of the input filter output impedance. 
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The following points regarding the feedforward 
compensation scheme presented are noteworthy: 

1. The feedforward design process is straightforward and 
it is seen that the feedforward loop gain is indepen- 
dent of the type of feedback control used and of the 
feedback loop transfer functions. 

2. An important characteristic of the feedforward con- 
trol scheme presented is its ease in implementation. 

3 . Extensive experimental data supported by analytical 
results show that significant improvement in perfor- 
meuice is.achieved with the use of feedforward in the 
following categories-- 

(a) open loop gain and phase margin; 

(b) audiosusceptibility; 

(c) output impedance; and 

(d) transient response. 

The data shows clearly that the peaking effect of the 
input filter output impedance- is eliminated with the 
use of feedforward. 

4. It is shown that the feedforward loop is independent 
of the input filter parameters and also. independent 
of the input filter configuration. 

Examination of the extensive experimental data con- 
firm that the use of the feedforward compensation 
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scheme does not have any detrimental effect on regu- 
lator performance. 

6. If a switching regulator is acquired as a 'black box' 
for use in a system, providing an inadequate input 
filter could lead to system instability [4,5]. This 
is shown analytically using a program that calculates 
the eigenvalues of the system, eind confirmed experi- 
mentally. It is then shown that the use of the feed- 
forward compensation scheme can stabilize a system 
that was made unstable due to input filter interac- 
tion with the regulator control loop. Experimental 
verification of the stabilizing action of the feed- 
forward loop is presented and this shows that an ar- 
bitrary input filter may be used with the feedfor- 
ward- loop controlled regulator. 

7. From the results presented it is logical to conclude 
that the feedforward loop can provide effective com- 
pensation for an unknown source impedance. For exam- 
ple, a preregulator which often has an unknown, dy- 
namic output impedance can interact with a switching 
regulator downstream and result in system . instabili- 
ty. The feedforward circuit effectively isolates the 
switching regulator from its source thus preventing 
interaction between the switching converter and 
equipment upstream. 
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The use of the feedforward compensation scheme re- 
moves some of the input filter design constraints 
thus making the input filter design process simpler 
and allows the input filter to be optimized. A high 
performance regulator system with an optimum input 
filter can thus be realized with the use of the feed- 
forward compensation scheme. 

9. The feedforward design concept presented can be easi- 
ly extended to types of control other than those used 
in this dissertation. 

10. Extension of. the concept of feedforward compensation 
to other types of switching regulators is possible. 
Such a scheme for a buck-boost regulator is presented 
and preliminary results indicate the validity of the 
compensation scheme. 

10.2 SUGGESTIONS FOR FUTURE WORK 

The following suggestions for future work are made: 

The concept of feedforward compensation for buck- 
boost switching regulators needs to be verified ex- 
perimentally. Data indicating the effect of using 
the feedforward scheme on performance categories like 
audiosusceptibility, output impedance, transient res- 
ponse and stability is also needed. 


1 . 
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2 . An investigation into the use of the feedforward com- 
pensation scheme for the boost converter should be 
made . 

3. The use of the feedforward loop is shown to isolate 
the regulator from its source, and this could be used 
to advantage in situations where regulators are put 
in series/parallel for load sharing. 

4. Current injected control is an important mode of con- 
trol and an investigation of the feasibility of using 
feedforward compensation for this type of control 
could be made. 

5. Single loop control using the constant frequency duty 
cycle type control is very widely used and test re- 
sults using feedforward compensation for this type of 
control are needed. 

6. It may be possible to simplify the adaptive feedfor- 
ward circuit by using the dc component of the input 
filter capacitor voltage instead of the supply vol- 
tage. This may have an effect on the transient res- 
ponse for a step change in supply voltage since then 
the duty cycle would change in accordance with chang- 
es in input filter capacitor voltage. 



Appendix A 


SOFTWARE FOR 

DISCRETE TIME SIMULATION OF A BUCK REGULATOR 
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***###*******-H*****-IHf*-H-*#****#***»*lf*1HHt*»******#***-tHf**'H-** 

* DISCRETE TIME SIMULATION OF A BUCK CONVERTER. * 

•tHHt1f*****#******************#**#*-M-*#****#»***H**#*****Ht*')t** 


THIS PROGRAM SIMULATES THE STEADY-STATE BEHAVIOR AS WELL AS 
THE TRANSIENT PERFORMANCE FOR A BUCK REGULATOR 
WITH INPUT FILTER, WITHOUT FEEDFORWARD. 


C 

C 


C 


DIMENSION X(6,9000). F1(6,6),F2(6,6).Gl(6,4).G2(6,t|),TIME(9000), 
1U(»4), PHI F(6,6), PHIN(6,6).DF(6.4),DN(6.4).TVEC1(6),TSWIT{2), 
2VISWIT(2),VO(9000), PHI FT( 6,6), DFT( 6,4), TVEC2( 9000) 

DIMENSION F3(6,6),G3(6,4),TVEC3(6),PHIFT3(6,6),DFT3(6,4), 
1PHIF3(6,6),DF3(6,4) 

REAL L1,L,K1,K2, IQMAX 

DATA RL1,L1,C1,L,RO,RN,RC,C/1,0,325.E-6,20.E-5,230.E-6,0.2,0.65* 
20.067, 300. E-6/ 

DATA RL,R11,R12,R13,R14,C2,R4,CP1/20.0,33.3E3, 16.7E3,2.E3, 
147.E3,0.01E-6.40.7E3,5600.E-12/ 

DATA ER,ET,VI ,EO/6.7,7.0,40.,20./ 

DATA TON,TOFMIN,TSWIT,VISWIT,VITON/2.20E-5,5.E-6,0.003, 10. , 
125.0,50, .0.88E-3/ 

DATA NIT, FRAC,EPS,TF/10,0.5,5.E-6,0.012/ 

DATA IQMAX,EQ,E0/8.0,0.2,0.7/ 

CL1S=0. 5291665 
VC1S=39. 48866 
ECS=7,0 

CLS=0. 09291243 
VCS=20.05715 
ERS=20. 05712 
VOS=19. 99638 

TON=VITON/VI 

K1=(R12/(R14*(R11+R12) ) )+( 1 . /R1 3 )-( RN/R4 ) 

K1=K1-( RN*RO*(RL+RC) )/(R4*RL*RC) 

K2=(RN/R4)-( 1 ./R13)-(R12/(R14*(R11+R12) ) ) 

DO 1 1=1,6 
DO 80 J=1,4 
G1 ( I , J )=0. 

G2(I,J)=0. 

DF( l,J)=0. 

DN( I , J )=0. 
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DFT( l,J)=0. 

DO 1 J=1,6 
FU l,J)=0. 

1 F2(I,J)=0. 

C 

DO 2 1=1,6 
DO 2 J=1,6 
PHIFT( l,J)=0. 

PHIF( 1,J)=0. 

2 PH1N(1,J)=0. 

C 

F1(1.1)=(-1.*RL1)/L1 

F1( 1,2) = -1 ./LI 

F1(2,l)=1./C1 

F1(2.4)=-1./C1 

F1( 3,2)=(-1 .*RN)/(CP1*R4) 

F1(3,4)=(-1 .*RL*RC*K1 )/ ( CPI *( RL+RC ) ) 

F1(3,5)=(RL*K2)/(CP1*(RL+RC)) 

F1(3,6)=1./(CP1*R13) 

Fl(4,2)=1,/L 

Fl(4,4)=(-1 .*R0)/L+(-1.*RC»RL)/(L*(RC+RL) ) 
F1(4,5)=(-1.*RL)/(L*(RC+RL) ) 
F1(5,4)=RL/(C*(RL+RC)) 

F1(5,5) = -1./(C<^(RL+RC)) 
F1(6,4)=(RL*RC)/(C2»R13*(RL+RC) ) 
F1(6,5)=RL/(C2*R13*(RL+RC)) 
F1(6,6)=-1./(C2»R13) 

Gl(l,l)=1./L1 

G1(3,2)=1./(CP1*R14) 

G1(3,3)=RN/(CP1*R4) 

G1(4,3)=-1./L 
DO 3 1=1,6 
DO 3 J=1,6 

3 F2( 1,J)=F1( l,J) 

F2(2,4)=0. 

F2(3,2)=0. 

F2(4,2)=0. 

G2(1,1)=G1(1,1) 

,G2( 3,2)=G1 (3,2) 

G2(3,4)=RN/(CP1*R4) 

G2{4,4)=-1 ./L 
DO 530 1=1,6 
DO 531J=1,6 



531 F3( l,J)=F2( I.J) 

DO 532 J=1,4 

532 G3( I ,J )=G2( I ,J ) 

530 CONTINUE 

DO 535 1=1,6 
F3(4,l)=0.0 

535 F3(l,4)=0.0 
DO 536 1=1,4 

536 G3(4,l)=0.0 
C 

NSWIT=0 
ITEMP=NSWIT+1 
TXX=TSWIT( I TEMP) 

NT=0. 1+1 ./frag 

TOFF=TON*(VI-EO)/EO 

TIME(1)=0.0 

ND=1 

M=0 

IT=0 

T=FRAC»TOFF 

C 

CALL STRAN(T,PHIF,DF,F2,G2,6,4,6,4,1,0,EPS.DIFMAX, ITER, IFLAG) 
C 

T=FRAC*TON 

CALL STRAN(T,PHIN,DN,F1,G1,6,4,6,4,1,0,EPS,D1FMAX, ITER, IFLAG) 
WRITE(6,704) 

704 FORMAT(//,' VALUES OF X ON THE NEXT PAGE') 

C 

X( 1 ,1 )=CL1S 
X(2, 1)=VC1S 
X(3,1)=ECS 
X(4,1)=CLS 
X(5,1)=VCS 
X(6,1)=ERS 
V0(1)=V0S 
U( 1 )=VI 
U(2)=ER 
U(3)=EQ 
U(4)=ED 
C 
C 

I F(X( 3,ND) .GE. ET) GO TO 60 
TEMPA=ABS(X(3,ND)-ET) 


o o o o o 


1F( TEMPA.lt. EPS) GO TO 40 
20 CONTINUE 

IF(X(3,ND).LT.ET) GO TO 28 

*************************************** 

MINIMUM OFF TIME CALCULATION * 

*************************************** 


T=TOFMIN 

CALL STRAN(T.PHIFT,DFT,F2,G2,6,4,6,4,1,0,EPS,DIFMAX, ITER, IFLAG) 

DO 6 1=1,6 

6 TVEC1( I )=X( l,ND) 

CALL FUTURE( TVEC1 ,U, PH I FT, OFT, 6, 6, 6, 4, 6, 6, 6, 4, 1 ) 

TEMP=TIME{ND) 

ND=ND+1 

TIME(ND)=TEMP+T 
DO 7 1=1,6 

7 X( I ,ND) = TVEC1( I ) 

VO( NO )=( RL*RC )*X( 4, ND )/( RL+RC)+RL*X( 5, ND)/( RL+RC ) 
IF(X(4,ND).GT.0.0) GO TO 570 
IT 1=0 

571 IF(ABS(X(4,N0)).LT.EPS) GO TO 579 
ITI=ITI+1 

IF(NIT.GE.NIT) GO TO 100 

SL0PE3=F2(4,4)*X(4,ND)+F2(4,5)*X(5,ND)+G2(4,4)*U(4) 

T3=-X(4,NO)/SLOPE3 

CALL STRAN(T3, PHI FT3,DFT3,F2,G2,6,4,6,4, 1,0,EPS,DI FMAX, ITER, I FLAG) 
DO 572 1=1,6 

572 TVEC3( I )=X( I ,N0) 

CALL FUTURE! TVEC3,U, PH I FT3 , DFT3, 6, 6, 6, 4, 6, 6, 6, 4, 1 ) 

TIME(ND)=TIME(ND)+T3 

DO 573 1=1,6 

573 X( I ,ND)=TVEC3( I ) 

V0( ND )=RL*RC*X( 4 , ND ) +RL*X( 5 , ND ) 

VO(ND)=VO(ND)/( RL+RC) 

GO TO 571 
579 CONTINUE 

CALL STRAN(T3,PHIF3,DF3,F3,G3,6,4,6,4,1,0,EPS,DIFMAX, ITER, IFLAG) 

DO 574 1=1,6 

574 TVEC3( I )=X( l,ND) 

X(4,N0)=0.0 

CALL FUTURE(TVEC3,U,PHIF3,DF3,6,6,6,4,6,6,6,4,1 ) 

TEMP=TIME(ND) 
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ND=ND+1 

TIME(N0)=TEMP+T3 
DO 575 1=1,6 
575 X( I ,ND)=TVEC3( I ) 

VO(ND)=(RL*X(5,ND) )/(RL+RC) 

570 CONTINUE 
GO TO 40 
C 

Q 

C OFF TIME CALCULATION * 

C 

28 T=FRAC*TOFF 
IT=0 

25 CONTINUE 

DO 8 1=1,6 

8 TVEC1( I )=X( I ,ND) 

IF(N0.LT.300) GO TO 5000 
WRITE(6,706) 

706 FORMAT(//,' STATE VECTOR IN OFF TIME') 
WRITE(6,9004)(TVEC1( I ), 1=1,6) 

9004 FORMAT(6F15. 1 1 ) 

WRITE(6,*) VO(ND) 

5000 CONTINUE 

IF(X(4.ND).LE.O) GO TO 502 
C 

CALL FUTURE(TVEC1,U,PHIF,DF,6,6,6,4,6,6,6,4,1 ) 
C 

TEMP=TIME(ND) 

ND=ND+1 

TIME(NO)=TEMP+T 
DO 9 1=1,6 

9 X( l,ND) = TVEC1( I ) 

VO( ND)=RL*RC*X( 4, ND)+RL*X( 5, ND) 
VO(ND)=VO(NO)/(RL+RC) 

GO TO 520 
502 CONTINUE 
C 

C ************************************* 

C * CALCULATION TO FIND WHEN X(4)=0 » 

C #***#*#***#*•*■#•«■*****#**********##*»# 

c 

IT 1=0 
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503 IF(ABS(X(4.ND)).LT.EPS) GO TO 509 
m=ITI+1 

IF( (Tl .GE.NIT) GO TO 100 

SLOPE3=F2(4.4)*X(4,ND)+F2(4,5)*X{5.ND)+G2{4,4)*U(4) 

T3=-X(4.ND)/SLOPE3 

CALL STRAN(T3, PHI FT3,DFT3, F2, G2, 6, 4, 6, 4, 1 . 0. EPS, D I FMAX, ITER, I FLAG) 
DO 504 1=1,6 

504 TVEC3( I )=X{ I ,ND) 

CALL FUTURE(TVEC3,U, PH I FT3 , DFT3, 6, 6, 6, 4, 6, 6, 6, 4, 1 ) 

TIME(ND)=TIME(ND)+T3 

00 505 1=1,6 

505 X( l,ND)=TVEC3( I ) 

VO(ND)=RL*RC*X(4,ND)+RL*X(5.ND) 

VO(ND)=VO(ND)/(RL+RC) 

GO TO 503 
509 CONTINUE 

CALL STRAN(T,PHIF3,DF3,F3,G3,6,4,6,4,1,0,EPS,DIFMAX, ITER, I FLAG) 

DO 5401=1,6 
540 TVEC3( I )=X( I ,ND) 

X(4,ND)=0.0 

CALL FUTURE(TVEC3,U,PH1F3,DF3,6,6,6,4,6,6,6,4,1 ) 

TEMP=TIME(ND) 

ND=ND+1 

TIME(ND)=TEMP+T 
DO 511 1=1,6 
511 X( I ,ND)=TVEC3( I ) 

V0(ND)=(RL*X(5,ND) )/(RL+RC) 

520 CONTINUE 

IF(X(3,ND).GE.ET) GO TO 30 
GO TO 25 

■If*********#**#*##*******#***#####****** 

CALCULATION TO HIT THRESHOLD WHEN * 

X(3)=ET. * 

*************************************** 


30 IF(ABS(X(3,ND)-ET).LT.EPS) GO TO 40 
IT=IT+1 

IF( IT. GE.NIT) GO TO 100 

SLOPE=F2(3,4)*X(4,ND)+F2(3,5)^>X(5,ND)+F2(3,6)*X(6,ND) 
SL0PE=SL0PE+G2{ 3,2)*U(2) 

T=(ET-X(3,ND))/SL0PE 



ooooooo 
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CALL STRAN(T,PHIFT,DFT.F2,G2,6.i|.6.4.1.0,EPS, DIFMAX. ITER, I FLAG) 
C 

DO 10 1=1,6 

10 TVEC1( I )=X( l,ND) 

CALL FUTURE(TVEC1,U,PHI FT,DFT,6,6.6,4.6,6.6.4, 1 ) 

TIME(ND)=TIME(ND)+T 

DO 11 1=1,6 

11 X{ l,ND)=TVEC1( I ) 

VO(ND)=RL*RC*X(4,ND)+RL*X(5,ND) 

VO(ND)=VO(NO)/(RL+RC) 

GO TO 30 
40 CONTINUE 

IF(ND.LT.300) GO TO 5001 
WRITE(6,440) TIME(ND) 

440 FORMAT(/, 'SWITCH ON T I ME= ' , F15. 1 1 ) 

5001 CONTINUE 

IF(TIME(ND).GE.TF) GO TO 120 

*«■»#•»*##»#»**#»»#»***#*****#»#«»******* 

CHANGE VI, TON, TOFF AT THE SWITCHING * 

INSTANT. ALSO CHANGE PHIF,DF, IF * 

NECESSARY. * 

*#****■»**#»«*****<»#*«*•»******»******#** 

IF(TIME(ND).LT.TXX) GO TO 60 
NSWIT=NSWIT+1 
ITEMP=NSWIT+1 
TXX=TSWIT( ITEMP) 

VI=VISWIT(NSWIT) 

TON=VITON/VI 
TOFF=TON*(VI-EO)/EO 
U(1)=VI 
t=frac*Toff 

C 

CALL STRAN(T,PHIF,DF,F2,G2,6,4,6,4,1,0,EPS, DIFMAX, ITER, I FLAG) 

C 

T=FRAC*TON 

C 

CALL STRAN(T,PHIN,DN,F1,G1,6,4,6,4,1,0,EPS,DIFMAX, ITER, IFLAG) 

C 

60 CONTINUE 
M=0 

T=FRAC*TON 

C 



nooooo o no, oooo 


ON TIME CALCULATION * 

65 CONTINUE 

DO 12 1=1,6 

12 TVEC1( I )=X( I.ND) 

I F(ND.LT.300) GO TO 5002 
WRITE(6.709) 

709 FORMAT(//,' STATE VECTOR IN ON TIME') 
WRITE(6,9005)(TVEC1( I ), 1=1,6) 

9005 FORMAT(6F15. 11 ) 

WRITE(6,*) VO(ND) 

5002 CONTINUE 

CALL FUTURE(TVEC1,U, PH I N, ON, 6, 6, 6, *1. 6, 6, 6, 4. 1 ) 

TEMP=TIME(ND) 

ND=ND+1 

TIME(ND)=TEMP+T 
DO 13 1=1,6 

13 X( I ,ND)=TVEC1( I ) 

V0( ND)=RL*RC»X( 4, ND)+RL*X( 5, NO) 
V0(ND)=V0(ND)/(RL+RC) 

M=M+1 

IF(X(4,ND).GE. IQMAX) GO TO 70 
IF(M.LT.NT) GO TO 65 

GO TO 999 

*»***#*#*»*#***#*#»**#**»***»**##****** 
CALCULATION TO FIND WHEN INDUCTOR * 

CURRENT HITS THE THRESHOLD VALUE. * 

*************************************** 


70 CONTINUE 
IT=0 

71 IF(ABS(X(4,ND)-IQMAX).LE.EPS) GO TO 999 
IT=IT+1 

IF( IT.GE.NIT) GO TO 100 

SL0PE2=X(2,ND)*F1(4,2)+X(4,ND)»F1(4,4)+X(5.ND)*F1(4.5) 

T=( IQMAX-X(4,ND) )/SLOPE2 

CALL STRAN(T,PHIFT,DFT,F1,G1,6,4,6,4,1,0,EPS,DIFMAX, ITER. I FLAG) 



no 


DO 15 1=1,6 

15 TVEC1( I )=X( I.ND) 

CALL FUTURE(TVEC1.U,PHIFT,DFT.6,6.6.«I,6.6.6.'<.1) 

TIME(ND)=TIME(ND)+T 

DO 16 1=1,6 

16 X( I ,ND)=TVEC1 ( I ) 

VO( ND)=RL*RC»X(4,ND)+RL*X(5,ND) 

VO{ NO )=VO( NO ) / ( RL+RC ) 

GO TO 71 
999 CONTINUE 
C 

IF(ND.LT.300) GO TO 5003 
WRITE{6,442) TIME(ND) 

442 F0RMAT(/, 'SWITCH OFF T I ME= ' , F15. 1 1 ) 

5003 CONTINUE 

1F(T1ME(ND).LT.TF) GO TO 20 
GO TO 120 
100 CONTINUE 

WRITE(6,204) TIME(ND) 

204 FORMAT(/, 'MAXIMUM ITERATION AT TIME=' , F15. 1 1) 

120 CONTINUE 

DO 104 1=1, NO 

WRITE(8,103) T1ME( 1 ),X(3, 1 ) 

103 FORMAT(2E20. 10) 

104 CONTINUE 

DO 105 1=1, ND 

WRITE(8, 103) T1ME( I ),X(4, I ) 

105 CONTINUE 

DO 106 1=1, ND 

WRITE(8,103) TIME( I ),VO( I ) 

106 CONTINUE 

DO 107 1=1, NO 

WRITE(8,103) TIME( I ),X(2, 1 ) 

107 CONTINUE 

DO 14 1=1, ND 
14 TVEC2( I )=VO( I ) 

CALL ABPLOT(NO, TIME, TVEC2) 

STOP 

END 


1 SUBROUTINE STRAN( TAU, PH I , THETA, A, B, NA, NB, MA, MB, MODE, NTERMS, 
1 TOL,DI FMAX, ITER, I FLAG) 



REAL PH I ( NA, NA ) . THETA( NA, NB ) , A( NA, NA ) , B( NA, NB ) , 

1 WORK! ( 10, 10),WORK2( 10, 10 ) , DUMMY( 1 , 1 ) , FAC1 , CONI , FAC2, C0N2 
C1 = l 
TS=TAU 
DO 4 1=1, MA 
DO 2 J=1,MA 
WORK1( l,J)=0 
W0RK2( I ,J )=0 
PHI{ l,J)=0 
2 CONTINUE 
WORK1 ( I , I )=1 
WORK2( I , I )=TAU 
PHI( I, I ) = 1 
4 CONTINUE 
DO 6 l=1,MA 
DO 6 J=1,MB 
6 THETA( I ,J)=0 
DIFMAX=1.E8 
ND0=50 

I F( NTERMS.GT.O) NDO=NTERMS 
FAC1=1 
I FLAG=0 

DO 1000 1=1, NDO 

CALL MXMUL(0UMMY,W0RK1.A,TS, 1, 1, 10, 10, NA, NA, MA, MA, MA, MA, 1 ) 

FAC1=FAC1*I 

C0N1=1 ./FAC1 

IF(NTERMS.EQ.O.AND. I,GE.4) CALL SERROR( PH I , W0RK1 , CONI , NA, NA, 

1 10, 10,MA,MA,DIFMAX) 

CALL MXADD( DUMMY, PH I , WORK! , Cl , CONI , 1 , 1 , NA, NA, 10, 10,MA,MA, 1 ) 

I F(MODE. EQ.2) GO TO 500 
FAC2=FAC1*( 1+1 ) 

C0N2=TAU/FAC2 

CALL MXA0D(DUMMY,W0RK2,W0RK1 ,C1,CON2, 1, 1, 10, 1 0, 1 0, 10, MA, MA, 1 ) 
500 CONTINUE 
ITER=I 

IF(NTERMS.GT.0.0R. I ,LT.4) GO TO 1000 
IF(DIFMAX.LE.TOL) GO TO 1100 
1000 CONTINUE 
1100 CONTINUE 

IF(MODE.EQ.I) CALL MXMUL( THETA, W0RK2. B, Cl , NA, NB, 10, 1 0, NA, NB, 

1 MA,MA,MA,MB,2) 

IF( ITER. EQ. NDO. AND. NTERMS.EQ.O) IFLAG=1 
RETURN 
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END 


1 SUBROUTINE SERROR( AMX, BMX, CCC, I A. JA, I B, JB. I DO, JDO, D I FMAX ) 
' DIMENSION AMX( IA,JA),BMX( IB, JB) 

Dl FMAX=1 . E-30 
DO 100 1 = 1, IDO 
DO 50 J=1,JD0 

I F(AMX( I , J ) . EQ.O.O ) GO TO 50 
CHANGE= ABS(BMX( I , J ) *CCG/AMX( l,J)) 

I F ( CHANGE. GT.D I FMAX) D I FMAX=CHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 


1 SUBROUTINE FUTURE(X,V, PH I , THETA, LP, MP, LT, MT, I P.JP, IT,JT,MODE) 
DIMENSION X(MP),V(MT), PHI(LP,MP),THETA( LT, MT ) , TEMP( 20 ) 

DO 20 1 = 1, IP 
SUM=0 

DO 10 J=1,JP 
10 SUM=SUM+PHI( l,J)*X(J) 

20 TEMP( I )=SUM 
DO 30 1 = 1 , I P 
30 X( I ) = TEMP( I ) ' 

I F(MODE. EQ.2) RETURN 
DO 60 1 = 1, IT 
SUM=0 

DO 50 J=1,JT 

50 SUM=SUM+THETA( I , J )*V( J ) 

X( I )=X( I )+SUM 
60 CONTINUE 
RETURN 
END 


1 SUBROUTINE MXADD( RMX, AMX, BMX, ACC, BCC, IR,JR, IA,JA, IB,JB, IDO, JDO, 
1 MODE) 

DIMENSION AMX( IA,JA),BMX{ IB,JB),RMX( IR,JR) 

IF( lA.LT. IDO. OR.JA. lt. JDO) GO TO 999 
IF( IB. LT. IDO. OR.JB.lt. JDO) GO TO 999 
GO TO ( 10, 100), MODE 
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10 CONTINUE 

DO 50 1=1. IDO 
DO 50 J=1.JD0 

50 AMX( l,J)=AMX( I , J )*ACC+BMX( l,J)*BCC 
GO TO 300 
100 CONTINUE 

IF( IR.LT. IDO. OR.JR.lt. JDO) GO TO 999 
DO 200 1=1, IDO 
DO 200 J=1,JD0 

200 RMX( l,J)=AMX( I , J ) *ACC+BMX( l,J)*BCC 
300 RETURN 
999 CONTINUE 
RETURN 
END 


1 SUBROUTINE MXMUL( RMX, AMX, BMX.CCC, IR.JR, lA.JA, IB.JB, IDA.JDA, IDB.JDB 
1.M0DE) 

DIMENSION AMX( IA,JA),BMX( IB,JB),RMX( I R, JR ) , TEMP( 20 ) 

I F( IDA*IDB*JDA*JOB.GT. IA»IB*JA*JB) GO TO 999 
IF(JOB.GT.JDA) GO TO 999 
GO TO ( 10, 210), MODE 
10 DO 100 1=1, IDA 
DO 20 L=1,JDA 
20 TEMP(L)=AMX( l,L) 

DO 80 J=1,JDB 
SUM=0 

DO 40 K=1,JDA 

40 SUM=SUM+TEMP( K)*BMX( K,J ) 

AMX( I , J )=SUM*CCC 
80 CONTINUE 
100 CONTINUE 
GO TO 600 
210 CONTINUE 

IF( IR.LT. IDA. OR.JR.lt, JOB) GO TO 999 
DO 400 1=1, IDA 
DO 380 J=1,JDB 
SUM=0 

DO 340 K=1,JDA 

340 SUM=SUM+AMX( I , K)*BMX( K, J ) 

RMX( I ,J ) = SUM*CCC 
380 CONTINUE 
400 CONTINUE 



600 RETURN 
999 CONTINUE 
RETURN 
END 
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■H4Ht**#*«*tt**-H******************»*********#******lHf»H4Hf-lt#-IHt* 

DISCRETE TIME SIMULATION OF A BUCK CONVERTER * 

THIS PROGRAM SIMULATES THE STEADY-STATE BEHAVIOR AS WELL AS 
THE TRANSIENT PERFORMANCE FOR A BUCK REGULATOR 
WITH INPUT FILTER AND WITH FEEDFORWARD. 


C 

C 


C 


DIMENSION X(7,9000). F1{7,7), F2(7,7),G1(7,4),G2(7,l4),TIME(9000), 
1U(4). PHIF(7.7). PHIN(7.7),DF(7.4).DN(7.4).TVEC1(7),TSWIT(2). 

2VI SWI T( 2 ) , VO( 9000 ) . PH I FT( 7, 7 ) . DFT( 7. 4 ) , TVEC2( 9000 ) 

DIMENSION F3(7,7),G3(7,4),TVEC3(7).PHIFT3(7,7),DFT3(7,4), 
1PHIF3(7,7),DF3(7,7) 

REAL L1,L,K1,K2, IQMAX.K3 

DATA RL1 , L1.C1 , L, RO, RN, RC, C/1 . 0. 325 . E-6. 22. E-5, 230 . E-6, 0 . 2. 0 . 65 . 
20.067.300.E-6/ 

DATA RL, R11,R12,R13,R14.C2,R4. CPI/20.0, 33. 3E3, 16. 7E3, 2. E3, 
147.E3,0.01E-6,40.7E3,5600. E-12/ 

DATA ER.ET.VI ,EO/6.7,7.0,40. ,20./ 

DATA TON,T0FMlN,TSWIT,VISWIT,VITON/2.20E-5,5.E-6,0.003, 10.0, 

125. ,50. ,0.88E-3/ 

DATA NIT, FRAC,EPS,TF/10,0. 1, 10,E-6,0.0120/ 

DATA IQMAX,EQ,ED/6.0.0,2,0.7/ 

DATA RF1,RF2.CF,VFS/5. 1E3,90.970,27.E-6,0.0507E-1/ 

CL1S=0. 69552390 
VC1S=39. 87245 
ECS=7.0 

CLS=0. 40158130 
VCS=20. 06141 
ERS=20. 06143 
VOS=20. 02125 
VFS=+0. 38948170 

TON=VITON/VI 

K1=(R12/(R14*(R11+R12) ) )+( 1 . /R1 3 )-( RN/R4 ) 

K1=K1-(RN*R0*(RL+RC) )/(R4*RL»RC) 

K2=(RN/R4)-( 1 ./R13)-(R12/(R14*(R11+R12) ) ) 

K3=RF2/(RF1+RF2) 

DO 1 1=1,7 , 

DO 80 J=1,4 
G1(L,J)=0. 

G2( l,J)=0. 
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DF( I.J)=0. 

DN( l,J)=0. 

80 DFT(l,J)=0. 

DO 1 J=1,7 
F1( I ,J )=0. 

1 F2(I,J)=0. 

C 

DO 2 1=1,7 
DO 2 J=l,7 
PHI FT( I , J )=0. 

PHIF( l,J)=0. 

2 PHI N( I , J )=0. 

C 

Fl( 1, 1 )=(-1 .*RL1 )/L1 

Fl(1,2)=-1./L1 

F1(2,1)=1./C1 

F1(2,3)=-1./(C1»(RF1+RF2)) 

F1(2,5)=-1./C1 
F1(3.1 )=1./C1 

F1(3,3)=(-1.*(CF+C1 ) )/(CF*C1»(RF1+RF2)) 

F1(3.3)=0 

F1(3,5)=-1./C1 

F1(4, 1 )=(-1 .»K3)/C1 

F1(4,2)=(-1.*RN)/(CP1*R4) 

FI (4. 3 )=( K3*(CF+C1 ) ) /( CF*C1 *( RF1 +RF2 ) ) 
F1(4,5)=( K3/C1 )-( (RL*RC*K1 )/( CP1*( RL+RC) ) ) 
F1(4,6)=(RL*K2)/(CP1*(RL+RC) ) 
F1(4,7)=1./(CP1*R13) 

Fl(5,2)=1./L 

F1(5,5)=( (-1 .*R0)/L)+( (-1 .*RC»RL)/(L*(RC+RL)) ) 
F1(5,6)=(-1.*RL)/(L*(RC+RL) ) 

F1(6,5) = RL/(C*(RL+RC) ) 

F1(6,6)=-1./(C*(RL+RC)) 
F1(7,5)=(RL*RC)/(C2*R13*(RL+RC) ) 
F1(7,6)=RL/(C2*R13*(RL+RC) ) 
F1(7,7)=-1./(C2*R13) 

G1(1,1)=1./L1 

G1(4,2)=1./(CP1*R14) 

G1(4,3)=RN/(CP1*R4) 

G1(5,3)=-1./L 
DO 3 1=1,7 
DO 3 J=1,7 

3 F2( l,J)=F1( I.J) 



F2(2.5)=0. 

F2(3,5)=0. 

F2(3,3)=0 

F2(4,2)=0. 

F2(4.5) = (-1 .*RL*RC*K1 )/( CP1»( RL+RC) ) 

F2(5,2)=0. 

G2(1.1)=G1(1.1) 

G2(4.2)=G1(4.2) 

G2(4,4)=RN/(GP1*R4) 

G2(5.4)=-1./L 
DO 530 1=1,7 
DO 531 J=1,7 

531 F3( 1,J)=F2( I.J) 

DO 532 J=1,4 

532 G3( l,J)=G2( I.J) 

530 CONTINUE 

DO 535 1=1,7 
F3(5.l)=0.0 

535 F3(l,5)=0.0 
DO 536 1=1,4 

536 G3(5.l)=0.0 
C 

NSWIT=0 

TEMP=NSWIT+1 

TXX=TSWIT(TEMP) 

NT=0. 1+1 ./FRAG 

TOFF=TON*(VI-EO)/EO 

TIME(1)=0.0 

ND=1 

M=0 

IT=0 

T=FRAC*TOFF 

C 

CALL STRAN(T,PHIF,DF,F2,G2,7,4,7,4,1,0,EPS,DIFMAX, ITER, IFLAG) 
C 

T=FRAC*TON 

CALL STRAN(T,PHIN,DN,F1,G1,7,4,7,4,1,0,EPS,DIFMAX, ITER, IFLAG) 
WRITE(6,704) 

704 FORMAT(//,' VALUES OF X ON THE NEXT PAGE') 

C 


X( 1, 1 )=CL1S 
X(2, 1 )=VC1S 
X(3,1 )=VFS 
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X(4. 1 )=ECS 
X(5.1 )=CLS 
X(6, 1 )=VCS 
X(7, 1 )=ERS 
V0( 1 )=VOS 
U( 1 )=VI 
U(2)=ER 
U(3)=EQ 
U(4)=E0 


IF(X(4,ND).GE.ET) GO TO 60 
TEMPA=ABS(X(4,ND)-ET) 

IF(TEMPA.LT.EPS) GO TO 40 
20 CONTINUE 

IF(X(4,n'd).LT.ET) go to 28 

»#**»*»#* 

MINIMUM OFF TIME CALCULATION * 

****•»***»#*****■«•*#***»*#****»##**#***#* 


T=TOFMIN 

CALL STRAN(T,PHIFT,0FT,F2,G2,7,4.7,4,1,0,EPS,DIFMAX, ITER, IFLAG) 

DO 6 1=1,7 

6 TVEC1 ( I )=X( I ,ND) 

CALL FUTURE(TVEC1 ,U, PH I FT, DFT, 7, 7, 7, 4, 7, 7, 7, 4, 1 ) 

TEMP=TIME(ND) 

ND=ND+1 

TIME(ND)=TEMP+T 
DO 7 1=1,7 

7 X( l,ND) = TVEC1( I ) 

VO(ND)=(RL*RC)*X(5,ND)/(RL+RC)+RL*X(6,ND)/(RL+RC) 

IF(X(5,ND).GT.0.0) GO TO 570 

ITI=0 

571 IF(ABS(X(5,ND)).LT.EPS) GO TO 579 
ITI=ITI+1 

IF( ITI .GE.NIT) GO TO 579 

SL0PE3=F2(5,5)*X(5,ND)+F2(5,6)*X(6,ND)+G2(5,4)*U(4) 

T3=-X(5,ND)/SLOPE3 

CALL STRAN( T3, PH I FT3 , DFT3 , F2, G2, 7,4, 7,4.1, 0, EPS, D I FMAX, I TER, I FLAG ) 
DO 572 1=1,7 

572 TVEC3( I )=X( I ,ND) 
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CALL FUTURE(TVEC3.U.PHIFT3,0FT3,7,7.7,4,7,7.7.4,1 ) 

TIME(ND)=TIME(ND)+T3 

DO 573 1=1,7 

573 X( l,ND)=TVEC3( I ) 

VO(ND)=RL«RC*X(5,ND)+RL*X(6,ND) 

VO(ND)=VO(ND)/(RL+RC) 

GO TO 571 
579 CONTINUE 

CALL STRAN(T3, PHI F3.DF3. F3, G3, 7, 4, 7, 4. 1 , 0, EPS, D I FMAX, ITER, I FLAG) 
DO 574 1=1,7 

574 TVEC3( I )=X( l,ND) 

X(5,ND)=0.0 

CALL FUTURE(TVEC3,U, PH I F3, DF3, 7. 7, 7.4, 7, 7, 7, 4, 1 ) 

TEMP=TIME( NO) 

N0=ND+1 

TIME(ND)=TEMP+T3 
DO 575 1=1,7 

575 X( I ,ND)=TVEC3( I ) 

V0(N0)=(RL*X(6,ND) )/(RL+RC) 

570 CONTINUE 
GO TO 40 

****#»*#»*•»*#*#*#*#*#*#**#**#*»*#*«***# 

OFF TIME CALCULATION * 

»«***********■»*******«**********#*#*#*# 

28 T=FRAC*TOFF 
IT=0 

25 CONTINUE 

DO 8 1=1,7 

8 TVEC1{ I )=X( l,ND) 

IF(ND.LT.300) GO TO 5000 
WRITE(6, 706) 

706 FORMAT(//,' STATE VECTOR IN OFF TIME') 

WRITE(6,9004)(TVEC1( I ), 1=1,7) 

9004 FORMAT(7F15. 11 ) 

WRITE(6,*) VO(ND) 

5000 CONTINUE 

I F(X( 5, ND) . LE.O.O) GO TO 502 
C 

CALL FUTURE(TVEC1,U, PH I F, DF, 7, 7. 7, 4, 7, 7. 7, 4, 1 ) 

C 

TEMP=TIME(ND) 
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ND=ND+1 

TIME(ND)=TEMP+T 
DO 9 1=1,7 

9 X( I.ND)=TVEC1( I ) 

VO( ND ) =RL*RC*X( 5 , ND ) +RL*X( 6. ND ) 
VO(ND)=VO(ND)/(RL+RC) 

GO TO 520 
502 CONTINUE 

************************************ 

* CALCULATION TO FIND WHEN X(5)=0 * 

************************************ 


ITI=0 

503 IF(ABS(X(5, NOn.LT.EPS) GO TO 509 
ITI=ITI+1 

IF( ITI .GE.NIT) GO TO 100 

SLOPE3=F2(5.5)*X(5,ND)+F2(5,6)*X(6,ND)+G2(5,4)*U(4) 

T3=-X(5,ND)/SLOPE3 

CALL STRAN(T3, PHI FT3,DFT3,F2,G2,7,4,7.4,1,0,EPS,DI FMAX, ITER, I FLAG) 
DO 504 1=1,7 

504 TVEC3( I )=X( I ,ND) 

CALL FUTURE(TVEC3,U,PHIFT3,DFT3,7,7,7,4,7,7,7,4,1 ) 

TIME(ND)=TIME(ND)+T3 

DO 505 1=1,7 

505 X( I ,ND) = TVEC3( I ) 

V0(ND)=RL»RC*X(5,ND)+RL*X(6,ND) 

VO(ND)=VO(ND)/(RL+RC) 

• GO TO 503 
509 CONTINUE 

CALL STRAN(T,PHIF3,DF3,F3,G3,7,4,7,4,1,0,EPS,DIFMAX, ITER, IFLAG) 

DO 540 1=1,7 
540 TVEC3( l)=X( I ,ND) 

X(5,ND)=0.0 

CALL FUTURE(TVEC3,U,PHIF3,DF3,7,7,7,4,7,7,7,4,1 ) 

TEMP=TIME(ND) 

ND=ND+1 

TIME( ND)=TEMP+T 
DO 511 1=1,7 

511 X( I , ND)=TVEC3( I ) ' 

V0(ND)=(RL*X(6,ND) )/{RL+RC) 

520 CONTINUE 

IF(X(4,ND).GE.ET) GO TO 30 
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Q *************************************** 

C CALCULATION TO HIT THRESHOLD WHEN * 

C X(3)=ET. * 

C *************************************** 

c 

30 IF(ABS(X(4,ND)-ET).LT.EPS) GO TO 40 

IT=IT+1 

I F( IT.GE.NIT) GO TO 100 

SL0PE=F2(4, 1 )*X( 1,ND)+F2(4,3)*X(3.ND)+F2(4,5)*X(5,ND)+ 
1F2(4,6)*X(6.ND)+F2(4,7)*X(7,N0) 
SL0PE=SL0PE+G2(4,2)*U(2)+G2(4,4)*U(4) 

T=(ET-X(4,ND) )/SL0PE 
C 

CALL STRAN(T,PHIFT,DFT,F2,G2,7.4,7,4,1,0,EPS,DIFMAX, ITER, I FLAG) 
C 

DO 10 1=1,7 

10 TVEC1( I )=X( I.NO) 

CALL FUTURE(TVEC1,U, PH I FT, DFT, 7. 7, 7, 4, 7, 7, 7, 4, 1 ) 

TIME(ND)=TIME(ND)+T 

DO 1 1 1=1,7 

11 X( l,ND)=TVEC1( I ) 

VO(ND)=RL*RC*X(5,ND)+RL*X(6,ND) 

VO(ND)=VO(ND)/(RL+RC) 

GO TO 30 
40 CONTINUE 

IF(ND.LT.IOOO) GO TO 5001 
WRITE(6,440) TIME(ND) 

440 FORMAT(/, 'SWITCH ON TIME=' , F15. 11 ) 

5001 CONTINUE 

IF(TIME(ND).GE.TF) GO TO 120 

****#*****«■*******#*#*#*#****»*#******* 

CHANGE VI,TON,TOFF AT THE SWITCHING * 

INSTANT. ALSO CHANGE PHIF,DF IF * 

NECESSARY. * 

*************************************** 


IF(TIME(ND).LT.TXX) GO TO 60 

NSWIT=NSWIT+1 

TEMP=NSWIT+1 

TXX=TSWIT(TEMP) 
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VI=VISWIT(NSWIT) 

TON=VITON/VI 
TOFF=TON*{ VI-EO)/EO 
U(1)=VI 
RF2=230.27 
K3=RF2/(RF1+RF2) 

F1(2,3) = -1 ./(C1*(RF1+RF2) ) 

F1(3,3)=(-1 .*(CF+C1 ) )/(CF*C1*(RFl+RF2) ) 

F1(4,1)=(-1.*K3)/C1 

F1{4, 3 )={ K3*(CF+C1 ) ) /{ CF*C1 »( RF1 +RF2 ) ) 

F1(4,5)=(K3/C1 )-( (RL*RC*K1 )/(CP1*(RL+RC) ) ) 

F2(2,3)=F1(2,3) 

F2(3,3)=F1(3,3) 

F2(4.1)=F1(4.1) 

F2(4.3)=F1(4,3) 

T=FRAC*TOFF 

C 

CALL STRAN(T, PHI F, OF, F2, G2, 7, 4, 7, 4, 1 , 0, EPS, D I FMAX, ITER, I FLAG) 
C 

T= FRACTION 
C 

CALL STRAN(T,PHIN,DN,F1,G1,7,4,7,4,1,0,EPS,DIFMAX, ITER, I FLAG) 
C 

60 CONTINUE 
M=0 

T=FRAC*TON 

*********#»***#***#»*********#**#**#*** 

ON TIME CALCULATION * 

**»***********##******»##****#*#******* 


65 CONTINUE 

DO 12 1=1,7 
12 TVEC1( I )=X( I ,ND) 

IF(ND.LT.IOOO) GO TO 5002 
WRITE(6,709) 

709 format;//, ' STATE VECTOR IN ON TIME') 
write; 6, 9005 )(TVEC1( I ), 1=1,7) 

9005 format; 7F15. 1 1 ) 
write;6,*) vo;nd) 

5002 CONTINUE 
C 

CALL FUTURE;TVEC1,U,PHIN,DN,7,7,7,4,7,7,7,4,1 ) 



oooooo 




TEMP=TIME(ND) 

ND=ND+1 

TIME{ND)=TEMP+T 
DO 13 1=1,7 
13 x( (,NO)=rvEci( n 

V0(ND)=RL*RC*X(5.N0)+RL*X(6.Nb) 

VO(ND)=VO(ND)/(RL+RC) 

M=M+1 

IF(X(5,ND).GE. IQMAX) GO TO 70 
IF(M.LT.NT) GO TO 65 
C 

GO TO 999 

CALCULATION TO FIND WHEN INDUCTOR * 
CURRENT HITS THE THRESHOLD VALUE. * 


70 CONTINUE 
IT=0 

71 IF(ABS(X|5,N0)-IQMAX).LE.EPS) GO TO 999 
IT=IT+1 

IF( IT.GE.NIT) GO TO 100 

SLOPE2=F1(5,2)*X(2,NO) + F1(5,5)^*X{5,ND)+F1(5,6)*X(6,ND) + 
1G1(5,3)*U(3) 

T=( IQMAX-X(5,ND) 1/SLOPE2 

CALL STRAN(T,PHIFT,DFT,F1,G1,7,4,7,4,1,0,EPS,DIFMAX, ITER, IFLAG) 
DO 15 1=1,7 

15 TVEC1( I )=X( I ,ND) 

CALL FUTURE(TVEC1,U,PHIFT,0FT,7,7,7,4,7,7,7,4, 1 ) 

TIME(ND)=TIME(ND)+T 

DO 16 1=1,7 

16 X( l,ND)=TVEC1( I ) 

VO( ND )=RL*RC*X( 5 , ND ) +RL*X( 6, NO ) 

VO(ND)=VO(ND)/(RL+RC) 

GO TO 71 
999 CONTINUE 
C 

IF(ND.LT. 1000) GO TO 5003 
WRITE(6,442) TIME(ND) 

442 FORMAT(/, 'SWITCH OFF TIME=' , F15. 1 1 ) 

5003 CONTINUE 

IF(TIME(ND).LT.TF) GO TO 20 
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GO TO 120 
100 CONTINUE 

WRITE(6,204) TIME(ND) 

204 FORMAT(/, 'MAXIMUM I TERAT I ON AT T I ME= ' , F15 . 1 1 ) 
120 CONTINUE 

111 DO 105 1=1, NO 

WRITE(8, 103) TIME( I ),X( 1. I ) 

105 CONTINUE 

DO 106 1=1, ND 

WRITE(8, 103) TIME( I ),X(2, I ) 

106 CONTINUE 

112 DO 104 1=1, ND 
WRITE(8,103) TIME( I ),X(3, I ) 

103 FORMAT(2E20. 10) 

104 CONTINUE 

DO 107 1=1, ND 

WRITE(8, 103) TIME( I ),X(4, I ) 

107 CONTINUE 

DO 122 1=1, ND 

WRITE(8, 103) TIME( I ),X(5, I ) 

122 CONTINUE 

DO 123 1=1, ND 

WRITE(8, 103) TIME{ I ),V0( I ) 

123 CONTINUE 

■ DO 124 1 = 1, ND 
WRITE(8, 103) TIME( I ),X(7, I ) 

124 CONTINUE 
109 CONTINUE 

DO 14 1=1, ND 
14 TVEC2( I )=VO( I ) 

STOP 
END 


SUBROUT I NE STRAN( TAU , PH I , THETA, A, B, NA, NB, MA, MB, MODE , NTERMS 
1 TOL,DIFMAX, ITER, IFLAG) 

REAL PH I ( NA, NA ) , THETAJ NA, NB ) , A( NA, NA ) , B( NA, NB ) 

1 WORK1 ( 10, 10),W0RK2( 1 0, 10 ) , DUMMY( 1 , 1 ) , FAC1 , CONI , FAC2, CON2 
C1 = 1 
TS=TAU 
DO 4 1=1, MA 
DO 2 J=1,MA 
WORKK l,J)=0 
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W0RK2( l,J)=0 
PHI( l,J)=0 
2 CONTINUE 
W0RK1(I,I)=1 
W0RK2( I . I )=TAU 
PHI ( I , I ) = 1 
4 CONTINUE 
DO 6 l=1,MA 
DO 6 J=1,MB 
6 THETA( I , J )=0 
Dl FMAX=1 . E8 
ND0=50 

IF(NTERMS.GT.O) NDO=NTERMS 
FAC1=1 
I FLAG=0 

DO 1000 1=1, NDO 

CALL MXMUL(DUMMY,WORK1,A, TS, 1. 1,10, 10, NA, NA, MA, MA, MA, MA, 1 ) 

FAC1=FAC1*I 

CON1=1 ,/FACI 

IF(NTERMS.EQ.O.AND, I .GE.4) CALL SERROR( PH I , WORK1 , CONI , NA, NA, 

1 10, 10,MA,MA,DIFMAX) 

CALL MXADD( DUMMY, PH I , WORK1 , Cl , CONI , 1 , 1 , NA, NA, 1 0, 1 0, MA, MA, 1 ) 
IF(MODE.EQ.2) GO TO 500 
FAC2=FAC1*( 1+1 ) 

C0N2=TAU/FAC2 

CALL MXAD0(0UMMY,W0RK2,W0RK1,C1,C0N2, 1,1, 10, 10,10, 10,MA,MA, 1 ) 
500 CONTINUE 
ITER=I 

IF(NTERMS.GT.O.OR. I.LT.4) GO TO 1000 
IF(DIFMAX.LE.TOL) GO TO 1100 
1000 CONTINUE 
1100 CONTINUE 

I F(MODE. EQ. 1 ) CALL MXMUL( THETA, W0RK2, B, Cl , NA, NB, 1 0, 1 0, NA, NB, 

1 MA,MA,MA,MB,2) 

IF( ITER.EQ.NDO.ANO.NTERMS.EQ^O) IFLAG=1 

RETURN 

END 


SUBROUTINE SERROR( AMX, BMX, CCC, IA,JA, IB,JB, I DO, JDO, D I FMAX ) 
DIMENSION AMX( IA,JA),BMX( IB,JB) 

Dl FMAX=1 . E-30 
DO 100 1=1, IDO 
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DO 50 J=1,JD0 

I F( AMX( I , J ) . EQ.0.0 ) GO TO 50 
CHANGE= ABS(BMX( I , J ) *CCC/AMX( I.J)) 
I F ( CHANGE . GT . 0 1 FMAX ) 0 1 FMAX=CHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 


SUBROUTINE FUTURE( X, V, PH I , THETA. LP, MP, LT, MT, IP.JP, IT.JT.MODE) 
DIMENSION X(MP),V(MT), PHI(LP.MP).THETA( LT, MT ) . TEMP( 20 ) 

DO 20 1=1, IP 
SUM=0 

DO 10 J=1,JP 
10 SUM=SUM+PHI( l,J)*X(J) 

20 TEMP( I ) = SUM 
DO 30 1=1, I P 
30 X( I )=TEMP( I ) 

IF(M0DE.EQ.2) RETURN 
DO 60 1 = 1 , IT 
SUM=0 

DO 50 J=1,JT 

50 SUM=SUM+THETA(I,J)*V(J) 

X( I )=X( I )+SUM 
60 CONTINUE 
RETURN 
END 


SUBROUTINE MXADD( RMX, AMX, BMX, ACC, BCG, IR,JR, IA,JA, IB,JB, IDO.JDO, 
1 MODE) 

DIMENSION AMX( IA,JA),BMX( IB,J6),RMX( IR.JR) 

I F( lA. LT. IDO.OR. JA. LT.JDO) GO TO 999 
I F( IB. LT. IDO.OR. JB. LT. JDO) GO TO 999 
GO TO ( 10, 100), MODE 
10 CONTINUE 

DO 50 1=1, IDO 
DO 50 J=1,JD0 

50 AMX( I , J )=AMX( I , J )*ACC+BMX( I , J )<»BCC 
GO TO 300 
100 CONTINUE 

IF( IR.LT. IDO.OR.JR.LT.JDO) 


GO TO 999 
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DO 200 1=1, IDO 
DO 200 J=1,JDO 

200 RMX( l,J)=AMX( I , J ) *ACC+BMX( l,J)*BCC 
300 RETURN 
999 CONTINUE 
RETURN 
END 


SUBROUTINE MXMUL( RMX, AMX, BMX, CCC. IR.JR. lA.JA. IB.JB, IDA.JDA, IDB.JDB 
1.MODE) 

DIMENSION AMX( IA,JA),BMX( IB,JB).RMX{ I R. JR ) . TEMP( 20 ) 

IF( I0A*IDB»JDA*JDB.GT. IA*IB*JA*JB) GO TO 999 
I F( JDB.GT. JDA) GO TO 999 
GO TO ( 10, 210), MODE 
10 DO 100 1=1 , IDA 
DO 20 L=1,J0A 
20 TEMPI L)=AMX( 1 ,L) 

DO 80 J=1,JDB 
SUM=0 

DO 1*0 K=1,JDA 

1*0 SUM=SUM+TEMP( K)*BMX( K, J) 

AMX( I , J ) = SUM*CCC 
80 CONTINUE 
100 CONTINUE 
GO TO 600 
210 CONTINUE 

IF( IR.LT. IDA. 0R.JR.lt. JDB) GO TO 999 
DO 400 1=1, IDA 
DO 380 J=1,JDB 
SUM=0 

DO 340 K=1,JDA 

340 SUM=SUM+AMX( I , K)*BMX( K, J ) 

RMX( I , J )=SUM*CCC 
380 CONTINUE 
400 CONTINUE 
600 RETURN 
999 CONTINUE 
RETURN 
END 
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c 

c ******************************************************************* 

C * THIS PROGRAM COMPUTES THE CLOSED-LOOP POLES OF A BUCK CONVERTER * 

C * WITH A SINGLE-STAGE INPUT FILTER, AS A FUNCTION OF THE INPUT * 

C * FILTER PARAMETERS. * 

C -(nut**************************************************************** 

C 

D I MENS ION A( 7 ) , B( 6 ) . CC1 ( 7 ) , CC2 ( 5 ) 

REAL NXEG, I XR. M, K1 . K2. K3 , K4. K5, K6, K7. L, LI , K8 

INTEGER NDEG, I ER 

COMPLEX Z(6), 21(5), Z2(6), 23(4) 

DATA VI , L,RO,RC,C,RL, VO/25. 0,230. OE-6, 0.2, 0.067. 300.0E-6, 10.0,20.0 
1/ 

DATA RN,R1 1 , R1 2 , R1 3 , R1 4, R4 . C2, CPI , M/0. 65. 33 . 3 E3, 16 . 7E3, 2. 0E5, 
147.0E3.40.7E3, 1 . OE- 10, 5600. OE-12,0. 88E-3/ 

C 

DO 2000 1=1,1000 
READ(5,*) RLI.LI.CI 
IF(RL1,LT.0. ) GO TO 2001 
0=0.8388579 
C 

Y1=RL*L1*C1*C*RC 

Y2=RL*L1*C1+C»RC*(C1*RL1*RL-(D**2)*L1 ) 

Y3=C1 *RL1 *RL- ( D*«2 )»L1 +C*RC*( RL-( 0**2 )*RL1 ) 

Y4=RL-(0**2)*RL1 

Y5=RL*L1*C1*C*RL 

Y6=RL*L1*C1+C*RL*(C1*RL1*RL-(D**2)*L1 ) 
Y7=C1*RL1*RL-(D**2)*L1+C*RL*(RL-(D**2)*RL1 ) 

Y8=L»C*RL*L1*C1 

Y9=( L*C*RL*C1*RL1 ) +( C*RL*L1 *C1 *( RO+RC+( L/(C*RL) ) ) ) 
Y10=L*C*RL+C*RL*C1*RL1*( RO+RC+( L/(C*RL) ) )+L1*C1*RL+( D**2 )*C*RL*L1 
Y11=C*RL*(RO+RC+( L/(C*RL) ) )+C1 *RL1 *RL+( 0**2 ) *( LI +C*RL*RL1 ) 
Y12=RL+(D**2)*RL1 
RX=(R11*R12)/(R11+R12) 

G=RX/R11 

K1=G/(R14+RX) 

K2=C2*( 1 .0+K1*R13) 

FM=(2.0*R4*CP1 )/(RN*M) 

K3=VO*RL*R4*FM 

K4=RN*L*C2*R13 

K5=V0*FM 

K6=CP1*C2*R13 

K7=D*RL*R4 


K8=C*(R0+RC+{L/(C*RL) ) ) 

C 

C ******************************************************** 

C * THIS CALCULATION IS FOR THE CLOSED LOOP POLES WITH * 

C * INPUT FILTER. * 

C 

A( 1 )=K7*K6*Y8 

A( 2 ) = K7*CP1*Y8+K7*K6*Y9+K5*KI|*Y5 

A( 3 )=K7*CP1*Y9+K7*K6*Y10+K3*K2*Y1+K5*RN*L*Y5+K5*K4*Y6 

A(4 )=K7*CP1*Y10+K7*K6*Y1 1+K3*K2*Y2+K3*K1*Y1+K5*RN*L*Y6+K5*K4*Y7 

A(5)=K7*CP1*Y11+K7*K6*Y12+K3*K1*Y2+K3*K2*Y3+K5*RN*L<»Y7+K5*K4*Y4 

A(6)=CP1*K7*Y12+K3*K1*Y3+K3*K2*Y4+RN*L*K5*Y4 

A( 7)=K3*K1*Y4 

NDEG=6 

CALL 2P0LR(A.NDEG.Z. lER) 

C 

WRITE(6, 101 ) RLI.LI.CI 

101 FORMAT! /, 10X, ' RL1 = ' , FI 5 . 1 0, ' L1 = ' , F15. 10, ’ C1 = ',F15.10) 
WRITE(6.103) 

103 FORMAT! /, 20X, ' THE CLOSED LOOP POLES WITH I. FILTER ARE-') 
WRITE!6, 102)!Z! IX), IX=1,6) 

102 FORMAT!/, 20X,2E15. 8) 

C 

B! 1 ) = K5*K4*Y5 

B! 2 )=K5*RN*L*Y5+K5*K4*Y6+K3*K2*Y1 

B! 3 )=K5*RN*L*Y6+K5*K4*Y7+K3*K2*Y2+K3*K1»Y1 

B!4)=K5*RN*L*Y7+K5*K4*Y4+K3*K1»Y2+K3*K2*Y3 

B! 5 )=RN*L*K5*Y4+K3*K1 *Y3+K3*K2*Y4 

B!6)=K3*K1*Y4 

NDEG=5 

CALL ZPOLR!B,NDEG,zr, lER) 

WRITE!6, 104) 

104 FORMAT! /, 20X, ' THE OPEN LOOP ZEROES ARE ') 

WRITE!6,102)!Z1! IX), IX=1,5) 

C 

CC1 ! 1 )=K7*K6*Y8 

CC1 ! 2 ) = K7*CP1 *Y8 + K7*K6»Y9 

CC1 ! 3 )=K7*CP1*Y9+K7*K6»Y10 

CC1 ! 4) = K7*CP1*Y10+K7*K6*Y1 1 

CC1 !5)=K7*CP1*Y1 1+K7*K6*Y12 

CC1!6)=CP1*K7*Y12 

CC1!7)=0.0 



CALL ZP0LR(CC1,6,Z2, lER) 

WRITE(6, 105) 

105 F0RMAT(/,20X. ' THE OPEN LOOP POLES ARE ') 

WR1TE(6,102)(Z21 IX). )X=1,6) 

C 

C ****************************************************************** 

C * THIS CALCULATION IS FOR THE CLOSED LOOP POLES WITHOUT INPUT * 

C * FILTER. * 

C ****************************************************************** 

c 

CC2( 1 ) = K7*K6*L*C 

CC2( 2 ) = K7*K6*K8+K7*CP1 *L*C+K5*K41*C*RL 
CC2( 3)=K7*CP1*K8+K7*K6+K3*K2*C*RC+K5*K4+K5*RN*L*C*RL 
CC2( 4 )=K7*CP1 +K3»K2+K3*K1*C*RC+RN*L*K5 
CC2(5)=K3*K1 

CALL ZPOLR(CC2,4,Z3, I ER) 

WRITE(6,785) 

785 FORMAT( /.20X, ' THE CLOSED LOOP POLES W/0 I. FILTER ARE-') 

WRITE(6, 102)(Z3( IX), IX=1,4) 

2000 CONTINUE 

2001 CONTINUE 
STOP 
END 
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*-H"IHHHt***********-(HHt******#*-H#*##H***'tt**-IHHHt#************l(-* 

STABILITY ANALYSIS OF A BUCK CONVERTER WITH INPUT FILTER * 

******-H***-tHHt********4Ht*******************4t***************H 


THIS PROGRAM COMPUTES THE EIGENVALUES OF A BUCK REGULATOR 
SYSTEM WITH INPUT FILTER. 

DIMENSION F1(6,6), F2(6.6),G1(6,4).G2(6,4).U(4), PHIF(6,6), 
1PHIN(6,6),DF(6,4),DN(6,4).X(6),ZT(6),AA(6.6),BB(6),TEMP1(6,4), 
2PHI FT(6,6).OFT(6,4),XA{6),ZN1(6).XB(6) 

DIMENSION PS I ( 6 , 6 ) . XN 1 ( 6 ) . fun 1 { 6 ) , FUN2 ( 6 ) , CON ( 4 ) , RVAL( 6.9) 
DIMENSION TEMP2(6,6),DX(6) 

DIMENSION WKAREA(70),WK(70) 

REAL L1,L,K1,K2 
COMPLEX W(6),Z(6,6),ZN 

DATA RL1 , L1.C1 , L, RO, RN, RC, C/O. 2000, 50 . OE-6, 220 . E-6, 230 . E-6, 0 . 2E0 
l,0.65EO,0.067EO,300.E-6/ 

DATA RL,R11,R12,R13,R14,C2,R4,CP1/10.0E0,33.3E3, 16.7E3,2.E5, 

147. E3, 0.01 E-8, 40. 7E3. 5600. E-12/ 

DATA ER , ET , V I , EO/6 . 7E0 , 7 . EO, 25 . EO, 20 . EO/ 

DATA VITON,Eq,EO/0.88E-3,0.2EO,O.7EO/ 

DATA NIT, EPS/90, 10. E-6/ 

DATA EPSl, ERROR/2. E-6, 0.5E-6/ 

DATA CON/0. 1E-1.0.05E-1 ,0.025E-1,0.0125E-1/ 

DATA TOLL/0.01/ 

K1=(R12/(R14*(R11+R12) ) )+( 1 . /R1 3 )-( RN/R4 ) 

K1=K1-(RN*RO*(RL+RC) )/(R4*RL*RC) 

K2=(RN/R4)-( 1 ./R13)-(R12/(R14*(R11+R12) ) ) 

TON=VITON/VI 
TOFFB=TON*( ( VI-EO)/EO) 

DO 1 1=1,6 
DO 2 J=1,4 
G1(I,J)=0. 

G2(I,J)=0. 

' DF( I . J )=0. 

DN( I , J )=0. 

CONTINUE 
DO 3 K=1,6 
F1( l,K)=0. 

F2(I,K)=0. 

AA( I , K)=0. 

PH1F(I,K)=0. 
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PHIN(l,K)=0. 

3 CONTINUE 

1 CONTINUE . 

DO 5 1=1,6 
ZT( I )=0. 

BB(l)=0. 

FUN1( I )=0. 

FUN2( I )=0. 

XA(l)=0. 

XB( I )=0. 

XN1(I)=0. 

ZN1( I )=0. 

5 CONTINUE 

F1( 1. 1 )=(-1 ,*RL1 )/L1 
F1(1,2)=-1.00/L1 
FI (2, 1 ) = 1 .00/C1 
Fl(2,3)=-1.00/C1 
F1(3,2)=1.00/L 

Fl( 3,3)=(-l .00*RO)/L+(-1 .00*RC»RL)/(L*(RL+RC) ) 
F1(3,4)=-1 .00*RL/(L*(RL+RC) ) 
F1(4,3)=RL/(C»(RC+RL)) 

F1(4,4)=-1 .00/(C*(RC+RL) ) 
F1(5,3)=(RL*RC)/(C2*R13*(RC+RL)) 
F1(5,4)=RL/(C2*R13*(RC+RL) ) 
F1(5,5)=-1.00/(C2»R13) 

F1(6,2)=(-1 .00»RN)/(CP1*R4) 
F1(6,3)=(-1.00*RL*RC*K1 )/(CP1*(RC+RL)) 
F1(6,4)=(RL*K2)/(CP1*(RL+RC) ) 
F1(6,5)=T.OO/(CP1*R13) 

Gl(1,1)=1. 00/LI 
G1(3,3)=-1.00/L 
Gl(6,2)=1 .00/(CP1*R14) 

G1(6,3)=RN/(CP1*R4) 

U( 1 )=VI 
U(2)=ER 
U(3)=EQ 
U(4)=ED 
DO 4 1=1,6 
DO 4 J=l,6 

4 F2( l,J) = F1( l,J) 

F2(2,3)=0.00 

F2(3,2)=0.00 

F2(6,2)=0.00 
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G2(1,1)=G1(1,1) 

G2(6,2)=G1(6,2) 

G2( 3,4)=-1 .00/L 
G2(6,4)=RN/(CP1*R4) 

C 

T=TON 

CALL STRAN(T,PHIN.DN.F1.G1,6,4,6,4,1,0,ERR0R.DIFMAX. ITER, I FLAG) 
T=TOFFB 

CALL STRAN(T,PHIF,DF,F2,G2,6,4,6,4,1,0,ERR0R,DIFMAX, ITER. I FLAG) 

************************************************* 

CALCULATION OF THE APPROXIMATE STEADY STATE * 
************************************************* 


DO 5000 1=1,6 
DO 5000 J=1,6 
AA(l,J)=0. 

DO 5001 K=1,6 

5001 AA( l,J)=AA( l,J)+PHIN( I , K )*PH I F( K, J ) 

AA( I ,J )=-1 .*AA( I ,J ) 

5000 CONTINUE 

DO 5002 1=1,6 

5002 AA( I , I ) = 1 ,+AA( 1,1) 

DO 5003 1=1,6 

DO 5003 J=1,4 
TEMPI ( I , J )=0. 

DO 5004 K=1,6 

5004 TEMP1( l,J)=TEMP1( l,J)+PHIN( l,K)*DF(K,J) 

5003 TEMP1( l,J)=TEMP1( l,J)+DN( l,J) 

DO 5005 1=1,6 

BB( I )=0. 

DO 5006 J=1,4 

5006 BB( I )=BB( I )+TEMP1( l,J)*U(J) 

5005 CONTINUE 
MM=1 
NN=5 

I AA=6 
I DGT=0 

CALL LEQT2F( AA,MM,NN, lAA.BB, IDGT.WKAREA, lER) 

DO 5007 1=1,5 

5007 X( I )iBB( I ) 

X(6)=ET-PH1 F(6, 1 )*X( 1 )-PHI F(6,2)*X(2)-PHIF(6, 3)*X( 3 )-PHI F(6,4)*X( 
14)-PHIF(6,5)*X(5)-DF(6. 1 )*U( 1 )-OF( 6,2 )*U( 2 )-DF( 6, 3 )*U( 3 )-DF( 6, 4 )* 



o o o o o 


2U(4) 

WRITE(6,5008) 

5008 FORMAT( 1H1 , 37X, ' APPROXIMATE STEADY STATE VALUES, X= ') 
WRITE(6.5009)(X( I ). 1=1,6) 

5009 F0RMAT(1H ,6F15.10) 

WRITE(6,5010) TOFFB 

5010 F0RMAT( 1H0, 37X, ' APPROXIMATE VALUE OF OFF TIME= ',F15.10) 

CALCULATION OF VhE EXACT STEADY STATE VALUES * 
************************************************* 


IT1=0 

900 DO 26 1=1,5 
ZT( I )=0.00 
DO 27 K=1,6 

27 ZT( I )=ZT( I ) + PHIF( l,K)*X(K) 

DO 28 J=1,4 

28 ZT( I )=ZT( I )+DF( I ,J )*U( J ) 

26 CONTINUE 

IT1=IT1+1 

SMAT=X(6)-PHIN(6, 1 )*ZT( 1 ) - PH I N( 6, 2 )*ZT( 2 )- PH I N( 6, 3 )*ZT( 3 ) - PH I N( 6, 4 
1 )*ZT(4)-PHIN(6,5)*ZT(5)-PHIN(6,6)*ET-ON(6,1 )*U(1 )-DN(6,2)*U(2)- 
2DN(6,3)*U(3)-DN(6,4)*U(4) 

I F( ABS( SMAT) . LE. EPS1 ) GO TO 70 
I F( IT1 .GT.60) GO TO 70 
DTOFF=0.01*TOFFB 
T=TOFFB+DTOFF 

CALL STRAN(T, PHI FT, OFT, F2, G2, 6, 4, 6, 4, 1 , 0, ERROR, D I FMAX, ITER, I FLAG) 
DO 29 1=1,6 
DO 29 J=1,6 
AA( I , J )=0.00 
DO 56 K=1,6 

56 AA( 1,J)=AA( l,J)+PHIN( I , K)*PH I FT( K, J ) 

AA( I ,J )=-1 .00*AA( I , J ) 

29 CONTINUE 

DO 31 1=1,6 

31 AA( I, I ) = 1.00+AA( I, I ) 

DO 32 1=1,6 
DO 32 J=1,4 
TEMPI ( I , J )=0.00 
DO 34 K= 1,6 

34 TEMP1( l,J) = TEMP1( l,J) + PHIN( I , K )*DFT( K, J ) 



o 


32 TEMP1( 1,J)=TEMP1( l,J)+DN( l,J) 

00 35 1=1,6 
BB( I )=0.00 
DO 36 J=1,4 

36 BB( I )=BB( I )+TEMP1( l,J)*U(J) 

35 CONTINUE 

MM=1 

NN=5 

IAA=6 

IDGT=0 

CALL LEQT2F(AA,MM,NN, lAA.BB, lOGT.WKAREA, lER) 

DO 37 1=1,5 

37 XN1(I)=BB(1) 

XN1 (6) = ET-PHI FT{6, 1 )*XN1 ( 1 )-PH I FT( 6, 2 )*XN1 ( 2 ) - PH I FT( 6, 3 )*XN1 ( 3 )- 
1 PHI FT(6,4 )*XN1 (4)-PHI FT(6,5)*XN1(5)-DFT(6, 1 )*U( 1 )-DFT(6,2)*U(2)- 
2DFT(6,3)*U(3)-DFT(6,4)»U(4) 

DO 38 1=1,5 
ZN1( I )=0.00 
DO 39 K=1,6 

39 ZN1( I )=ZN1( I )+PHIFT( l,K)»XN1(K) 

DO 40 J=1,4 

40 ZN1( I )=ZN1( I )+DFT( I ,J)*U(J) 

38 CONTINUE 

SMATN=XN1(6)-PHIN(6,1 )*ZN1( 1 )-PH I N( 6.2 )*ZN1 ( 2 )- PH I N( 6, 3 )*ZN1 ( 3 ) - 
1 PHIN(6,4)»ZN1(4)-PHIN(6,5)*ZN1(5)-PHIN(6.6)*ET-DN(6. 1 )*U( 1 )-DN( 6, 2 
2)*U(2)-DN(6.3)*U(3)-DN(6,4)*U(4) 

DSMAT=SMATN-SMAT 

TOFFB=TOFFB+( ( -1 . 00*SMAT )/( DSMAT/DTOFF) ) 

T=T0FFB 

CALL STRAN(T,PHIF,DF,F2,G2,6,4,6,4,1,0,ERROR,DIFMAX, ITER, I FLAG) 

DO 42 1=1,6 
DO 42 J=1,6 
AA( I , J )=0.00 
DO 43 K=1,6 

43 AA( l,J)=AA( l,J)+PHIN( I , K)*PH I F{ K, J ) 

AA( I , J )=-1 .00*AA( I , J ) 

42 CONTINUE 

DO 55 1=1 ,6 

55 AA( I, I )=AA( I, I )+1. 

DO 44 1=1,6 
DO 44 J=1,4 
TEMPI ( I , J )=0.00 
DO 46 K=1,6 
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46 TEMPI ( l,J)=TEMP1{ l,J)+PHIN{ l,K)*DF(K,J) 

44 TEMP1( l,J)=TEMP1( l,J)+DN( I.J) 

DO 47 1=1,6 
BB( I )=0.00 
DO 48 J=1,4 

48 BB( I )=BB( I )+TEMPl ( I , J )*U( J ) 

47 CONTINUE 
MM=1 
NN=5 
IAA=6 
IDGT=0 

CALL LEQT2F(AA.MM,NN, lAA.BB. lOGT.WKAREA, lER) 

DO 49 1=1,5 

49 X(1)=BB(1) 

X(6) = ET-PHI F(6, 1 )*X( 1 )-PHI F(6,2)*X(2)-PHIF(6. 3)*X( 3 )-PHI F(6,4 )* 
1X(4)-PHI F(6.5)*X(5)-DF(6.1 )*U( 1)-DF(6,2)*U(2)-DF(6,3)*U(3)- 
2DF(6,4)*U(4) 

GO TO 900 
70 CONTINUE 

WRITE(6,50) 

50 FORMAT( 1H0, 37X, ' EXACT STEADY STATE VALUES, X= ') 

WRITE(6,199)(X( I ), 1=1,6) 

199 F0RMAT(1H ,6F15.10) 

WRITE(6,198) TOFFB 

198 FORMAT( 1H0,37X, ' EXACT OFF TIME= ’,F15.10) 

DVAL=TON/( TON+TOFFB) 

WRITE(6,624) DVAL 

624 FORMAT( 1H0, 37X, ' VALUE OF 0=',F15.10) 

WRITE(6,51) IT1 

51 FORMAT( 1H0,37X, ' NO. OF ITERATIONS REQD.= ’,14) 

******************************************************** 

* CALCULATION OF THE MATRIX PS I AND ITS EIGENVALUES * 
##*####*******#*■»**##***#*#»«##»»****»*■»#*#*»#**»*#»##■*•# 


DO 2000 1=1,6 
DO 2010 IRICH=1,4 
DO 61 IZ=1,6 
XA( IZ)=X( IZ) 

61 XB(IZ)=X(IZ) 

DX( I )=CON( IRICH)*ABS(X( I ) ) 
IF(ABS(X( I )).LE.TOLL) DX(I)=0.01 
XA( I )=X( I )+DX( I ) 





XB( I )=X( I )-DX( I ) 

TFFC1=T0FFB 

IT=0 

T8=T0FFB 

CALL STRAN(T8.PHIF.0F,F2.G2,6,4,6,4,1,0.ERR0R,DIFMAX, ITER. I FLAG) 

1101 ZETA=PHIF(6.3)*XA(3)+PHIF(6,4)* 

1XA(4)+PHI F(6,5)*XA(5)+1 .0*XA(6)+DF(6,1 )*U(1 )+DF(6,2)*U(2)+ 
2DF(6.3)*U(3)+DF(6,4)*U(4)-ET 
I F(ABS(ZETA) . LE. EPS) GO TO 1102 
I F( IT.GE.NIT) GO TO 1003 
DTOFF=CON( IRICH)*TFFC1 
T=TFFC1+DT0FF 

CALL STRAN(T, PHI FT, OFT, F2, G2, 6. 4, 6, 4, 1 , 0, ERROR, D I FMAX, ITER, I FLAG) 
ZETAN=PH I FT( 6, 3 )*XA( 3 )+PH I FT( 6. 4 

1 )*XA(4)+PHI FT(6,5)*XA(5)+1 .0»XA(6)+DFT(6, 1 )*U( 1 )+DFT(6,2)* 
2U(2)+DFT(6,3)*U(3)+DFT(6,4)*U(4)-ET 
DZETA=ZETAN-ZETA 
SLOPE=DZETA/DTOFF 
TFFC1=TFFC1-(ZETA/SL0PE) 

T=TFFC1 

CALL STRAN(T,PHIF,DF,F2,G2,6,4,6.4,1,0,ERR0R,DIFMAX, ITER, I FLAG) 

IT=IT+1 

GO TO 1101 

1102 DO 13 J=1,6 
TEMP=0.00 
DO 14 K=1,6 

DO 81 MENT=1,5 
81 PHI F(MENT,6)=0.0 
PHI F(6,6)=1 .0 

14 TEMP=TEMP+PHIF(J,K)*XA(K) 

DO 15 K=1,4 

15 TEMP=TEMP+DF(J,K)*U(K) 

FUN1 ( J )=TEMP 

13 CONTINUE 

TFFC2=T0FFB 

IT=0 

T=T0FFB 

CALL STRAN(T,PHIF,DF,F2,G2,6,4,6,4,1,0,ERR0R,DIFMAX, ITER, I FLAG) 

1103 ZETA=PHI F(6,3)*XB(3)+PHI F(6,4)* 

1XB(4)+PHIF(6,5)*XB(5)+1 .0#XB(6)+DF(6.1 )*U( 1 )+DF( 6, 2 )*U( 2 )+ 
2DF(6,3)*U(3)+DF(6,4)*U(4)-ET 

I F( ABS(ZETA) . LE. EPS) GO TO 1104 
I F( IT.GE.NIT) GO TO 1003 
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17 CONTINUE 

WRITE(6,92) 

92 FORMAT( 1H0, 37X, ' MATRIX PSI— *) 

DO 30 1=1,6 

30 WRITE(6.93)(PSI( l,J),J=1,6) 

93 F0RMAT(6F15. 10) 

GO TO 1004 

1003 WRITE(6,100) I 

100 FORMAT(//,' CONVERGENCE NOT OBTAINED FOR X(l),l= ’,14) 

1004 CONTINUE 

CALL EIGRF(PSI,6,6,2,W,Z,6,WK, lER) 

WRITE(6,94) 

94 FORMAT( 1H0.37X, ' THE EIGENVALUES ARE— ') 

WRITE(6,103)(W( I ), 1=1,6) 

103 FORMAT( 1H ,6F15. 10) 

WRITE(6, 104) WK{ 1 ), lER 

104 FORMAT( 1H0,37X, ' CONVERGENCE TOLERANCE=* , F15. 10, ' IER=',I4) 
905 CONTINUE 

STOP 

END 


SUBROUTINE STRAN(TAU, PH I , THETA, A, B, NA, NB, MA, MB, MODE, NTERMS, 
1 TOL,DIFMAX, ITER, IFLAG) 

REAL PH I { NA, NA ) , THETA( NA, NB ) , A( NA, NA) , B( NA, NB ) , 

1 W0RK1 (10, 10),WORK2( 10, 1 0 ) , DUMMY( 1 , 1 ) , FAC1 , CONI , FAC2, C0N2 
C1 = 1 
TS=TAU 
DO 4 1=1, MA 
DO 2 J=1,MA 
W0RK1( l,J)=0 
W0RK2( I , J )=0 
PHI( I.J)=0 
2 CONTINUE 
W0RK1( I , I ) = 1 
W0RK2( I , I )=TAU 
PHI ( I , I )=1 
4 CONTINUE 
DO 6 l=1,MA 
DO 6 J=1,MB 
6 THETA( l,J)=0 
D I FMAX= 1 . E8 
NDO=50 
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IF(NTERMS.GT.O) NDO=NTERMS 
FAC1=1 
I FLAG=0 

DO 1000 1=1, NDO 

CALL MXMUL( DUMMY, WORK1 , A, TS, 1 , 1 , 10, 10, NA, NA, MA, MA, MA, MA, 1 ) 

FAC1=FAC1*I 

CON1=1 ./FAC1 

IF(NTERMS.EQ.O.AND. I .GE.4) CALL SERROR( PH I , WORK1 , CONI , NA, NA, 

1 10, lO.MA.MA.DI FMAX) 

CALL MXADD( DUMMY, PH I , WORK1 , Cl , CONI , 1 , 1 , NA, NA, 10, 10,MA,MA, 1 ) 

I F(MODE. EQ.2) GO TO 500 
FAC2=FAC1*( 1+1 ) 

CON2=TAU/FAC2 

CALL MXADD(DUMMY,W0RK2,W0RK1,C1,C0N2, 1, 1, 10, 10,10, 10,MA,MA, 1 ) 
500 CONTINUE 
ITER=I 

I F(NTERMS.GT.O.OR. I .LT.4) GO TO 1000 
I F( Dl FMAX. LE.TOL) GO TO 1100 
1000 CONTINUE 
1100 CONTINUE 

I F(MODE. EQ. 1 ) CALL MXMUL( THETA, W0RK2, B, Cl , NA, NB, 10, 10, NA, NB, 

1 MA MA MA MB 21 

I F( i TER. EQ. NDO. AND. NTERMS. EQ.O) I FLAG=1 
RETURN 
END 


SUBROUTINE SERROR( AMX, BMX, CCC, IA,JA, IB,JB, 1 DO, JDO, D I FMAX) 
DIMENSION AMX( IA,JA),BMX( IB,JB) 

Dl FMAX=1 . E-30 
DO 100 1=1 , IDO 
DO 50 J=1,JDO 

I F( AMX( I , J ) . EQ.O.O ) GO TO 50 
CHANGE= ABS{ BMX( I , J ) *CCC/AMX( I , J ) ) 

I F(CHANGE.GT.DI FMAX) D I FMAX=CHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 


SUBROUTINE MXADD( RMX, AMX, BMX, ACC, BCC, IR,JR, IA,JA, IB,JB, IDO, JDO, 
1 MODE) 
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DIMENSION AMX( IA,JA),BMX( IB,JB),RMX( IR,JR) 

I F( lA. LT. IDO.OR.JA. lt. JDO) GO TO 999 
IF( IB. LT. IDO. OR.JB.lt. JDO) GO TO 999 
GO TO ( 10, 100), MODE 
10 CONTINUE 

DO 50 1=1, IDO 
DO 50 J=1,JD0 

50 AMX( l,J)=AMX( I ,J )*ACC+BMX( l,J)*BCC 
GO TO 300 
100 CONTINUE 

IF( IR.LT. IDO. OR.JR.lt. JDO) GO TO 999 
DO 200 1=1, IDO 
DO 200 J=1,JD0 

200 RMX( l,J)=AMX( l,J)*ACC+BMX( l,J)«BCC 
300 RETURN 
999 CONTINUE 
RETURN 
END 


SUBROUTINE MXMUL( RMX, AMX, BMX, CCC, IR,JR, IA,JA, IB,JB, IDA,JDA, IDB,JDB 
1,M0DE) 

DIMENSION AMX( IA,JA),BMX( IB,JB),RMX( IR,JR),TEMP(20) 

IF( IOA*lDB»JDA*JDB.GT. IA*IB»JA*JB) GO TO 999 
I F( JDB.GT.JDA) GO TO 999 
GO TO ( 10, 210), MODE 
10 DO 100 1=1, IDA 
DO 20 L=1,JDA 
20 TEMP( L)=AMX( I , L) 

DO 80 J=1,JDB 
SUM=0 

DO i*0 K=1,JDA 

40 SUM=SUM+TEMP( K)*BMX( K,J ) 

AMX( I , J )=SUM*CCC 
80 CONTINUE 
100 CONTINUE 
GO TO 600 
210 CONTINUE 

I F( I R. LT. IDA.OR. JR. LT. JDB) GO TO 999 
DO 400 1=1, IDA 
DO 380 J=1,JDB 
SUM=0 

DO 340 K=1,JDA 


340 SUM=SUM+AMX( I , K)*BMX( K, J ) 
RMX( I , J )=SUM*CCC 
380 CONTINUE 
400 CONTINUE 
600 RETURN 
999 CONTINUE 
RETURN 
END 
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**************************************************** 

* STABILITY ANALYSIS OF A BUCK CONVERTER * 

*************************************************** 


THIS PROGRAM COMPUTES THE EIGENVALUES OF A BUCK 
REGULATOR SYSTEM WITH INPUT FILTER AND 
FEEDFORWARD. 

DIMENSION F1(7,7), F2(7,7).G1(7.4),G2(7,4),U(4). PHI F(7,7), 
1PHIN(7,7).DF(7,4),DN(7.4).X(7).ZT(7),AA(7.7).BB(7),TEMP1(7,4), 

2PHI FT(7,7).DFT(7.4).XN1(7).ZN1{7) 

D I MENS ION PS I ( 7, 7 ) , XA( 7 ) , FUN1 ( 7 ) , FUN2( 7 ) , XB( 7 ) , CON( 4 ) , RVAL( 7,9) 
DIMENSION TEMP2(7.7),DX(7) 

DIMENSION WKAREA(90),WK(90) 

REAL L1,L,K1,K2,K3 
COMPLEX W(7),Z(7,7),ZN 

DATA RL1 , LI ,C1 . L. RO. RN, RC, C/O . 20E0, 1 . 425E-3 . 220 . E-6,230. E-6,0.2E0, 
10.65E0,0.067E0.300.E-6/ 

DATA RL,R11,R12,R13,R14,C2,R4,CP1/10.0E0,33.3E3, 16.7E3,2.E5, 

1 47. E3, 0. 0 IE-8, 40. 7E3, 5600. E-12/ 

DATA ER,ET,VI,EO/6.7E0,7.EO,25.EO,2O.EO/ 

DATA VITON,EQ,ED/0.88E-3,0.7E0,0.2E0/ 

DATA NIT, EPS/150,5. E-6/ 

DATA EPS1/0.5E-6/ 

DATA ERROR/0.5E-6/ 

DATA CON/0.2E-1,0.1E-1.0.05E-1,0.025E-1/ 

DATA CF,RF1,RF2,TOLL/2.7E-5,5. 1E3,220.0EO,0. 11/ 

K1=(R12/(R14*(R1 1+R12) ) )+( 1 . /R1 3 )-( RN/R4 ) 

K1=Kl-( RN*RO*( RL+RC) )/(R4*RL*RC) 

K2=(RN/R4)-( 1 ./R13)-(R12/(R14*{R11+R12) ) ) 

K3=RF2/( RF1+RF2) 

TON=VITON/VI 
TOFFB=TON*( (VI-EO)/EO) 

DO 1 1=1,7 
DO 2 J=1,4 
G1 ( I , J )=0. 

G2(I,J)=0. 

DF(l,J)=0. 

DN( l,J)=0. 

2 CONTINUE 

DO 3 K=1,7 
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F1(I,K)=0. 

F2( l,K)=0. 

AA(l,K)=0. 

PHIF(l,K)=0. 

PHIN(l,K)=0. 

3 CONTINUE 

1 CONTINUE 

DO 5 1=1,7 
2T( I )=0. 

BB( I )=0. 

FUN1( I )=0. 

FUN2( I )=0. 

XA( I )=0. 

XB( I )=0. 

XN1( I )=0. 

ZN1( I )=0. 

5 CONTINUE 

F1(1,1 )=(-!. *RL1 )/L1 

F1(1,2)=-1.00/L1 

F1(2,l )=1. 00/Cl 

F1(2,3) = -1./(C1*(RF1+RF2) ) 

Fl(2,4)=-1 .00/Cl 
F1(3,1 )=1./C1 

Fl(3,3)=-1 .*(C1+CF)/(C1*CF*(RF1+RF2)) 
Fl(3,4)=-1 ./Cl 
FI (4,2)=1 .00/L 

FI (4,4) = ( -1 .00*RO)/L+(-1 .00*RC*RL)/(L*(RL+RC) ) 
FI (4,5)=-l .00*RL/( L*( RL+RC) ) 
F1(5,4)=RL/(C*(RC+RL) ) 
F1(5.5)=-1.00/(C*(RC+RL)) 
F1(6,4)=(RL*RC)/(C2»R13*(RC+RL)) 
F1(6,5)=RL/(C2*R13*(RC+RL) ) 

F1(6,6)=-1 .00/(C2*R13) 

F1(7, 1 )=-1 .*K3/C1 

F1(7,2)=(-1 .00*RN)/(CP1*R4) 

F1(7,3)=K3*(C1+CF)/(C1»CF*(RF1+RF2)) 

FI (7,4) = ( K3/C1 )-( ( RL*RC*K1 )/( CPI *( RL+RC) ) ) 

F1(7,5)=(RL*K2)/(CP1*(RL+RC)) 

F1(7,6)=1.00/(CP1*R13) 

G1 ( 1, 1 )=1 .00/LI 
G1(4,3)=-1.00/L 
Gl(7,2)=1 .00/(CP1*R14) 

G1(7,3)=RN/(CP1*R4) 
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U(1)=VI 
U(2)=ER 
U(3)=EQ 
U(4)=ED 
DO 4 1=1,7 
DO 4 J=1,7 
4 F2( l,J) = F1( I.J) 

F2(2.4)=0. 

F2(3.4)=0. 

F2(4.2)=0. 

F2(7.2)=0. 

F2(7,4)=(-1 .»RL*RC*K1 )/( CP1*( RL+RC) ) 

G2(1,1)=G1(1,1) 

G2(7.2)=G1(7.2) 

G2(4,4)=-1 .00/L 
G2(7.4)=RN/(CP1*R4) 

C 

T=T0N 

CALL STRAN(T,PH1N.DN,F1,G1,7,4.7,4,1,0,ERR0R,DIFMAX, ITER, I FLAG) 
T=TOFFB 

CALL STRAN(T,PHIF,DF,F2,G2,7.4,7.4,1,0,ERR0R,DIFMAX, ITER, I FLAG) 

*#»*****************#**###**#****##***#*#*##******* 

* CALCULATION OF THE APPROXIMATE STEADY STATE * 
w************************************************** 

DO 5000 1=1,7 
DO 5000 J=1,7 
AA(l,J)=0. 

DO 5001 K=1,7 

5001 AA( l,J)=AA( l,J)+PHIN( I , K ) *PH I F( K, J ) 

AA( l,J)=-1.*AA( l,J) 

5000 CONTINUE 

DO 5002 1=1,7 

5002 AA( I, I ) = 1.+AA( 1,1) 

DO 5003 1=1,7 

DO 5003 J=1,4 
TEMPI ( I , J )=0. 

DO 5004 K=1,7 

5004 TEMP1( l,J)=TEMP1( l,J)+PHIN( l,K)*DF(K,J) 

5003 TEMPI ( l,J)=TEMP1( l,J)+DN( l,J) 

DO 5005 1=1,7 

BB(l)=0. 
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DO 5006 J=1,4 

5006 BB( I )=BB( I )+TEMP1( l,J)*U(J) 

5005 CONTINUE 

MM=1 
NN=6 
I AA=7 
IDGT=0 

CALL LEQT2F(AA,MM,NN, lAA.BB, IDGT.WKAREA, I ER ) 

DO 5007 1=1,6 

5007 X( I )=BB( I ) 

X(7) = ET-PH1F(7.1 )*X( 1 )- PH I F( 7 , 2 )*X( 2 )-PH I F( 7 , 3 )*X( 3 ) - PH I F( 7 , 4 )*X( 
14)-PHIF(7.5)*X(5)-DF(7.1 )*U(1)-DF(7.2)*U(2)-DF(7,3)*U(3)-DF(7.4)* 
2U(4)-PHIF(7,6)*X(6) 

WR1TE(6,5008) 

5008 FORMAT( 1H1,37X. ' APPROXIMATE STEADY STATE VALUES, X= ') 
WRITE(6,5009)(X( I ), 1 = 1,7) 

5009 FORMAT( 1H ,7F15. 10) 

WRITE(6,5010)TOFFB 

5010 FORMAT( 1HO,37X, ' APPROXIMATE VALUE OF OFF TIME=' , F15. 10) 

*************************************************** 

* CALCULATION OF THE EXACT STEADY STATE VALUES * 
***■»*###*#****»******#*#******#***#*■»**»#**#******* 


IT1=0 
900 DO 26 1=1,6 
ZT( I )=0.00 
DO 27 K=1,7 

27 ZT( I )=2T( I ) + PHIF( l,K)*X(K) 

DO 28 J=1,4 

28 ZT( I )=ZT( I )+DF( I ,J )*U( J ) 

26 CONTINUE 

IT1=IT1+1 

SMAT=X(7)-PHIN(7,1 )*ZT( 1 ) - PH I N( 7, 2 )*ZT( 2 ) - PH I N( 7, 3 )*ZT( 3 ) - PH I N( 7, 4 
1 )>^ZT(4)-PHIN(7,5)<fZT(5)-PHIN(7,7)*ET-ON(7,1 )*U( 1)-DN(7,2)*U(2)- 
2DN(7,3)*U(3)-0N(7,4)*U(4)-PHIN(7,6)*ZT(6) 

IF(ABS(SMAT).LE.EPS1) GO TO 70 
I F( IT1 .GT.60) GO TO 70 
DTOFF=0.01*TOFFB 
T=TOFFB+DTOFF 

CALL STRAN(T, PH I FT, DFT, F2, G2, 7, 4, 7, 4, 1 , 0, ERROR, D I FMAX, ITER, I FLAG) 
DO 29 1=1,7 
DO 29 J=1,7 



AA( l,J)=0.00 
DO 56 K=1,7 

56 AA( l,J)=AA( l,J)+PHIN( I , K)*PH I FT( K, J ) 

AA( 1,J)=-1.00*AA( I.J) 

29 CONTINUE 

DO 31 1=1,7 

31 AA( I . 1 )=1 ,00+AA( 1,1) 

DO 32 1=1,7 

DO 32 J=1,4 
TEMP1( l,J)=0.00 
DO 34 K=1,7 

34 TEMPI! I.J)=TEMP1( I , J ) +PH I N( 1 , K )*DFT( K, J ) 

32 TEMPI! I,J)=TEMP1( l,J)+DN( l,J) 

DO 35 1=1,7 

BB! I )=0.00 
DO 36 J=1,4 

36 BB! I )=BB! I )+TEMP1! I,J)*U!J) 

35 CONTINUE 
MM=1 
NN=6 
IAA=7 
IDGT=0 

CALL LEQT2F!AA,MM.NN, 1AA,BB, IDGT,WKAREA, lER) 

DO 37 1=1,6 

37 XN1!I)=BB!I) 

XN1!7)=ET-PHI FT!7, 1 )*XN1! 1 )-PHIFT!7,2)*XN1!2)-PHIFT!7,3)*XN1!3)- 
1PHI FT!7,4)*XN1!4)-PHIFT!7,5)*XN1!5)-DFT!7,1 )*U! 1 )-DFT!7,2)*U!2)- 
2DFT!7,3)*U! 3)-DFT!7,4)*U!4)-PHIFT!7,6)*XN1!6) 

DO 38 1=1,6 
ZN1!I)=0.00 
DO 39 K=1,7 

39 ZN1! I )=ZN1! 1 )+PHlFT! I,K)*XN1!K) 

DO 40 J=1,4 

40 ZN1 ! I )=ZN1 ! 1 )+DFT! I ,J )*U! J ) 

38 CONTINUE 

SMATN=XN1 !7)-PHlN! 7, 1 )*ZN1!1)-PHIN!7,2)*ZN1!2)-PHIN!7,3)»ZN1!3)- 
1PHIN!7,4)*ZN1!4)-PHIN!7,5)*ZN1!5)-PH1N|7,7)*ET-DN!7,1 )*U!1 )-DN!7,2 
2)*U!2)-DN!7,3)*U!3)-DN!7,4)*U!4)-PHIN!7,6)*ZN1!6) 

DSMAT=SMATN-SMAT 

TOFFB=TOFFB+! !-1 ,00*SMAT)/!DSMAT/DTOFF) ) 

C 

T=TOFFB 

CALL STRAN!T, PHI F,DF, F2, G2, 7, 4, 7, 4, 1 , 0, ERROR, D I FMAX, ITER, I FLAG) 
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DO 42 1=1,7 
DO 42 J=l,7 
AA( I . J )=0.00 
DO 43 K=1,7 

43 AA( l,J)=AA( l,J) + PHIN( I.K)*PHIF(K,J) 

AA( I ,J )=-1 .00*AA( I ,J ) 

42 CONTINUE 

DO 55 1=1,7 

55 AA( I, I )=AA( I, I )+1. 

DO 44 1=1,7 
DO 44 J=1,4 
TEMPK l,J)=0.00 
DO 46 K= 1,7 

46 TEMPI ( l,J)=TEMP1( l,J)+PHIN( l,K)*DF(K,J) 

44 TEMPK l,J)=TEMP1( l,J)+DN( l,J) 

DO 47 1=1,7 

BB(l)=0.00 
DO 48 J=1,4 

48 BB( I )=BB( I )+TEMP1( l,J)*U(J) 

47 CONTINUE 
MM=1 
NN=6 
IAA=7 

I DGT=0 

CALL LEQT2F(AA,MM,NN, IAA,BB, 1DGT,WKAREA, lER) 

DO 49 1=1,6 

49 X( I ) = BB( I ) 

X(7)=ET-PHIF(7,1 )*X(1 )-PHIF(7,2)*X(2)-PHIF(7,3)*X(3)-PHIF(7,4)* 
1X(4)-PHIF(7,5)*X(5)-DF(7,1)*U(1)-DF(7,2)*U(2)-DF(7,3)*U(3)- 
2DF(7,4)*U(4)-PHI F(7,6)*X(6) 

GO TO 900 
70 CONTINUE 

WRITE(6,50) 

50 FORMAT( 1H0,37X, ' EXACT STEADY STATE VALUES, X=') 

, WRITE(6, 199)(X( I ), 1 = 1,7) 

199 F0RMAT(1H .7F15.10) 

WRITE(6,198) TOFFB 

198 FORMAT( 1H0,37X, ' EXACT OFF TIME= ’,F15.10) 

WRITE(6,51) IT1 

51 FORMAT( 1H0, 37X, ' NO. OF ITERATIONS REQD,= ’,14) 

-K- -H- # -H- # *H> -H- -M- # # -(f # -H- ## -H- -H -M- ## -M- ««#«««##«##### ^ -H- H 4f 

* CALCULATION OF THE MATRIX PSl AND ITS EIGENVALUES * 

***-H-4f**'M”H-****‘H-*-»f*4f*4f**********###***-»HHf*#****'H-*-H-#-M-*'H-*Hf-M‘* 



DO 2000 1=1,7 
DO 2010 IRICH=1,4 
DO 61 IZ=1,7 
XA(IZ)=X(IZ) 

61 XB(IZ)=X(IZ) 

DX( I )=CON( 1RICH)*ABS(X( I ) ) 

I F(ABS(X( 1 ) ) . LE.TOLL) DX( I )=CON( I RICH) 

XA( I )=X( I )+DX( I ) 

XB( I )=X( I )-DX( I ) 

TFFC1=T0FFB 

IT=0 

T=TOFFB 

' CALL STRAN(T, PHIF,DF, F2, G2, 7, 4, 7, 4, 1 . 0, ERROR, D I FMAX, ITER, I FLSG) 

1101 ZETA=0.0*XA( 1 )+0.0*XA(2)+PHIF(7,3)*XA(3)+PHIF(7,4)* 

1XA(4)+PHI F(7,5)*XA(5)+1 .0*XA(7)+0.0*U(1 )+DF(7,2)*U(2)+ 
2DF(7,4)*U(4)-ET+PH1 F(7,6)*XA(6) 

I F( ABS(ZETA) . LE. EPS) GO TO 1102 
IF( IT.GE.NIT) GO TO 1003 
DTOFF=CON( IRICH)*TFFC1 
T=TFFC1+DT0FF 

CALL STRAN(T, PHI FT,DFT, F2, G2, 7, 4, 7, 4, 1 , 0, ERROR, D I FMAX, ITER, I FLAG) 
ZETAN=0.0*XA( 1 )+0 . 0*XA{ 2 )+PH I FT( 7, 3 )*XA( 3 )+PH I FT( 7, 4 
1 )»XA(4)+PHI FT(7,5)*XA(5)+1 .0*XA(7)+0.0*U(1 )+DFT(7,2)* 
2U(2)+DFT(7,4)*U(4)-ET+PHI FT(7,6)*XA(6) 

DZETA=ZETAN-ZETA 

SLOPE=DZETA/DTOFF 

TFFC1=TFFC1-(ZETA/SL0PE) 

T=TFFC1 

CALL STRAN(T,PHIF,DF,F2,G2,7,4,7,4,1,0,ERR0R,DIFMAX, ITER, IFLAG) 

IT=IT+1 

GO TO 1101 

1102 DO 13 J=1,7 

TEMP=0.00 . 

DO 14 K= 1.7 
DO-81 MENT=1,6 

81 PHI F(MENT,7)=0.0 
PHI F(7,7) = l .0 
PHI F(7, 1 )=0.0 
PHI F(7,2)=0.0 

14 TEMP=TEMP+PHI F( J, K)*XA( K) 

DO 15 K=1,4 

DF(7, 1 )=0.0 

15 TEMP=TEMP■^DF(J,K)»U(K) 
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FUN1(J)=TEMP 
13 CONTINUE 

TFFC2=TOFFB 

IT=0 

T=TOFFB 

CALL STRAN(T,PHIF,DF.F2,G2,7,4,7,4,1,0.ERR0R.DIFMAX, ITER. IFLAG) 

1103 ZETA=0.0*XB(1 )+0.0*XB(2)+PHIF(7.3)*XB(3)+PHIF(7,4)* 
1XB(4)+PHIF(7,5)*XB(5)+1 .0*XB(7)+0.0*U(1)+DF(7.2)*U(2)+ 
2DF(7.4)*U(4)-ET+PHIF(7.6)*XB(6) 

I F(ABS(2ETA) . LE. EPS) GO TO 1104 
I F( IT.GE.NIT) GO TO 1003 
DTOFF=CON( IRICH)*TFFC2 
T=TFFC2+DT0FF 

CALL STRAN(T. PH I FT. OFT, F2, G2, 7. 4, 7, 4. 1 , 0, ERROR, D I FMAX. ITER, I FLAG) 
2ETAN=0. 0*XB( 1 )+0. 0*XB( 2 )+PH I FT( 7. 3 )*XB( 3 )+PH I FT( 7, 4 
1 )*XB(4)+PHI FT(7,5)*XB(5)+1 .0*XB(7)+0.0*U( 1 )+DFT(7,2)* 
2U(2)+DFT(7,4)*U(4)-ET+PHI FT(7,6)*XB(6) 

02ETA=2ETAN-2ETA 

SL0PE=D2ETA/DT0FF 

TFFC2=TFFC2-(2ETA/SL0PE) 

T=TFFC2 

- CALL STRAN(T.PHIF.DF,F2,G2,7,4,7,4,1,0.ERR0R,DIFMAX. ITER, IFLAG) 
IT=IT+1 
GO TO 1103 

1104 DO 62 J=1,7 
TEMP=0.00 
DO 63 K=1,7 

DO 82 MENT=1,6 
82 PHI F(MENT,7)=0.0 
PHIF(7,7)=1.0 
PHI F(7, 1 )=0.0 
PHI F(7,2)=0.0 

63 TEMP=TEMP+PHI F( J, K)*XB( K) 

DO 64 K= 1,4 

DF(7,.1 )=0.0 

64 TEMP=TEMP+DF(J,K)*U(K) 

FUN2( J )=TEMP 

62 CONTINUE 

DO 2011 J=1,7 

2011 RVAL(J. IRICH)=(FUN1(J)-FUN2(J))/(2.0*DX(I )) 

2010 CONTINUE 

DO 65 J=1,7 

RVAL( J,5)=(RVAL( J,2)-0.25*RVAL(J, 1 ) )/0.75 



RVAL(J.6)=(RVAL(J,3)-0.25*RVAL(J.2))/0.75 

65 RVAL(J,7)=(RVAL(J,4)-0.25*RVAL(J,3))/0.75 
DO 66 J=1,7 

RVAL( J . 8 ) = ( RVAL( J , 6 ) -0 . 0625*RVAL( J . 5 ) )/0 . 9375 

66 RVAL(J,9)=(RVAL(J,7)-0.0625*RVAL(J,6) )/0.9375 
DO 67 J=1,7 

67 TEMP2(J, I )=(RVAL(J,9)-0.015625*RVAL(J,8))/0. 984375 
2000 CONTINUE 

DO 17 1=1,7 
DO 17 J=1,7 
TEMP=0.00 
DO 19 K=1,7 

19 TEMP=TEMP+PHIN( I , K)*TEMP2( K, J ) 

PSI ( 1 . J )=TEMP 
17 CONTINUE 

WRITE(6,92) 

92 FORMAT( 1H0.37X, ' MATRIX PSI= ') 

DO 30 1=1,7 

30 WRITE(6,93)(PSI( l,J),J=1,7) 

93 EORMAT(7F15. 10) 

GO TO 1004 

1003 WRITE(6,100) I 

100 F0RMAT(//,' CONVERGENCE NOT OBTAINED FOR X(l),l= ’,14) 

1004 CONTINUE 

CALL EIGRF( PS I , 7, 7, 2, W, Z, 7, WK, I ER ) 

WRITE(6,94) 

94 FORMAT( 1H0, 37X, ' THE E I GENVALUES ARE— ' ) 

WRITE(6, 103)(W( I ), 1 = 1,7) 

103 F0RMAT(1H ,7F15.10) 

WRITE(6, 104) WK( 1 ), I ER 

104 FORMAT! 1H0, 37X, ' CONVERGENCE TOLERANCE= ', FI 5 . 1 0, ' I ER= ',14) 
905 CONTINUE 

STOP 

END 

C 

C 

SUBROUT I NE STRAN( TAU, PH I , THETA, A, B, NA, NB, MA, MB, MODE, NTERMS, 

1 TOL,DI FMAX, ITER, I FLAG) 

REAL PH I ( NA, NA ) . THETA( NA, NB ) , A( NA. NA ) , B( NA, NB ) , 

1 W0RK1 (10,10), W0RK2( 10,10), DUMMY! 1,1), FAC1 , CONI , FAC2, C0N2 
C1 = 1 
TS=TAU 
DO 4 1 = 1, MA 
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DO 2 J=1,MA 
W0RK1(I,J)=0 
W0RK2( I , J )=0 
PHI( I.J)=0 
2 CONTINUE 
WORK1( I. I ) = 1 
W0RK2( I , I )=TAU 
PH1(1,1)=1 
4 CONTINUE 
DO 6 l=1,MA 
DO 6 J=1,MB 
6 THETA( l,J)=0 
Dl FMAX=1 . E8 
NDO=50 

IF(NTERMS.GT.O) NDO=NTERMS 
FAC1=1 
I FLAG=0 

DO 1000 1=1, NDO 

CALL MXMUL( DUMMY, W0RK1 , A, TS, 1 , 1 , 10, 10, NA, NA, MA, MA, MA, MA. 1 ) 

FAC1=FAC1*I 

C0N1=1 ./FAC1 

IF(NTERMS.EQ.0.AN0. I .GE.4) CALL SERROR( PH I , WORKl , CONI , NA, NA, 

1 10, 10,MA,MA,0I FMAX) 

CALL MXA0D< DUMMY, PH I , WORKl , Cl , CONI , 1 , 1 , NA, NA, 1 0, 1 0, MA, MA, 1 ) 
IF(MODE.EQ.2) GO TO 500 
FAC2=FAC1*( I+1 ) 

C0N2=TAU/FAC2 

CALL MXADD( DUMMY, W0RK2, WORKl, Cl, C0N2, 1,1, 10,10, 10, 10,MA,MA, 1 ) 
500 CONTINUE 
ITER=I 

IF(NTERMS.GT.O.OR. I .LT.4) GO TO 1000 
IF(DIFMAX.LE.TOL) GO TO 1100 
1000 CONTINUE 
1100 CONTINUE 

I F(MODE. EQ. 1 ) CALL MXMUL( THETA,WORK2, B, Cl , NA, NB, 10, 10, NA, NB, 

1 MA,MA,MA,MB,2) 

IF( ITER.EQ.NDO.AND.NTERMS.EQ.O) I FLAG=1 
RETURN 
END 
C 
C 

SUBROUTINE SERROR( AMX, BMX, CCC, IA,JA, IB,JB, I DO, JDO, D I FMAX ) 
DIMENSION AMX( lA, JA),BMX{ IB,JB) 
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Dl FMAX=1 . E-30 
DO 100 1=1, IDO 
DO 50 J=1,JD0 

I F(AMX( I , J ) . EQ.0.0 ) GO TO 50 
CHANGE= ABS( BMX( I . J )*CCC/AMX( I , J ) ) 
I F(CHANGE.GT.DI FMAX) D I FMAX=GHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 


SUBROUTINE MXADD( RMX, AMX, BMX. ACC, BCC, IR,JR, 1A,JA, IB,JB, IDO,JDO, 
1 MODE) 

DIMENSION AMX( IA,JA),BMX( IB, JB),RMX( IR,JR) 

IF( lA.LT. IDO. 0 R.JA.lt. JDO) GO TO 999 
IF( IB. LT. IDO. OR.JB.lt. JDO) GO TO 999 
GO TO ( 10, 100), MODE 
10 CONTINUE 

DO 50 1=1, IDO 
DO 50 J=1,JD0 

50 AMX( l,J)=AMX( l,J)*ACC+BMX( I,J)*BCC 
GO TO 300 
100 CONTINUE 

IF( IR.LT. IDO. OR.JR.lt. JDO) GO TO 999 
DO 200 1=1, IDO 
DO 200 J=1 , JDO 

200 RMX( l,J)=AMX( I , J ) *ACC+BMX( l,J)*BCC 
300 RETURN 
999 CONTINUE 
RETURN 
END 


SUBROUTINE MXMUL( RMX, AMX, BMX, CCC, IR,JR, IA,JA, IB,JB, IDA,JDA, IDB,JDB 
1,M0DE) 

DIMENSION AMX( I A, JA ) , BMX( I B, JB) , RMX( I R, JR ) , TEMP( 20 ) 

I F( IDA*IDB*JDA*JDB.GT. IA*IB*JA*JB) GO TO 999 
I F( JDB.GT.JDA) GO TO 999 
GO TO ( 10, 210), MODE 
10 DO 100 1=1, IDA 
DO 20 L=1,JDA 
20 TEMP( L)=AMX( 1 , L) 
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DO 80 J=1,JDB 
SUM=0 

DO 40 K=1,JDA 

40 SUM=SUM+TEMP( K)*BMX( K.J ) 

AMX( I , J )=SUM*CCC 
80 CONTINUE 
100 CONTINUE 
GO TO 600 
210 CONTINUE 

IF( IR.LT. IDA.ORi.JR.LT.JDB) GO TO 999 
DO 400 1=1, IDA 
DO' 380 J=1.JDB 
SUM=0 

DO 340 K=1,JDA 

340 SUM=SUM+AMX( I , K)*BMX( K.J ) 

RMX( I ,J )=SUM*CCC 
380 CONTINUE 
400 CONTINUE 
600 RETURN 
999 CONTINUE 
RETURN 
END 



DTOFF=CON( IRICH)*TFFC2 
T=TFFC2+DT0FF 

CALL STRAN(T,PHIFT,DFT.F2.G2,6.4.6,4,1,0.ERR0R,DIFMAX, ITER, IFLAG) 
ZETAN=PHIFT(6,3)*XB(3)+PHIFT(6,4 

1 )*XB(4)+PHI FT(6,5)*XB(5)+1 .0*XB(6)+DFT(6.1 )*U( 1 )+DFT(6.2)* 
2U(2)+DFT(6,3)*U(3)+DFT(6,4)*U(4)-ET 
DZETA=ZETAN-2ETA 
SLOPE=DZETA/DTOFF 
TFFC2=TFFC2-(ZETA/SL0PE) 

T=TFFC2 

CALL STRAN( T, PH I F. OF, F2, G2, 6, 4. 6, 4, 1 , 0, ERROR, D I FMAX, I TER, I FLAG ) 
IT=IT+1 
GO TO 1103 
1104 DO 62 J=1,6 
TEMP=0.00 
DO 63 K=1,6 
DO 82 MENT=1,5 
82 PHI F(MENT,6)=0,0 
PHI F(6,6)=l .0 

63 TEMP=TEMP+PHI F(J,K)*XB( K) 

DO 64 K=1,4 

64 TEMP=TEMP+DF(J,K)*U( K) 

FUN2(J )=TEMP 

62 CONTINUE 

DO 2011 J=1,6 

2011 RVAL(J, IRICH)=( FUN1 ( J ) - FUN2( J ) )/( 2, 0*DX( I ) ) 

2010 CONTINUE 

DO 65 J=l,6 

RVAL(J,5)=(RVAL(J,2)-0.25*RVAL(J,1 ) )/0.75 
RVAL( J , 6 ) = ( RVAL( J , 3 ) -0 . 25*RVAL( J , 2 ) )/0 . 75 

65 RVAL(J,7)=(RVAL(J,4)-0.25*RVAL( J,3))/0.75 
DO 66 J=1,6 

RVAL(J,8)=(RVAL(J,6)-0.0625*RVAL(J,5) )/0.9375 

66 RVAL(J,9)=(RVAL(J,7)-0.0625*RVAL(J,6))/0.9375 
DO 67 J=l,6 

67 TEMP2(J, I ) = (RVAL(J,9)-0.015625*RVAL(J,8))/0. 984375 
2000 CONTINUE 

DO 17 1=1,6 
DO 17 J=1,6 
TEMP=0.00 
DO 19 K=l,6 • 

19 TEMP=TEMP+PHIN( I , K )*TEMP2( K, J ) 

PSI ( I , J ) = TEMP 
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